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Abstract 

We consider the classical three-dimensional motion in a potential 
which is the sum of n attracting or repelling Coulombic potentials. As- 
suming a non-collinear configuration of the n centres, we find a universal 
behaviour for all energies E above a positive threshold. 

Whereas for n = 1 there are no bounded orbits, and for n = 2 there 
is just one closed orbit, for n > 3 the bounded orbits form a Cantor set. 
We analyze the symbolic dynamics and estimate Hausdorff dimension and 
topological entropy of this hyperbolic set. 

Then we set up scattering theory, including symbolic dynamics of the 
scattering orbits and differential cross section estimates. 

The theory includes the n-centre problem of celestial mechanics, and 
prepares for a geometric understanding of a class of restricted n-body 
problems. 

To allow for applications in semiclassical molecular scattering, we in- 
clude an additional smooth (electronic) potential which is arbitrary except 
its Coulombic decay at infinity. Up to a (optimal) relative error of order 
1/E, all estimates are independent of that potential but only depend on 
the relative positions and strengths of the centres. 

Finally we show that different, non-universal, phenomena occur for 
collinear configurations. 
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1 Introduction 



The n-body problem of celestial mechanics is the problem of solving the Hamilton 
equation for the Hamiltonian function 



.2 z Z 

on the phase space T*M over the configuration space 

M := {(gi, . . . , qn) G {R^T \ % for i ^ k] . 

In celestial mechanics v = —\ and the Zi coincide with the positive masses 
mj > (in units where the gravitational constant equals one). However in an 
electrostatical context v = +1, the Zi are interpreted as charges and may be 
positive or negative. 

Whereas the one-body problem corresponds to free motion and the two-body 
problem was solved by Newton, the n > 3-body problem is analytically non- 
integrable. 

If one of bodies is much faster than the others then one may approximate 
its motion by considering the n-centre problem. There the Hamiltonian function 
H : T*M ^ M on the phase space T*M over the configuration space 

M ■.= R'\{s^,...,Sn} 

of that body is given by 

n „ 

H{p,q):=lp' + V{q) with V{q) := -J^ T^:^- U-l) 

1=1 1^ ^'1 

The 2-body problem thus reduces after separation of the centre of mass motion 
to the 1-centre problem, that is the Kepler problem. 

Moreover, the 2-centre problem is integrable and has been solved by Jacobi 
(see Appendix B). This solution is particularly relevant for calculating the motion 
of artificial satellites, since the gravitational field of the earth can be approximated 
by the one of two centres (analytically continued, since the earth is oblate and 
not prolate). This has been used by Vinti in [Vi]. See |GKM] for an application. 

However, like the n-body problem, the n-centre problem for n > 3 is ana- 
lytically non-integrable, see Bolotin [Boj . 
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The analogy between the Coulombic and the gravitational interaction is not 
perfect, since in electrostatics one has repulsive as well as attractive forces. This 
allows us to approximate molecules by static configurations of nuclei surrounded 
by an electronic cloud of opposite charge. Thus the electronic potential V of the 
molecule has the form 

V{q} = j2T^\+W{q), (1.2) 

si G being the position, and —Zi ^ the charge of the Zth nucleus, multiplied 
by the test charge. The smooth electronic potential 14^ : M'^ ^ M may partially 
shield the nucleonic charge. We model this by assuming the existence of a net 
asymptotic charge Z^o G M with 

V{q} = ^ + 0{\q\-'-') (g-oo). (1.3) 

Thus we consider the generalized n-centre problem (II. ip with these Coulombic 
potentials V. 

The understanding of the motion in an n-centre potential, n > 3, is very 
limited if one considers negative energies E. There one expects a complicated 
mixture of ergodic components and motion on KAM tori. 

However, by combining known techniques of celestial mechanics, we show 
in this article that the tiigti energy motion allows for a more or less complete 
qualitative and even quantitative description. 

Up to an error of relative order 1/E, the quantitative aspects treated here 
do not depend on the precise form of the smooth potential W but only on the 
charges Zi and positions si of the nuclei. 

The qualitative structures do not even depend on these data but only on the 
number n of nuclei. Thus the case n = 1 of an atom resembles the Rutherford 
case, whereas the n = 2-atomic molecule is similar to the two-centre problem 
solved by Jacobi. Here we are mainly interested in the case n > 3 where the 
dynamics is no longer analytically integrable. 

One motivation for this work is to establish the basis for a geometric un- 
derstanding of certain restricted (n + l)-body problems of celestial mechanics, 
where one is interested in the motion of a fast test body in the force field of n 
bodies, whose motion is assumed to be known. 

Such an analysis should be based on perturbation theory for the n-centre 
problem. In a joint work [DK| with T. Dierkes, we intend to show that indeed 
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the (2-dim.) n-centre problem is structurally stable in the following sense. For 
a sufficiently small local perturbation there exists a homeomorphism conjugating 
the two phase portraits and leaving the asymptotic initial and final directions of 
the scattering orbits invariant. This homeomorphism is unique up to changes in 
the flow direction. 



Survey of Results We now describe the techniques and results of this paper. 



The class of potentials V under consideration is introduced, together with 
other basic definitions, in Sect. [21 An important consequence of the fall-off of 
V is the existence of an interaction zoneXZ, a ball in configuration space which 
contains the points si and which, once left, cannot be reentered by an orbit, see 
( I2.17p . In Def. 12.41 we formulate the standing assumption that the centres are 



non-collinear, i.e. no three centres are on a line. 

If V contains a negative singularity {Zi > for some /), then the Hamiltonian 
flow generated by (11.11) is incomplete. In the Kepler problem it is well-known 
that the only sensible way to continue a collision orbit is to reflect it at the 
singularity in configuration space, for this is the limit behaviour for the Keplerian 
conic sections in the limit of vanishing angular momentum. 

But for our purposes we need to control smoothness and energy dependence 
of the resulting flow, for potentials of the form (II. 2p . Due to the singularity at 
si the usual comparison techniques for o.d.e. are not effective for proving such a 
result directly. 

Instead, near si we apply in Sect. [3] the so-called Kustaanheimo-Stiefel 
transformation. The Hopf map 



(z,uiz) 



with the Pauli matrices ai, extends the Hopf fibration — ^ S"^. A cognate map 
of the phase spaces (cotangent bundles) is known [KSllStS ] as the Kustaanheimo- 
Stiefel (KS) transformation. It relates the positive energy Keplerian dynamics 
with the one of a particle in an inverted harmonic potential. 

The KS transformation was used in celestial mechanics (see, e.g. Aarseth 
and Mikkola [AaMi] and articles by Aarseth and by Heggie cited therein), applied 
to spectral problems of Schrodinger operators (see, e.g. [HK5W] ). and shown in 
[GKJ to regularize the semiclassical dynamics. 
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The perturbation estimates of Sect. |4] work in the covering phase space of 
the KS transformation. 

Although the leading contribution to V near si is the Kepler potential 

Zi 



\q- si\' 



(1.4) 



the following argument indicates that we cannot simply approximate a collision 
orbit of (II. ip by the Kepler hyperbola for the potential (11.41) which has the same 
initial conditions {po,qo). 

The effective scattering region of (ll.4p is a ball of radius 0{1/E) around si, 
see Lemma 19.11 On the other hand, the time the particle needs to reach the 
singularity at si from its initial position go is of order 1/^/E. 

Due to the (bounded) difference 

Viq}- 



of forces exerted on the two particles, within this time the distance between the 
true orbit and the Kepler orbit can grow to something of the order 1/E. So by 
the above the Kepler orbit starting at {po,qo) may miss the effective scattering 
region. 

But this means that shortly after collision the distance between the two 
trajectories does not necessarily decrease as oo. 

So instead we effectively show that for every orbit with initial condition (po, go) 
there is a Kepler orbit, whose initial condition differs only by 0{l/E) (in a natural 
metric) from (poj'fo). snd which remains in a C(l/_E')-neighbourhood through 
the whole collision process. In Prop. 14.11 the corresponding statements about 
perturbations of the inverted harmonic oscillator are formulated. More precisely, 
it is shown that these perturbed solutions are C^-near to the ones of the harmonic 
oscillator. 

Moreover the true scattering process is approximated by Kepler scattering in 
the sense, so that we may use the Kepler solutions if we want to linearize the 
true flow (see Prop. 18. 6i) . 

Whereas the KS transformation is particularly suited to understand the topo- 
logical and geometrical structures, its disadvantages consist in its local nature, 
the introduction of a higher dimensional phase space and the reparameterization 
of time. The first problem, the locality near one singularity, should not be insur- 
mountable. In fact Helffer and Siedentop found in [HS] a generalization of the 
KS transformation to two centres. 
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Time reparameterization, however, is unwanted since we are not only inter- 
ested in the orbits but also in quantities like the time delay of scattering orbits. 
Thus we describe in Sect. [5] another regularization method. In Theorem 15.11 
we show that the incomplete Hamiltonian system 

{T*M,uJo,H) 

may be uniquely extended to a smooth complete Hamiltonian system 

iP,u;,H) 

whose phase space P D T*M is a smooth six-dimensional manifold, u \s a 
smooth symplectic two-form on P restricting to the canonical symplectic two- 
form = Yl^=i ^ (^Pi °" T*M, and H : P ^ R \s a smooth Hamiltonian 
function with H\rp^j^ = H. 

In fact we linearize the flow near each negative singularity and then add a 
copy of M X S*^ to T*M. Here M parametrizes the energy and S'^ the incoming 
(and outgoing) direction of the collision orbit. 

Thus we may henceforth work with the complete smooth flow 

<l>:MxP^P (tGM) 

generated by H. 

By Theorem 15.11 we may think of the Coulomb singularity as an artefact of 
the use of the inappropriate phase space coordinates {p,q) which, however, leads 
to a non-trivial topology of P. 

In Sect. [6] we introduce the M0ller transformation by comparing the flow 
$* on P with the flow $^ : Poo Poo generated by the Kepler Hamiltonian 
function 

Hooip,q) ■■= - -j^ 

on its phase space P^o, with Z^o G M defined by (11.30 . 
The M0ller transformations are then given by 

■= lim $-*oIdo$^, 

t— »itoo 

where Id canonically identifies the two phase spaces outside a region near the 
singularities. 
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The fi^ are continuous and, under mild conditions on V near infinity, smooth 
canonical transformations, see Theorem 16. 31 a nd 1^751 So the flow is asymptotically 
complete, see Corr. 16.41 

However, we are not primarily interested in these typical properties, but in 
the specific traits of the multi-scattering dynamics. 

It is a general fact that to a large extend the unbounded motion, i.e., scat- 
tering is determined by the bounded orbits of positive energy. We denote by 6^ 
the set of such bound states of energy E. 

To control these orbits, we combine the above perturbative results for the 
single scattering process with an analysis of the motion inside the interaction 
zone IZ, but away from the singularities (Sect. [71). This turns out to be a 
C^-perturbation of relative order 0{1/E) of the free motion. Using this result 
and the one of Sect.|8l we may approximate the true dynamics by a combination 
of free and of Keplerian motion. 

In Sect. [5] the perturbative results of Sect. |4]for the flow in KS space are 
used to obtain Prop. I8.6] which says that the single scattering process is C^-near 
to pure Kepler scattering, up to a relative error 0{1/E). 

Now in Sect. M it is shown that if a trajectory is not strongly reflected in 
uniformly bounded time by the singularities, it soon leaves the interaction zone 
(Prop.Ea). 

The bound states are then analyzed in Sect. [15] by a Poincare section 
technique. We erect Poincare surfaces which, in their projections to configuration 
space, sit between pairs of singularities and then define in (110.41) a Poincare map 

Ve. 

In Proposition 111.21 of Sect. [Til we estimate the linearization o^Ve, which, 
up to a relative error 0{1/E), only depends on the scattering angle at si and 
the charge Zi. 

This allows us to establish in Sect. [T21 Thm. 112.81 a symbolic dynamics 
for bE, E > Eth- The bounded orbits are described by their sequence of near- 
collisions, (which is well defined, since we assumed the configuration of centres 
to be non-collinear). 

Thus for n = 1 there are no bounded orbits, for n = 2 the set Be consists of 
one closed orbit, which is closed. For the case n > 3 of primary interest, Be is 
locally homeomorphic to the product of a Cantor set and an interval. 

All bounded orbits are hyperbolic. Be has measure zero w.r.t. Liouville mea- 
sure on the energy shell. The Morse index of a periodic trajectory equals its 
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number of near-collisions. 

In Sect. [T3] the fractal dimension of this set Be of bounded orbits is esti- 
mated. More precisely, we consider its HausdorfF dimension diiiiT^ and its upper 
box counting dimension dimg, since it is known that most other definitions lead 
to numerical values between these two. 

Theorem 113.51 says that for energies E above the threshold Ei^ and n > 3 
centres 

1 + 2d{E) ■ (1 - 0{{ElnE)-^)) < dimn{bE) < 

< dimeibE) < 1 + 2d{E) ■ (l + 0{{E\nE)-^)) , 



with the solution d{E) of a finite matrix eigenvalue problem ( 113. lOp . In particular 
they meet the rough estimate 

dim„(6,,) = 1 + + O ((lnE)-2) = dims(6^) 

In Sect. [T4] the topological entropy of the flow of energy E is estimated. 

Topological entropy is a quantity which, roughly speaking, measures the in- 
formation loss per time unit about the state of the system. Here for positive 
energy E the energy shell is non-compact, so that we use Bowen's definition 
of entropy. Prop. 114.51 states that 

^top($s) = ^top($]jr6^), 

i.e. that the source of information loss is the intricate structure of the set hE 
of bounded orbits, whereas the scattering orbits only have a sub-exponential 
dependence on initial conditions. 

So the estimates in the proof of Thm. 114.41 can be based on symbolic dy- 
namics. It states that for E large /itop(^ij) = for n = 1 or 2 centres, whereas 
for n > 3 

/^top($^) = /^t% ■ ■ (l + + 0{l/E)^ . 

Here the constants and Ch^^^ are determined by solving a finite matrix 
eigenvalue problem. 

Whereas the factor \/2E is of kinematical nature, in the simplest case of 
an equilateral triangle [n = 3) resp. tetrahedron [n = 4) of side length d the 
constant h'^^ equals \ia.(n — l)/d. 



9 



This divergence of topological entropy is not in contradiction with integrability 
(in the sense of independent C°° integrals of motion) of the hamiltonian dynamics 
above the energy threshold (compare with Bolsinov and Taimanov [BTJ). 

Sect. [15] is devoted to the classification of the scattering orbits. After 
excluding asymptotic directions in cones of aperture around the axes 

through two nuclei, and near the forward direction, one obtains such a universal 
classification. Thm. [T575] states that the orbits of given energy and asymptotic 
directions are enumerated by the succession of the nuclei they visit. In particular 
they form a Cantor set if n > 3. 

The differential cross section -^{E,9^ ,9^) of the scattering process is an- 
alyzed in Sect. [TSl Roughly speaking, this experimentally accessible quantity 
is the 'probability' that a particle of energy E and initial direction 9^ has final 
direction 9^ . For general potential scattering in one cannot even expect the 
cross section measure on S'^^^\9^ (defined in Def. 116.21) to be absolutely contin- 



uous w.r.t. Haar measure, see [Kn3]. Here, however, this is the case for energies 
E > EtYy, and ■^{E,9~ ,9^) is defined as the Radon-Nikodym derivative. 

Moreover (after excluding cones of .E-independent aperture) by Thm. 116.41 
it differs from the well-known Rutherford cross section in 

1 \ / ry \ d-1 



of the single Coulomb potential with charge Z = ^/Yl^=i o^'v by 
da ,„ ; , , /da 



d9+ 



iEj-,9^)=(^iE,9-,9^)) .(l + Oil/E)). 
\d9+ /ru 



So for these asymptotic data the intricate structure of the scattering orbits is 
not showing up in the cross section. In fact, Thm. 116.41 also states that the 
differential cross section is approximated by the Rutherford cross section outside 
the (much smaller) cones of aperture 0{1/\/E), though with less accuracy. 

In the final Sect. [T7]we show by counterexamples that many results of the 
paper do not generalize if we drop the assumption of non-collinearity. 

The first Appendix is devoted to a comparison between the two-dim. case 
treated in [KK] and the three-dim. situation of the present paper. Basically, the 
analysis of |KK] is the one of geodesic motion on smooth many-handled surfaces 
of negative curvature, whereas here we apply perturbation techniques around 
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infinite energy. 

In the second Appendix we describe the (known) solution of the purely Coulom- 
bic two-centre problem, and its bifurcation diagramme. 

Notation. Transposed vectors are used according to typographical, not mathe- 
matical needs. 

Acknowledgement. I thank the Max Planck Institute for Mathematics in the 
Sciences (Leipzig), where a part of this paper was written, for hospitality and 
support. 

2 Basic Definitions 

We consider the time evolution generated by a Hamiltonian function 



with n Coulombic singularities of the potential V situated at the points 



To control the asymptotic behaviour, we assume that V decomposes into the 
sum of a purely Coulombic potential and a short range potential. By this we 
mean the following: 

Definition 2.1 A smooth, real-valued function V on the configuration space 



H{p,q) + 



(2.1) 



(si 7^ 4 for ij^k). 



M := R% 



\{si,..., 



Sn} 



(2.2) 



is called Coulombic if 



1. There exist Zi 0, I e {1, . . . , n}, such that 





(2.3) 
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2. The potential vanishes at infinity, i.e. 

\imViq) = 0, (2.4) 

and there exist Z^o € K, called the asymptotic ctiarge, e G (0, 1] and 

> 2max(|si|, . . . , |s„|) 

such that for some Ci > 



VV{q)-Z^^ 



< m > i^min) (2.5) 

1^ 



and 



mm(|gi| , \q2\y^^ 

(2.6) 

Example 2.2 For the class of purely Coulombic potentials 

n 

the asymptotic charge = Ym=i ^i- Here V meets (12. 5p and (12.61) with e = 1 
(and Ci = ^nZ^^^). 

The question of dynamics in a Coulombic potential is called the n-centre problem 
of classical mechanics. 

Remarks 2.3 1) For Zi > and W = this corresponds to the idealization 
of a celestial body in the force field of n other bodies of masses Zi with fixed 
positions si. 

But our definition also covers the physical situation of (classical) scattering 
by the potential of partially ionized, neutral {Zoo = 0) or negatively charged 
quantum molecules. There Zoo does not coincide with J2^=i^i- Again, one ex- 
pects potentials with e = 1, due to exponential decay of the quantum mechanical 
eigenfunctions of the bound electrons, which in turn leads to an electronic charge 
distribution which decays exponentially (see Agmon |Ag| ). 
2) Note that V{q) and the asymptotic potential —Zoo/ \q\ appear symmetrically 
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in (12.51) . and that (12.41) and (12.61) are met by the asymptotic potential (the last 
one even with e = 1). We choose Ci large enough so that with the given con- 
stant e (12. Sp is met by both potentials. This will be used in Sect. [6]for simplifying 
the existence proof of the M0ller transformations. 

3) By a suitable translation of the origin, one could for n > 2 assume that -Rmin 
is, say, twice the maximal distance between the centres (but this assumption will 
not be used). 

For n > 2 we set 

dmin := inin d'^'' and rfmax := max d'^'' with d^'^ := \sk — Si\ . (2.8) 

In the case of a single atom [n = 1) we fix dmm by setting rfmm := 2i?min, say. 
A threshold of the energy is 

K.ax := sup (|v(g-) I qe M,VV{q) = } U {0}) . (2.9) 

One has < Vmax < oo, and in many cases Kiax = 0. 

We will generally assume that there are no more than two nuclei on one line: 

Definition 2.4 The configurations of singularities in 

ATC := { {s^, ...,Sn)e (M|)" \ ^ j ^ k t : {s, - s,) x {s, - 4) o} 
are called non-collinear (NC) configurations. 

In three dimensions this is a weaker assumption than that the nuclei are in general 
position, since the latter also means that there are no more than three nuclei in 
one plane. 

The space J\fC of configurations is a smooth connected manifold. 

Most statements in this article are shown for values of H above some thresh- 
old energy Et^. For given W (in particular for the purely Coulombic case W = 0) 
and charges Zi this is a function Eth ■ -NC ^ R of the positions (si, . . . , s„) of 
singularities. Et\i diverges near the set \ A/C of partly collinear configura- 

tions. This is not an artefact of our methods of proof. Instead, we will see in 
Sect. [T7]that several assertions proved in previous sections become wrong if one 
drops the NC condition. 
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One parameter measuring the degree of non-collinearity of a configuration of 
n > 3 nuclei is the minimal angle 

ttmin : -A/C ^ (0,7r/3] amin := mill a{i,j,k), (2.10) 

where for i ^ j ^ k 

a{z,j,k) G [0,7r), cos{a{z, j, k)) := {P'\s^'''), (2.11) 
with the directions 

._ ~ 
^ • I ^ —ill 



I 



is the angle between s,; and s^, seen from sj. For n = 2 we set a^i^ '■= tt/S. 
The Hamilton equations 

p = —'VV{q) , q=p 

lead to solutions t i-^ {p{t),q{t)) = $*(a;o) of the initial value problem with 
initial values xq = {po, go) which exist uniquely up to eventual collisions with the 
si. In Sect. [5] we extend uniquely to a smooth complete flow. 
The virial identity 

I {m,m) = - vm)) - im^^^vim)) (2.12) 

holds true for any trajectory t {p{t), Q{t)) with energy E := H{xq) (whenever 
q{t) 7^ si). Let us choose a function R^i^ : (0, oo) — * M of the energy, called the 
virial radius, with 



max 



(in9)l,^)<f and \{q,VV{q))\<El2 > i?.,.(E)). 

(2.13) 

As a consequence of part 2 of Def. l2.1l of Coulombic potentials, such a function 
exists. Without loss of generality (w.l.o.g.) we assume R^-^^ to be continuous, 
nonincreasing and constant for energies E > Et^ above some threshold. Property 
(I2.13P already implies that R^ij.{E) > 2Zoo/E. Technically we choose i?vir with 

R,,,{E)>ma.x{2R^^,C2/E) with C2 := 31{1 + l/e)R]^:^C, (2.14) 
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Then by (IZT^ and (IZBil 

I {m,m) > f > if \m\ > RvUe). (2.15) 

We conclude that a configuration space trajectory t i-^ q(t) of energy E leaving 
the ball of radius i?vir(-E) cannot reenter this ball in the future but must go to 
spatial infinity: 

Namely assume w.l.o.g. that (g(0),p(0)) > 0. By ( lZl5i) 

(P d 

—q\t) = 2^^mJit))>E (t>0) 

so that 

f{t) >qo+ lEt" (t > 0). (2.16) 

We shall mainly deal with energies E > Eth so that we may consider the E- 
independent interaction zone 

IZ := {qe M| I \q\ < i?vir(^th)}. (2.17) 

Sometimes we will use outside JZ the shorthand 

Kr:=V-Vo, with Vo,{q):=-^ (2.18) 

In addition to the condition (12.50 on VV^r this potential is of short range in the 
sense V^riq) = 0{\q\^^~'), too, as follows from (JZl) and (lZ5j) . 

Due to collisions with the nuclei situated at si, the flow on the phase space 
T*M of the particle is incomplete. There are several ways to regularize the 
collision orbits which are all essentially equivalent. In Section [3]we now introduce 
the Kustaanheimo-Stiefel (KS) regularization method, and use it in Sect. |4] for 
comparison estimates which control the deviation from the Keplerian motion 
near a singularity. Later, in Sect. [5l the flow is then regularized without a time 
change. 

3 The Kustaanheimo-Stiefel Transformation 

in |KS] Kustaanheimo and Stiefel related the Kepler motion in three spatial 
dimensions to the motion of a resonant harmonic oscillator in four dimensions, 
thus linearizing the dynamics, see the book jStS] by Stiefel and Scheifele. 
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Although the authors had used spinor theory in order to derive their results, 
the emphasis of their article was more on the application to perturbation theory 
then on the geometry of the problem. In [Ku] Kummer presented this aspect of 
the KS-transformation and related it to the approach [Mo] by Moser. 

Our presentation will be based on the quaternion algebra over M 



Z2 Zl 



Zl, Z2 e 



with matrix multiplication (see, e.g., [KR]). We use the basis 

= (11, ZO-i, ZCr2, ^0-3) 

of H, the ai being the Pauli matrices. The direct sum decomposition 

H = M ■ 11 © ImH 

with 

ImH := {Z eM \ = \ -l with A < 0} 

= Spanjj(Ji, 12,13) 

into real and imaginary space is orthogonal w.r.t. the inner product 

HxH^M , (X,F) := |tr(XF*), 

1 

X X* := X* being the conjugation. The norm |X| := (X, X) 2 is multiplica- 
tive: 

|XF| = |x||r| {x,Yem). 

The vector product x : ImH x ImH ImH is given by 

X X F := |(Xr - FX), 

and we have 

XF = - (X,r) 1 + X X F (X,yGlmH). 
We consider the Hopf map 

TTo : H ^ ImH, 7ro{Z) := Z*hZ = % ( ^^^^ -2^i22_ \ ^3 -^^ 



-2ZiZ2 Z2Z2 — ZiZi 
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which is a surjection — > whose preimages are the orbits of the isometric 
group action 

ao-.S^^ Aut(e), ao{<^){Z) := exp{l3if)Z. 
This action is free on H := ]HI\{0}. We have the canonical symplectic one-forms 

e := |tr(P*dg) = 3f?(pidgi + P2dq2) ((P, Q) G T*U) 

and 

^ := ltr{P*dQ) ((P, g) e T*ImM) 

on the cotangent bundles, and denote by d :— d\j.^^, resp. 9 :— 0\j.^j^ their 
restrictions. 

The restricted Hopf map ttq := ttq \^ is then related to the KS-transformation 

Tf : T*m -> T*ImH 

^(P^Q) = {^{Q*hP + P*hQ),Q*hQ) (3.2) 

of the cotangent bundles. 

We consider the quadric surface 

S :=I-\Q) C T*H 

for the bilinear form 

J : T*H ^ M, X(P, Q) := |tr(Q*/3P) = + ^2^2) 

and its restriction S -.^ S (1 T*'E.. Then 

0r5 = ^*^^r5, (3-3) 

since on S 

Q*hP = P*hQ e ImH. (3.4) 
T, d and tt are all invariant w.r.t. the group action 

a:S^^ Aut(T*H), a{^){P,Q) := (exp(/3(^)P, exp(/3(^)g), (3.5) 
and T is a Hamiltonian function generating that time-</7-flow. 
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Lemma 3.1 Restricted to the quadric surface S C T*M, the lift of a Hamilto- 
nian function H{p, q) = ^Ip]'^ — Z / \q\ + W{q} equals 

H o 7r(P, Q) = \Q\-' (||P|2 -Z)+ WiQ*hQ) ((P, Q)eS). (3.6) 

For the angular momentum L{p, q) = qx p 

L o 7r(P, Q) = \{Q*P - P*Q) ((P, Q)eS), (3.7) 

whereas the Runge-Lenz vector F(p, q) = p x L — Z-^ transforms into 

F o 7r(P, Q) = -lP*hP + {H o 7r(P, Q) - W{Q*hQ)) ■ Q*hQ ((P, Q) e S). 

(3.8) 

Proof, (see also |Ku]). 

By (13. 4p the restricted KS-transformation 7i\g maps |QP = \Q*l3Q\ onto \q\ 
and \2Q\-' ■ \P\' = \{2\Q\^)-'Q*hP\^ onto \p\\ implying B^. 
Since g,p G ImH, 

L = q X p = i(gp -pg). 
The KS-transformation vr gives, using (13. 4p 

f {Q*hP + P*hQ)Q*hQ - Q*hQ{Q*hP + P'-Jsg) 

= m 

P*hQQ*hQ - Q*hQQ*hP , ^ 

4|g|2 -iW^ 

This shows (13.71) . To prove (13.81) . we use the identity 

p X {q X p) = \{pqp — qpp) = \{pqp + q\p\'^) {p, q € ImH) 
which implies that 

Fip, q) = \pqp + q ■ iH{p, q) - W{q)). 
The term pqp lifts to 

{2\Q\')-\P*hQ){Q*hQ){Q*hP) = -\P*hP, 
finishing the proof. □ 
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4 Application of the KS-Transformation 



In a configuration space ball 

Bi{r) := {gGM|| \q- Si\ < r} 

of radius r := Cq < Idmin around g = we want to regularize the motion 
generated by the restricted Hamiltonian function 

H{p, q) = ibT - + Wi{q), ((p, q) G T*mc,) \ {0})) , 

Wi{q}:=J2—^ + W{q} (4.1) 

being a smooth function on Bi{cq) C M'^ = ImH. The radius Cq G (0, ^rfmin) is 
chosen so that 



\Wm\ < i (/ = 1, . . . , n , |g - .11 < Cq). (4.2) 

1^ I 



For simplicity of notation we assume si = 0. For regular values E of the en- 
ergy H the orbits on H^^{E) coincide with the ones of the zero surface of the 
Hamiltonian function 

\q\-iHip,q)-E), (4.3) 

since they are already determined by the form of the submanifold H~^{E). 
By (13. 6p the lift of (14. 3 p with the KS-transformation equals 

He{P, Q) := I (P, P) + (g, Q) i~E + Wi{Q*IsQ)) - Z^. (4.4) 

The orbits of its Hamiltonian vector field on S* fl (H£;)~^(0) project to the ones 
of H on H^^(E), since the symplectic one-form 9 on the phase space region 
T*{Bi{cq) \ {0}) lifts according to (lO) . 

For Wi = the Hamiltonian function H^; is the one of a four-dimensional 
harmonic oscillator (with negative potential for E > 0). In the general case the 
additional potential Q t-^ Wi{Q* I3Q) will be a small perturbation if E is large, 
at least for \Q\ < ^JcTq. 

From now on we assume initial conditions to lie in S* C T*HI. 

Instead of considering the Hamilton equations 

= 2EQ + 2R{Q) 
TsQ = \P 
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of (14. 4p with perturbation 

R{Q) := -QWi{Q*hQ) - l\Q\^VWi{Q*hQ), (4.5) 
we want the parameter E to appear only in the perturbative term. Thus we set 

X := ( ^ ) with P := P/VSE 



Q 



and use the time variable r := ■ s. Then 

with R(X) - 

The initial conditions Xq := meet the energy constraint 

\Po\' - \Qo\\l - Wi{QlhQo)/E) = |. (4.7) 



So with 
and 

the region 



|p|2_|g|2 ^ ^ 



T>i := |(P,g) G M'* X \Q\ < Cq,S{P,Q) > E,^ 



of phase space points meeting (14.70 for some E > E^h is bounded. 
It is natural to compare the solution X{t) of ( 14. 6 p with the solution 

^(-)=f"tMl! "ImUVo (4-9) 
^ ' \ smn(r)ll cosn(r)ll J ^ ' 



of the linear initial value problem 

dr V 11 



^-^={1 l]-^ ■ nO) = Xo. (4.10) 
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We denote this linear flow (corresponding to Wi = 0), restricted to the invariant 
phase space domain 



th 



'bL,i := { (P, Q) G M'* X \P\ ^\Q\, Sl{P, Q) > \E, 
with 

8dP.Q):=Zi/{\P\'-\Q\') (4.11) 

by 

^L,i ■■ M X 'bL,i ^ X>L,Z, (4.12) 

and set 

VL,i:={{P,Q)ei>L,i\\Q\<CQ}. 

Our condition (14.20 on the radius Cq implies that 

T>L,i 3 'Di- 

Let 

:X>; -^MU{±oo}, T;^(Xo) :=±inf{t >0 | |Q(±t, Xq) | > Cq}. 

be the exit times from X>/, 

^r.Vi^ T>i, Ui := {(t,Xo) eRxVilte [Tr(Xo),T+(Xo)]} 

the maximally extended KS flow on X>/, and 

^fi^ :T>i^dT>i , X^^i{Tf{X),X) 

the map to the exit points. 

Basically we are interested in the Poincare map 

f */+(x) ,(p,q) < 

: dVi ^ dVi, X = {P,Q)^ { > / (4.13) 

[ *r(x) ,(p,g)>o 

that permutes incoming and outgoing data, but up to now we do not even know 
whether this is defined everywhere. 

Therefore we compare with the linear flow (I4.12p and thus introduce in anal- 
ogy its exit times T^^ : T>l; — > M. They are finite and smooth, and we set 

^%iX) :=$^/Tf;(X),X). 
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The analog 

of (14.131) for the linear flow is a involutive diffeomorphism: 



with (4.15) 



QJ \ Q-uP J ■ |p|2 + |g|2 

since |n| < 1 on T>i ;. 

Ideally, one would like to prove that (14.130 is approximated by the map ( 14.141) 



so that the C^-norm of o - Id is of order 0{l/£). This would be 

true if the trajectories of ^i, resp. would spend a uniformly bounded time 
inside X>i, resp. T>l,i- But this is not the case, and in order to keep the error 
terms small we will compare the two flows for initial conditions on the pericentric 
hypersurface 

Til := {{P,Q)eVi I {P,Q) = o'j cnL,i := {{P,Q)e'bL,i I (p,g) = o}. 

Note that by transversa lity of the linear flow to that hypersurface the pericentric 
time 

Tly'i>L,i^^ with *°/X):=$^,KT°/X),X)G?iz.,^ (4.16) 
is uniquely defined and smooth. 

Proposition 4.1 For E^h large there is a unique pericentric time Tj^ : 2?^ — > M 
witti 

:= $KT?(X),X) G Hi, (X e Vi). 
Ttie functions < Tj* < T;+ : X>/ ^ M are smootti, and 

|T±(Xo) - T^/Xo)| = Oil/SiXo)) (Xo G Tii). (4.17) 

Ttie exit times T^^ of ttie linear flow are estimated by 

exp{±2Tl,iXo)) = -pf—SiXo) + 0{S%Xo)) (Xo G 7^,), (4.18) 
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e(Xo) := y 1 + ^-^l-^^^")!^ being the eccentricity of the corresponding Kepler 
hyperbola. The diffeomorphism 



o*['(X) ,(P,Q)>0 



onto its image which conjugates the maps ("^i = (S;) ^ o ^j^ i o S^j is C^-near 
to the identity in the sense that 

\Ei{X) -X\= 0{1/S{X)), (X G d-Di), (4.19) 

and the solution ^^^i of the linear problem is C^-near to in the sense 

\<fi{X)-<fL,ioEi{X)\ = 0{1/S{X)) (4.20) 
\\D<fi{X)-D<fL,ioSi{X)\\ = 0{l/e{X)). (4.21) 

Remarks 4.2 1) The energy-independent estimate (I4.2ip may seem to be poor 
but is in fact optimal in its energy dependence since, relative to the optimal 
estimates 

= 0{£iX)/eiX)) = \\D<bL,iiX)l (4.22) 

it is of order 0{1/E) (The r.h.s. of (14.221) is obtained by inserting the time bound 
(I4.18P into the linearization of PTOi) ). 

2) For pericentric initial data Xq e Tii the total time spent inside the ball equals 
T^;(Xo) - Tli{Xo) = ±2Tf i{Xo). Estimate (l438l) for that time is presented 
in a form needed to evaluate the term D^l,i ° 3;(X) in ( 14.211) . 

In polar coordinates (r, 99) the Kepler hyperbola has the parametric form (see 
e.g. [Th], Chapter 4.2) 



\Z\ecos{ip — (fo) + Z 

The denominator has the zeros (p^, and A(p := — (p^ is the angle under 
which the hyperbola is seen from the origin. Thus 

cos{^A(p) = , 

so that Aip e (tt, 2n] for Z > and Aip e [0, tt) for Z <0. 
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On the other hand, the total change in direction of the velocity vector 
equals 

Alp = sign(Z) ■ {Aif - tt) = 2arcsin(l/e(Xo)) G (0,7r]. 
Thus (14.181) can be rewritten as 

exp (±2T£,(Xo)) = 4c.sm(jAV^) ^^^^^ ^ OiS\Xo)) (Xo G Ti,). 

(4.23) 

This equation will be useful for the study of orbit instability, since exp (t) equals 
the expansion of the unstable manifold of (14. 9 ft after time t. 

Proof. We set {P{t),Q{t)) := $/(t,Xo) for Xq = (Po,Qo) e X>z. For ^th 
large the squared distance t |Q(t)P is a strictly convex function of time, since 



2,JQiW = -Am.Qit)} (4.24) 



= |P(t)p + |Q(t)p - (i?(g(t)), Q{t)) /S{Xo) 

> \P{t)\' + m)\'.{l~C^cJS{Xo)) 

> i(|p(t)p + ig(t)p) >o 

with Lipschitz constant 

A:=sup| '^^^^|~^^|^^^' |Qi|,|Q2| <CQ,Qi^Q2}, (4.25) 

(one notes from inspection of (14.51) that R(0) = 0). 

We can bound (I4.24p more precisely from below by using the inequality 

|Xp ^ \P\' + \Q\' > l\\P\' - \Q\\1 - W,{Q*hQ)/S{X))\ = 

(4.26) 

which follows from (14.71) and is valid for E^h > maXg'gB;(cq) 

Thus Xf and T° are uniquely defined finite functions. By transversality of 
Til w.r.t. the flow the pericentric time Tj* is smooth. The hypersurface d'Di is 
transversal to the flow, too, except at dT>i fl 7ii. Thus it is only there that we 
have to control smoothness of the exit times T^. 

The maps 

Vi^RxHi, X^(TO(X),*°(X)) 
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and 

RxHi^W^x M^ {t, X) ^ $L,i(-^> X) 

are diffeomorphisms onto their images, since X>i ; does not contain the (single) 
stationary point of the linear flow Hence the composition 

pr.Vi^ V^^u Pi{X) := $z.,K-T°(X), (4.27) 

of these diffeos is a diffeomorphism onto its image. 

Thus in order to compare the flows and it suffices to compare the 

trajectories 

X{t) ■= {P{t),Q{t)) := Mt,Xo) and Y{t) := (Pi(t), Qz.(t)) := *L/t,Xo) 

for pericentric initial conditions Xq = {Po,Qo) G Tii. 

We partition Hi into the regions Hi"^ := {(-P, Q) E Hi \ \Q\ < ^cq} and 

Hi> ■.= {{P,Q)eHi\'^CQ<\Q\<CQ}. 

The flow through Hi^ is uniformly transversal to dl^L^i, whereas the exit times 
are uniformly bounded on Hi^ . 

1) On Hi'^ the second equation in (I4.24p together with (14.71) yields for 



g,{t) := \Qit)\' '2^3iit) - 2g,{t) = h,{t)/8{X,) (4.28) 

with 

h,{t) := Zi - {R{Q{t)), Qm - \Q{t)\'Wi{Q*{t)hQ{t)) {\t\ < T+(Xo)), 
^i(O) = IQoP and ^^(O) = so that 




9iit) = \Qor + / / (4^71 (m) + 2hiiu)/£iXo))duds. 



JQ 



As Cs := sup|Q|<,^ \Zi - {R{Q), Q)) - \Q\''Wi{Q* hQ)\ < oo, for E,^ large and 
times t between and T^(Xo) the integrand is bounded below by 3|(5oP and 
above by 5cq so that 
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As g,{Tf{Xo)) = cl, 

|y/i-|golV4^±Tf(Xo) < 2^1 - |golV4 (^o = (^o,go) e ni>) 

(4.29) 

and the same estimate holds for Est. (14.291) implies in particular the uniform 
bound \Tf\, |T^ J < 2 on Tii^ . 

The difference between these times is much smaller: 

\Tf{X,)-T%{Xo)\ = O (y/l-|go|V4/^™) (^0 G ni>). (4.30) 

Namely setting g2{t) := \Qit)\^ - |gL(t)P, so that ^2(0) = ^2(0) = and, by 

(E2H]), 

we get 

g2{t)= [ [ {4g2{u) + 2h2{u)/S{Xo))duds 
Jo Jo 

or 

b(t)|<^sinh2(t) 

for C4 := C3 + Zmax- 

As 4 = |gz.(T^/Xo))p = sinh2(T±/Xo))|goP and 

^72(Tf (Xo)) = 4 - |gi(Tf (Xo))P = 4 - sinh2(T±(Xo))|go 



2^^,f72(t)-2(72(t) = 7^ with /i2(t) := /ii(t) - ^i, 



smh(Tf (Xo)) ■ ^/l - < sinh(T^/Xo)) 



±fv \\ /1 , 



<sinh(Tf(Xo))- Wl 



igoP^(Xo)' 



In view of ( KTQh this gives ( l43Qll . 
In turn (I4.30p implies 

|$i,KTt(Xo),Xo) -*L,KTf(^o),Xo)| + 

||D$z.,KTt(Xo),Xo)-Z)$L,z(T±(Xo),Xo)|| 

= O (y/l - |golV4/^(Xo)) . (4.31) 
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A comparison between the initial value problems (14.60 and (14.101) using a 
Gronwall estimate on the uniformly bounded time interval \t\ < 2 yields 

l^jiXo) - ^UXo)\ + \\D^liXo) - D^l,iXo)\\ < (4.32) 

(XoG?i,>,tG [Tr(Xo),T+(Xo)]). 

Setting t := Xf (Xq) in (ICTI) . the triangle inequality and (14.311 ) leads to 

|*,±(Xo) - *t(Xo)| + ||/^*^^(Xo) - /^*t^(Xo)|| = 
|$KTf(Xo),Xo)-$i/T^/Xo),Xo)| + 

||D$KTf (^o),Xo) - D$^,KT^,(Xo),Xo)|| = o (y/i - |golV4/^) • 

This does not only prove the estimate KW} on U C 9X>i, 

but also shows that S is continuously differentiable at the submanifold 

Til n = ^i+{ni>) n ^r{ni>) 

of phase space points where the flows and ^^^i are tangential to dT>i (for X 
in this set H(X) = X and DH(X) = 11). 

2) • On T-Ci"^ we begin with a rough estimate. For time 

r = r(Xo) :=ln(3cQ/|Xo|) (4.33) 



\QL{±r)\ = I cosh(r)Qo + sinh(r)Po| = J smh'{T)\Po\^ + cosh'(r)|QoP 



> ^i|Xo|2 . (cosh(2r) - 1) > ^i|Xo|2 • (iexp(2r) - 1) 

V7 



A^l 2 



\W > (4.34) 



2 



since by (14. 7 p 

_ .„ ,. ~ . _ . Z, - \OnPW,(0%TM) 

< c„. 



\Xo? = 2|QoP + (iPoP - IQcH < k + - WWjmM) 



2^1 ^ ^ — ^1 

Therefore any trajectory of the linear i\o\N $l ^ with these initial conditions leaves 
the region T>i before time r (and enters it after time — r). 
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• We now do perturbation theory around these linear solutions, setting 

Zit) ■.= X{t)-Y{t). 
Z(t) meets the integral equation 

m= f iil l\z{s)+E-'R{X{s))\ds 



_ _ ]1 0^ 
with E := £(Xq). Gronwall's inequality says that 

f{t)<Aexp(^J^ g{s)ds^ (4.35) 

if f{t) < A + J^g{s)f{s)ds. Applied to f{t) := \Z{t)\, we get for < t < 
T/-(Xo) 

fit) < (^f{s) + E-'\R{X{s))-R{Y{s))+R{Y{s))\)ds (4.36) 

< E-' j\R{Y{s))\ds + 1^ (^l + ^y{s)ds 

< E-^C, j\QL{s)\ds + j'^[l + ^^f{s)ds 
with the Lipschitz constant Ci from (14.251) . But 

/ \QL{s)\ds = I I sinh(s)Po + cosh(s)(5oM'5 
Jo Jo 

< \Xo\ [ e'ds < \Xo\e\ 
Jo 

Thus the constants in ( 14.350 can be chosen as 



and we obtain for \t\ < r(Xo) with X{t) G T>i 



\Z{t)\ =/(t) <Ci\Xo\e 



l,| exp((l+A/E)-r(Xo)) 



E 

1,1 1 
1*1 



|Xo|el*l ■ 0{E~2) = 3cq ■ 0{E-2) (4.37) 
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for EtY, large since by ( 143^ and i^TM 

r(Xo) = ln(3cQ/|Xo|) < In (Sv^cq/E^^) = \n{cE\). 
In particular we conclude from (I4.34p that 

0<±T±(Xo) <r(Xo), 

since otherwise by (14.37^ 

|Q(±r)| > \Ql{±t)\ - \Z{±r)\ > ^cq - 3cq ■ O (^E"!^ > cq. 

We use (I4.37P as an input for a refined estimate which will imply \Z(t)\ = 
0{l/E). To that end we note that by (143711 

m)\ < \QL{t)\ + \Z{t)\ < el*l|Xo| + \Z{t)\ < 2el*l|Xo| < r(Xo)). 

We write ^ 

Z{t) = E-' [ exp (( 2 S ) (t - .)) R{X{s)) ds. (4.38) 

W.l.o.g. we may assume that Wi{si) = 0, since otherwise we may shift the energy 
E by that constant, producing an error term of relative order 0{1/E). Then 
instead of the Lipschitz estimate (14.250 for R we use 

R{Q)<C2\Q\' {\Q\<cq). 
Inserting this into ( l438i) we get for < t < r(Xo), using ( l433l ). 

\Z{t)\ <^l\xp{t-s)e'^\X,\'ds<^-^{e'-l)\X,\ < ^-^\QL{t)\, 
since 



\QLit)\ = ^smh^(t)|PoP + cosh^(t)|QoP >sinh(t)|Xo|. 

A similar estimate holds for > t > —t{Xq). 
But \Q{t)\ > |Qi(t)| - \Z{t)\, so that 

I^WI< (^-i) '\Qit)\, 
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showing that the diffeomorphism pi of (14.271) onto its image is C°-near to the 
identity in the sense 

= 0{1/£{X)) {X = (P, Q) e Vi). (4.39) 

If we assume X E dT^i so that \Q\ = cq, then (14.390 shows that 

\Ql\-cq = 0{1/S{X)) for {Pl,Ql):=Pi{X). (4.40) 

Since X = ^'^^(Xo) with Xq G Tii^, and dT^i^i is uniformly transversal to the 
flow <^L^i through 'H{^ we obtain from (14.400 

|Tf (Xo) - T±/Xo)| = 0{l/S{X,)) (Xo G -nn- (4.41) 

In turn, this and ( l439l ) imply ( 09l ). 

• We also get ( 14.170 from the time estimates ( 14.300 and ( 14.410 in the two regions. 

• Estimate ( 14.180 for the exit time of the linear flow is derived as follows. By 
( 14.110 the 'linear' energy parameter equals 



lAP-igoP 

The (lifted) angular momentum (see (13. 7p of the pericentric initial data Xq = 
(Po,Qo) equals 

|L(Xo)| = \\QlP,-P*Qo\ = i|Im(QSPo)| = ||Po| |Qo| = V^£lpQ\Po\ \Qo\, 
since {Po,Qo) = 0, see ( 13.71) . Thus 
\Zi\-e{Xo) 



£l{Xo) 



(Zi/SLiXoW + 4|PoP|goP = \Po\' + IQo 



On the other hand the exit times T^^(Xo) and T^,(Xo) = -T^;(Xo) are 
implicitly given by the equation 

|gi(T+/Xo),Xo)| =CQ 
with QL{t,Xo) = sinh(t)Po + cosh(t)(5o, whence 

exp(2T+/Xo)) = I 5 i/'^i^ I, + OiSUXo)) 

I-mdI + Ivol 



4c,^i(Xo) 



\Zi\e{Xo) 
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+ Oi£liXo)). 



This proves ( KTEh . 

• The C°-estimate (14^ follows immediately from (l4l^ . 

• In order to obtain the C^-estimate (14.211) . one considers 

DZ(t) solves the integral equation 



DZ{t) 



11 

1 



DZ{s) + E'^DR{X{s))DX{s) ds 



W.l.o.g. we may again assume that Wi{si) = 0, since otherwise we may shift the 
energy E by that constant, producing an error term of relative order 0{1/E). 
Then by inspection of ( 14.50 one notes that the matrix DR((P, Q)) = at Q = 0. 
We may thus estimate 



||M((p,g))|| </:3|g| ((P,g)GX>, 



(4.42) 



where £3 > is the Lipschitz constant of DR for \Q\ < cq. 

From (14.42ft we obtain (14.211) by a Gronwall estimate similar to (I4.37p . which 
we apply to 



DZ{t) 



1 ' 



3) + E-'DR{X{s)) DZ{s) + E-'DR{X{s))DY{s 



ds. 



witht<T+(Xo),/(t):= \\DZit)\\, 

g{s) := 1 + CsCq/E > ||( « g ) + E-'DRiXis) 

and 

A :-- 



T+(Xo) 



DR{X{s))DY{s)ds 



To obtain the last estimate, the time bounds (14.171) and (14.181) are inserted. 
Thus for t < T+(Xo) 



f{t) < Aexp / g{s)ds )=0( E-2e{Xo)-2 • O E2e{Xoy2 



□ 
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5 Regularization by Phase Space Extension 



In Sect. |4] the motion in configuration space ]R| near a singularity at si was 
regularized using the KS transform. This will enable us in Sect. [8] to compare 
that motion with the motion in the Kepler potential q i— —Zi/\q — si\. 

However, we would like to apply these local estimates to a complete Hamilto- 
nian flow on a phase space which arises by a completion of T*M. Therefore in the 
case of attracting singularities {Zi > 0) we now employ a different regularization. 

To preserve continuity of the motion with respect to the initial conditions, a 
particle colliding with a nucleus at si E M| must be reflected backwards. Then 
we parametrize the state of the colliding particle by its energy and by its incoming 
(or outgoing) direction. That is, we complete phase space by adjoining manifolds 
M X S"^, one for each attracting singularity. 

Note that for all energies E the energy surfaces H^^{E) could be completed 
topologically using only the Kustaanheimo-Stiefel construction (by taking the 
quotient of the quadric surface X^^(O) w.r.t. the circle action a defined in (13.50 ). 
This, however would lead to a time change which is unwanted here. 

Theorem 5.1 There exists a unique smootti extension {P,lu,H) of the Ha- 
miltonian system {T*M,uj,H), where the phase space P is a smooth six- 
dimensional manifold with 

P := T*M U y (M X S^) 

l<l<n 

Zi>0 

as a set, u is a smooth symplectic two- form on P with 

3 

uj = a; := ^ cig. A dpi, 

and H : P —>■ M. is a smooth Hamiltonian function with H\j,,£j = H. 
The smooth Hamiltonian flow 

^■.RxP^P (5.1) 

generated by H is complete (and we often write $*(x) instead of ^{t, x) ). 
For all energies E > V^ax (defined in p.9\] ) the energy shell 

:= {x G P I H{x) = E} (5.2) 

is a smooth, five-dimensional manifold. 
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Proof. It is clear that the particle cannot escape to spatial infinity in finite time. 
If < there is no need here to regularize, since then by (14.20 the minimal 
distance is bounded by |g — s;] > \Zi\/{2E) if > resp. \q — si\ > Cg if < 0. 

For the remaining case of an attracting singularity [Zi > 0) we linearize the 
motion near collision by using as phase space coordinates the angular momentum 
components, the direction of the Runge-Lenz vector Fi, energy and the time 
passed since the pericentre of the orbit. The first five of these six functions are 
constant on the Kepler orbit. Then we add the collision manifold of phase space 
points with time and angular momentum both equal to zero. This manifold is 
then parametrized by energy and by Fi/\Fi\ and is thus diffeomorphic to M x S'^. 

In [KK] the case of d = 2 dimensions is treated, the calculations being more 
detailed than here. Comparing with the Delaunay coordinates (see, e.g., (AM], 
Chapter 9.3), instead of the semi-major axis we use energy, and instead of the 
mean anomaly (which is only defined for negative energies) we use time. 

If the potential V is not centrally symmetric around si, then at collision some 
of the former constants of motion cease to be smooth functions of time. To 
remedy this is we redefine them by using their value at the pericentre of the 
orbit. 

1) More specifically, we introduce adapted coordinates to regularize the flow in 
the phase space neighbourhood ^7f, < e < Cq, of the /th nucleus, with 



On U[, the Hamiltonian function has the form H{p, q) = Hi{p, q) + Wi{q} with 

Hi{P,q):=y'--^ (5.4) 



and the smooth additional potential Wi on Bi{cg) 

\n- 



One basic estimate on f/f , valid for e small, is 

at |g - si 



{{q-si)-p) > l^-^-{q-si)-VWi{q)>0. (5.5) 
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Every collision orbit with si enters JJ[, as 
as q approaches si. 

2) We first treat the Keplerian case Wi = 0. The angular momentum L/ : f/f — > 
relative to the position of the /th nucleus equals 

Li{p, q) ■■= {q - si) X p. (5.6) 

Let Ti : U[ ^ M. be the time elapsed since the closest encounter of the 
Kepler solution with the nucleus. By (15. SD there is only one such pericentre of 
the orbit, with distance rmm- Ti is given by 

Ti{p, q) ■■= I ■ sign((g - si) ■ p) (5.7) 

Jr^Um ^2r^Hi{p, q) + 2Zir - Lf{p, q) 

with 

r -Zi + ^Zf+2Hi (p,q)Lf {p,g) fj /n. 

r^Mq) ■■= { . ^^M'f) '"^"'^^ . (5.8) 

I mp,q)/2Zi ,Hi = 

Ti is a smooth function, which can be seen by explicit evaluation of the integrals: 



r 



V2r^E + 2Zr - L"^ 



dr = (5.9) 



Z L'^ Z 

1 + 



In (et + lZ+ ^E{r^E + Zr - |L2)^ 



y/2E\ rE 2r^E {2Efl^ 

for > and Z > (see Thirring [Th], for more information). 

The Runge-Lenz vector F; : IJf relative to the /th centre is given by 

Fiip,q) ■.= pxLi{p,q)-Zi-^^. (5.10) 

On its domain f/f of definition Fi is non-zero: \Fi\^ = 2\Li\^Hi + Zf > Zf/A. 
Thus we may define the pericentral direction Fi : Uf — * S"^ by Fi := Fi/\Fi\. 
The angular momentum vector L/ is perpendicular to that direction: Lr Fi = 0. 
The map 

y-.U^^ T*(R X S') \ 0, (p, q) ^ {fu U- Hu F,) 
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is a diffeomorphism onto its image, (0 denoting the zero section of the cotangent 
bundle T*(R x S^)). 

3) The Poisson brackets between the above variables are given by 

[fu Hi} = 1, , {Li, Hi} =6. {Fi, Hi} = 6. {Fi, fi} = and {Li, fi} = 0, 

so that in particular the angular momentum and the pericentral direction are 
constants of the Kepler flow. The components of angular momentum and the 
asymptotic direction have the Poisson brackets 

{{Fi),, {Fi)j} = 0, {(Li),, iLi)j} = e,jk{Li)k, {{U)^, {Fi)j} = ei,k{Fi)k, 

using the Poisson brackets {(F/)^, {Fi)j} = -2Hieijk{Fi)k. 

4) Because of the above Poisson brackets with Hi, by introducing the above 
coordinates, we obtain a chart in [/f which explicitly linearizes the (incomplete) 
Kepler flow. 

The motion is then regularized in the following way. One defines a completion 
of JJf by setting f/f := ?7f U (M x S"^) as a set, and one introduces a topology 
on f/f by extending the map y = {Ti, Lf, Hi, Fi) to 

y := {Ti, Li; Hi, F,) : f/f ^ T*(M x S^) 

by mapping {h,f) G M x S'^ onto the point (0,0; /i,f) of the zero section. By 
that procedure we obtain the topological manifold P and, by taking limits, we 
extend the Hamiltonian H to a continuous function H : P ^ M. The topology 
of P is thus determined by the purely Coulombic local Hamiltonians Hi, and, by 
taking limits for the collision orbits, we are able to extend the flow generated by 
to a complete continuous flow $* on P. 

Moreover, the calculation of all the Poisson brackets shows that we may 
continuously extend the symplectic form (X]i=i ^ ^Pi) lu^: ^""^ obtain 

a nondegenerate two-form, which is smooth in the new coordinates. 

5) To generalize the construction to the case of the flow generated by H which 
is of the local form Hi + Wi, we define similar canonical coordinates in U[ which 
linearize the //-flow $*. 

By (15. 5p for £ > small enough, the orbits are transversal to the hypersurface 
Sf. We extend that hypersurface to the topological submanifold 

5f := U (M X 5^) c f/f . (5.11) 
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If Wi = 0, we define the differential structure near the /th collision manifold 
M X S*^ by pull-back with 3^. In general we define a map 

y := {fi, Li- Hi, Fi) : f/f ^ T*(M x S^) 

by letting Hi := H\jje be the energy and 

fi{^\x)):=t (xeSf) (5.12) 



the time passed since the passage of the pericentre. Note that (I5.12p defines T; 
everywhere on U[ , since every orbit in U[ passes Sf exactly once. 



Li{<l>\x)) := Li{x) and Fi{<l>\x)) := Fi{x) {x e 



are then the angular momentum and the asymptotic direction at the pericentre 
X of the orbit. 

Clearly, by fiat, Hi,Li and Fi are constant on one orbit ^*{p,q), whereas 
Ti{^^{p, q)) is an affine function of time t. Therefore we have linearized the full 
motion. We must show that the functions (Ti, Lf, Hi, Fi), restricted to U[ , are 
indeed smooth coordinates. But this follows from the smoothness of the lifted 
functions in Lemma IXTl smoothness of time change, and smoothness of the KS 
flow (14. 6 p (which of course exists for all energies E E M). 

Since y = (Ti, Lf, Hi, Fi) defines a homeomorphism of U[ onto its image, 
we use it to define a differential structure on the whole of U[ and thus on P. 

Smoothness of the energy shells for > Vmax follows by noticing that 
these E are regular values of H : P ^ M.. □ 



6 M0ller Transformations 

Next we define the M0ller and scattering transformations which compare the 
asymptotics of the motion with a 'free motion'. We base ourselves on the articles 
|Hu] of Hunziker and [Sim] of Simon. The recent monograph [DGj by Dereziriski 
and Gerard treats these questions in the context of classical and quantum me- 
chanical n-body scattering. 

Due to the long-range character of the Coulomb interaction we cannot in 
general use the flow generated by the Hamiltonian function as 'free motion'. 
Instead, we compare with the Kepler motion generated by 

Hoo{p,q):=lp' + VUq} with V^q) = (6.1) 
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Thus we consider the smooth complete flow 

C --Poo^Poo (6.2) 

generated by (16. ID . 

• If Zoo = 0, $^ is the free flow on the phase space Poo := T*R^. 

• If Zoo > 0, we regularize T*(R^ \ {0}) in the way described in Thm. [5?I] 
to obtain Poo. 

• If Zoo < 0, then <l>^ is already complete on Poo := T*{R^ \ {0}), since 
particles of finite energy cannot meet the origin at g = 0. 

Thus we are to compare motions $* and $^ on the different phase spaces P 
and Poo. We cannot just identify P with Poo by neglecting the measure zero sets 
projecting to the singularities, since later on we will be interested in certain sets 
of measure zero like the bound states (moreover, these bound states will turn 
out to be crucial in our analysis of scattering, too). 

We overcome the above difficulty by observing that it suffices to identify Poo 
with P in a neighbourhood of spatial infinity. More precisely, let 

Poo,+ ■■= {xePool Hooix) > 0} (6.3) 

be the set of phase space points with positive 'free' energy. Then the orbit $^(a;) 
starting at x G Poo,+ goes to spatial infinity for large positive and negative times. 
The ball {^G M| | |g| < Pmin} contains all singularities of V and the singularity 
at the origin of the Kepler Hamiltonian Hoo- Therefore, we can canonically 
identify points (p, G Poo with points {p,q)EP if |g| > Pmin, and we denote 
this identification by Id. Thus the M0ller transformations 

■= lim $"*oIdo$*^ (6.4) 

are formally maps Cl"^ : Poo,+ P- They exist as pointwise limits, see Thm. 1^751 
below. 

First some standard definitions (see [Hu] ): 
Definition 6.1 

' , n I ^/ , TTT,+ \ • 1 111 7+ 1+ ^ ^ 

E 



:= {x 


G P 


b := b+ 


n 6" 


:= {x 


G P 


s := 


n 



, x) is bounded } , 6^ := 6 fl S 

(the bound states) , bE '■= b HT^e 

i b^ and E[x) > 0} , s| := n 

(the scattering states) , se '■= s ^T^e- 
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We shall show that = f2=''(Poo,+) so that the term 'scattering states' is really 
justified. 

Remarks 6.2 1) By continuity of $*, 6^ can be represented as the union b"^ = 
U^]^6^''^ of compact sets b^'''. Hence 6^ and are measurable w.r.t. Liouville 
measure ^ 

A := —u Auj Au 

on the symplectic manifold {P,uj)- 

2) The sets b"^ and are also $*-invariant. 

For all > the sets 6^ are closed, and bE is compact, being a subset of 
the compact region in T^e projecting to the ball {g G | |g| < -Rvir(-E')} in 
configuration space. 



Theorem 6.3 The limits 



= lim o Id o 



t^±oo 



exist pointwisely on Poo,+ C Poo 3nd ttius define the M0ller transformations 
: Poo,+ — ^ s^. These are measure-preserving homeomorphisms and intertwine 
and <l>^; 



" oo 



The asymptotic limits : — > and L^*" : — > of the momentum and 
the angular momentum L{p, q) := q x p 

P^{xq):= lim po<|)*(xo) and L^{xq) := lim L o $*(xo) (6.5) 



t— >itoo t— >itoo 



are continuous functions. If {po,qQ) = xq with go := \qo\ > Rvir{E) and 
± {%,Po) > 0, then 

p^ixo) =po + (l/(goV^)) , L^ixo) = L{x,) + O (l/(g^v^)) , (6.6) 

and for (Pq, Qo) := fi'^(R), gb) and E > Pth 

Po = + O (l/(gi+^V^)) , go = go + 0(l/(g^P)). (6.7) 
If e = 1 in Def. \2.1\ of Coulombic potentials, then for all energies E > 

VERmin C2 



< 



qo Ego 



Qo - qo 



(remark that by Def p.l4\} of C2 one always has < I). 



Eqo 
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Proof. First we show (16.51) . Since the system is reversible, we consider only the 
case t — > +00. For initial conditions Xq G there exists a time to with 

\q{to)\ > Rvir{E) and (g(to),p(to)) > 

for (p(t),q(t)) := $*(a;o), since otherwise \q{t)\ would be uniformly bounded as 

t ^ 00. 

W.l.o.g. we assume to = 0. Then by (12.160 



Thus by (EH]) 



|VF(g(t))| < 



■^00 I + CrR 



mill 



it > 0), 



the limit p^{xo) = \imt^^p(t) exists, and for this choice of xq 

p^ixo)~Po = Oil/{qoVE)). 

Being a locally uniform limit of the continuous functions xq i-^ p{t,xo), p^ 
is continuous. Next we show that the asymptotic limit Z+(xo) of the angular 
momentum exists. For t > we can estimate 



|VV(g(t))xg(t)| 



VViqit)) - Zo 



m 



xq{t) 



< C,R^^ \m\~ < Cl^min (go' + lEf)' 



using (12. 5 p and (12.160 . which shows the existence of the limit Z+(xo), and (16.60 . 
Continuity of L"*" follows as above. 

We now seek a Kepler hyperbola {P{t),Q{t)) := $^(Xo) which is positive 
asymptotic to (p(t),q{t)) and write 

r{t) ■.= q{t)-Q{t). 
Then r is a solution of the differential equation 



f{t) 



\m-mf 



W{q{t)) with lim f{t) = 0. (6.8) 

t—*oo 
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Setting for > 



'G C([0,oo),R3) 



|f|| := sup |f(t)| < min(i?^i„ C2/E) \ (6.9) 



(with C2 = 31(1 + l/e)Rl^^^Ci from (1Z141) ). by (IZTBI) and (1Z141) 



(^r)(t) := 



(is 



dr Z, 



q{r) - f(r) 



- Vl^(g(r)) 



{T)-f{T)f 

is well-defined for r G Ce, noting that by (12.141) 

|g(r) -r(r)| > R^-,^. 

We estimate |(JFr)(t)| as follows. The integrand of (I6.10p is bounded by 

g(r) - f(r) 



(6.10) 



|g(r) -r(r) 



- Vl^(g(r)) 



< 



< 



VV(g(r)) - Z. 



+ 



r r 



ClRn 



\q{T) - r{r)\-' |g(r)r 

9Cii?mm 



r 



\2+e 



min(|g(r)|,|g(r)-r(r)|)2+^ 



< 



r 



|2+e 



by the decay assumptions (12. Sp and (12. 6p (which are valid for the asymptotic 
potential, too, see Remark [Z3l 2). and Def. (16. 9p . For the last inequality we used 
|r(r)| < \ |g(r)|, following from (lZl4p . 
Estimating |g(r)| with (12.160 gives 



dr Z, 



q{T) - r{T) 



VV{q{r)) 



(6.11) 



J s 

^ 9(2 + V2)CiR^i„ f r- 
< = max I V -C/S, go 



dr 



E 



Hence a second integration yields 



llJ^f^l < 9(2 + V2)CiR^..^-^ < 



Eq'o 



(6.12) 
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and thus for e = 1 or > C2/Rmm or |go| > RmriE) large JF maps Ce into 
itself. Similarly, for fi E Ce 

II JFfi - J'f2|| 

Inir) - r2{r)\ 



oo /"OO 



Jo ^^min(|g(r) -ri(r)|,|g(r) -f2(r)|)2H 

POO POO 

< 8Ci||fi-f2|| ds dT\q(T)r^~' 







< 8(2 + v/2)C.i±^||f, - < ?^||f, - r-Al (6.13) 



By (1Z141) for e = 1 or E > C2/ Rnnn or |go| > ^vir(^) large the r.h.s. of 



is bounded above by |||ri — r2||, so that the map JF : ^ C^; is a contraction, 
and thus has a unique fixed point r G Ce- The first estimate in ( 16.71) follows 
from (16. lip and the second from (16.121) . 



After having shown unique existence of : — > Poo,+, we show unique 

existence of the M0ller transforms by the same method of integral equations, 
interchanging the roles of the two potentials. This is indeed possible with the 
same constants (see Remark IZ3l 2). 

The proof of the remaining statements is the same as in |Sim] . □ 

By energy conservation and (12.40 the modulus asymptotic momenta p^{x) equals 
\j2H[x), which is non-zero by Def. l6.1l of s^. So the asymptotic directions 

P^:s^^S' , p^{x):=^^ (6.14) 

are well-defined. 

Corollary 6.4 1. s"^ = f2^(Poo,+)- ttiat is, every positive energy orbit wtiich 
is unbounded in positive (negative) time is positively (negatively) asymp- 
totic to a Kepler hyperbola. 

2. The motion $* generated by the Hamiltonian function H : P ^ W is 
asymptotically complete, that is, up to a subset of Liouville measure 
zero the phase space consists of bound states and scattering states: 

X{P\{bUs)) = 0. 

Similarly for the Liouville measure Xe on the energy shell 

XE{^E\{bEUSE)) = {E>0). 
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3. The scattering transformation 

S ■.= n+ on- : D ^ (6.15) 

with domain D := C Poo,+ snd range ^^^(s) is continuous and 

^l^-invariant, i.e. 

5 o = o S. 



Proof. 

1. The equality of with i7^(Poo,+) is a consequence of Theorem 16.31 since 
nf = (n^)-^ is defined on 

2. This is a consequence the fact that C b'^''^ for t > but 
A(<I>*(6+''')) = A(6+'^), since $* is canonical (see [Hu]). 

3. Follows from Theorem 16.31 □ 

Theorem 6.5 Let V be a Coulombic potential whose partial derivatives decay 
at infinity according to 

{v{^ + ^) '='^ O (igT'"'-'-^) (/? e nl) (6.16) 

for some < e < 1. T"/)en t/?e M0ller transformations : Poo,+ — ^ are 
diffeomorphisms and canonical transformations. 

If iPo,Qo) = xo with go := |go| > RviviE) and ±{qo,po) > 0, then for 
multi-indices G Nq combined in 7 := {a, (3), 82^ := c^po^go 

(^^(^0) - L{xo)) = O (^g-l^l-^E-i(l"l+^))) , (6.18) 
and for {Pq, Qq) := n^{po, %) and E > ^th 

dliPo-Po) = O (g-l^l-^-i?-i(l"l+^)) , dliQ, - g-o) = O (g^ l^l-^-l"!/^-^ 

(6.19) 



42 



Remarks 6.6 1) The above decay condition (16.160 is met by purely Coulombic 
potentials ( 12.71) with e = 1. 

2) It is more natural to rescale the momentum by setting v := p/^/2E (and 
going to the rescaled angular momentum L/\^2E). In terms of these variables 
all estimates in Thm. [6751 show an energy dependence of the order 0{l/E). 

Proof. Note that all estimates of the theorem for 7 = coincide with the ones 
of Thm.O 

1) We estimate the derivatives of 

q{t, xo) =qo + tpo- / W{q{T, Xo))dT ds. 
Jo Jo 

We have for g := > 1 

dlq{t, Xo) = dliqo + tpo) - £ Vd2,V{q{T, Xo))dr ds 
= dliqo + tpo)- (6.20) 
E E r f D'^VViqir^Xo)) [dt^ir, Xo), . . . , dt;' q{T, Xo)) dr ds 

|7'^)|>0 

or 

(1 + Q)idim = dliqo + tpo)- (6.21) 
E E r f D'^VViqir^Xo)) {dfyir^Xo), . . . ^df^ q{r,Xo)) dr ds 

|7(')|>0 

with the linear operator Q = Q^^ given by 

Q{w){t):= [ [ DVV{q{T,Xo))w{T)dTds (t > 0). (6.22) 
Jo Jo 



We note that on the r.h.s. of (16.210 only partial derivatives or order |7(*)| < g 



appear, so that we can perform an induction in g, if we are able to invert the 
operator 1+ Q. 

We assume w E C with 



C := |w G C ([0,oo),M^) 



\w\\ X := sup \w{t)\/{t)x< 00 )■ (6.23) 
t>o 
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for 

(t)A := a/1 + (xty 

(note that C is independent of the choice of A > 0). 

Q maps C into itself, and we want to prove that for all positive energies the 
operator norm of Q is strictly smaller than one. 

If we assume 

qo>RviriE) and (go,Po) > 0, (6.24) 
then estimate (I2.16p holds: 



\qit)\ > qo ■ {t) X for all t > 0, with A := ^/Ws/go- (6.25) 
( 16.1611 implies that 

d^,V{q) = 0(|gr'^'-')- (6.26) 
So we find cat > such that 

p^VV(g(r))|| < c^(go(r)A)-^-' {N G N). (6.27) 

For estimating the norm of Q, we may restrict ourselves to w E C with 

||^||a = 1. (6.28) 

Inserting (K27\i and (KM into (K22\\ . we get 

\Q{wm\ < c^qf I I {r)fdTds<-^^Jj^ {l + u')-'duds 



^0 



CiVr CiTT , , 



^ql Eqo 



Assuming (16.241) . a suitable choice of the E'-dependence of i?vir consistent with 
assumption (12.140 then gives 

||Q.ollA<i and \\Q.,Jx = 0{l/{qoH{xo))) (6.29) 

for the operator norm HQ^^^oIIa of Q w.r.t. the norm || ■ ||a in (I6.23p . 

We return to (I6.2ip and estimate its r.h.s. Using (16.250 and (16. 27^ . we get 

D''VV{q{T,Xo)) (aj;^'g(r,Xo),...,aj;"'g(r,Xo)) dr ds (6.30) 

pt PS ^ 1 

<CNqo''-' / / {r)fdrds-]l\\di''qi-,Xo)\\x< const. {t)xqo^^^E-2\''\-\ 
Jo Jo 
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assuming 



\\dtqi;Xo)\\x = O for < \i\ < g. (6.31) 



Inserting (16.290 into (I6.2ip . we see that (I6.3ip holds for the start of the induction 



i.e. |7'| = 1. Also, (lOUj) is consistent with (Oil) , so that (OB holds for all 
multi-indices 7' G Nq. 

2) Now similar to (16.211) . the derivatives of the momentum meet the recursion 



dl{p{t,x,)~p,) = - (6.32) 
f D'^VViqir, xo)) (s^^^^glr, Xo), • • • , dt""' q{r, x,)) dr. 



|7(')|>0 



We insert (16.270 and (16.310 into (16.320 and performing the time t ^ 00 limit. 
Then by Lebesgue's Dominated Convergence Theorem ^^^^^^(xo) exists, and 
estimate (16.170 is valid, dj^^p^ As a locally uniform limit of continuous functions, 
d2gP^ is continuous. 

3) Similar to the proof of (16.31) we use the ansatz 

^xo-^L{x{t, Xo))dT 

t 

92o (Vl^sr(g(r, xo)) X g(r, xq)) dr 
J2 E r (VV;r(g(r, xo)) X g(r, Xo)) 



d2o^Qi^, Xo),---, d2'o'^qir, xq)) dr, 



in order to estimate the smoothness of the difference between the actual and the 
asymptotic angular momentum. 

Note that the short range potential Vsr = V — Voo appears in assumption 
(I6.16p . Thus for all t > the above expression is bounded above by 



l/3hi-.^-2l-l / {r)-,'-^dT < const. 
Jo 
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So by dominated convergence d2^L~^{xo) exists, and is estimated by (I6.18p . 
4) As in Thm. [Q] we now consider the Kepler hyperbola 

(p(t;xo),g(t;xo)) , = {Po,Qo) = n+{xo) 

which is positive asymptotic to $*(a;o) = {p{t, xq), q{t, xq)) and write 

r{t) = r{t; Xq) := q{t, Xq) - Q{t\ Xq). 
Then r is the solution of the integral equation ( I6.10p . Thus formally 

dim = P P dlv{v^{Q{T-xo))-V{q{T,xo)))dTds 

= P,{t;xo)+I]{t;xo) (6.33) 
with the integrals 

d2^VV,Mr,Xo))dTds 

and 

mt-xo) := /"/"^xoV {vMir;xo))-VUq{T,Xo))) drds. 

We estimate the first integral in a similar fashion as the r.h.s. of (I6.20p . but use 
the stronger estimate (16.160 for Vsr instead of (16.261) for V. We thus obtain 
(with g = |7|) 

\int;xo)\<i2 E 

N=l ^(i) + ...+^(JV)^^ 
|7(*'|>0 

r>00 POO , ^ 

drds 



D'^VVUqir, xo)) [0^' q{r, a;o), • • • , Xo)) 

noo ^ 
{T)f~'dTds-Y[ 
i=l 



N 

\ Ho 



< const, gg 2' I 
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Unlike J^, I2 depends on r. So we split it into 



with the linear operator 



00 /"OO 




DVVoo{Q{t] xo))w{T)dT ds {t > 0) 



N=l ^(i)+...+^(JV)^^ 
|7(i)|>0 




and 



d2n^Qi^, a;o), . . . , d^f^qir, xq) ] dr ds 



-'4 



E E 

N=2 ^(i)+...+^(]V), 
|7('-)|>0 

f 00 /"OO 




v(l) 



(iV) 



drds. 



On the space of bounded w, V = Vx^ is bounded, since for sup^ \w{t)\ = 1 



nCXD 
{T)lHrds 



1 - 



u 



du 



Assuming (16.241) . a suitable choice of the .E-dependence of _Rvir consistent with 
assumption (I2.14p thus gives 

\\Vx,\\<\ and \\Vx,\\=0{l/{q,H{x,))), 

so that we can invert 1 — P for all energies > 0. 

Il and Ij only contain derivatives d2lr{T;xo) with |7'| < (7. 
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Majorizing 14o(Q) — Kx.(g) by 2\DVooiq)'r\, the terms in I] are estimated by 

{'r)x^\riT;xo)\dTdsT\\\d2''fq{-,Xo)\\x 

since \ f{t,xo)\ < Rm\n (see (16. 9p ). and using (I6.3ip . 

Finally, inserting Q = q — r into the first A^-linear form of Ij and expanding, 
we see that each term contains at least one factor of the form df^^f{-, xq). So 

n = O [q-\'\-'E-\\-\~'^ , 

too. 

Together with the above estimates for and this finally shows that 

dlA-^: xo) = (11 - v)-\r, + u + m = o{q,\'\~^E-\\-\~\ 

which is equivalent to the second assertion in (I6.19p . By a parametrized version of 
the Banach Fixed Point Theorem (see, e.g. [DGj, Prop. A. 2. 2), xq ^ d2^^f{-;xo) 
is continuous. 

The first assertion in (I6.19p follows similarly, since 

dl{Po-Po) = - d2,V (vUQiT;xo)) - Viqir^xo))) dr. □ 

7 The Flow between Near-Collisions 

We remind the reader of the relation (I2.13p for the virial radius Rviy-{E). We 
assume that Eth > Vmax so that we may assume -Rvir(-E) = -Rvir for all E > Eth- 

By the virial identity (I2.12p a trajectory of energy E leaving the interaction 
zone IZ = {q E M^i \ \q\ < i?vir(-Eth)} cannot reenter it. 

We show first that a particle cannot stay inside the interaction zone for a 
long time without having close encounters with the nuclei. Then we will control 
long trajectories within IZ by Poincare section techniques. 

To quantify this, let 

n 

XZ{r) ■=IZ\[jmi{Bi{r)), 

1=1 



48 



Bi{r) = {g G ]R| I \q — si\ < r} being the ball of radius r around the Zth nucleus. 

Instead of considering the restriction of the flow generated by H, it is 

technically convenient to consider the flow $^ generated by 

Hj,{v,x):=l\v\' + ^ on {He)-\1). 

We then have for {p{t),q{t)) = $*(po,go) and {v{s),x{s)) = ^%{vo,xo) 
with 

{vo,Xo) := {po/V2E,qo) 

= (^V2Ev{V2Et),xiV2Et)) . (7.1) 
In the lemma below we use the standard Euclidean metric on ]R| x M|. 

Lemma 7.1 For E > Eth the rescaled flow in IZ{cg) is C^-near to the 
linear flow generated by the Hamiltonian function = in the sense that 

sup {\<!>%{vo,xo) - {vo,xo + svo)\ + \\D<^%{vo,xo) - 
se[o,T] 

= 0{1/E) (7.2) 

x{[0,T])cIZ{c,) (7.3) 
Moreover, condition ( 17.31) is never satisfied ifT > 3-Rvir- 

Proof. We assume that T < 3-Rvir. prove estimate (17. 2p and then the necessity 
of that assumption. The constants 

Lo:= sup \V{q)\ , L, := sup \VV{q)\ 

and 

L2:= sup \D^V{q)\ 
are finite, since XZ{cq) is compact. We bound the differences 

Zi{s) := {v{s) - vq, x{s) - (fo + svq)) 

and 

Z2(.):=D<l>^(t;o,Xo)-(AS) 
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in (17.21) using their integral equations 

ns pu 

x{s) - (fo + svq) = -{2E)~^ j I VV{x{t)) drdu 

Jo Jo 

and 

_i r fo d^v{x{t)) ' 



and thus obtain 







T + T'^/2 T 
sup < Li , sup |Z2(s)| < L2— . 

se[o,T] se[o,T] 

This is indeed of order 0{l/E). 

Moreover, for E > Eth > max(4Lo, Li_Rvir) we have 



l-V{xo)/E> 2 
and \x{s) — {xq + svo)\ < -Rvir/4 for s < T < 3-Rvir so that for T := 3/2vir 

Thus for T = 3-Rvir we obtain a contradiction with our assumption (17.30 . since 
the diameter of equals 2i?vir so that xq and x(T) cannot be both G IZ. □ 



8 The Single Scattering Process 

We now consider motion inside the ball Bi(cg) near the /th singularity, using the 
KS estimates of Prop. I4~ll For initial conditions in the phase space region 

Vi := tt{Vi) = {xeP \ \qix) - si\ < Cg, H{x) > Ett} 

over the ball the exit times 

Ti^-.Vi^R , Ti^ixo) := ± inf{t > | \q{±t, Xq) - Si\ = cj 

are well-defined and smooth. We compare them with the exit times T^^ for the 
purely Keplerian motion in Bi{cq) generated by (15.41) . and set 

^f:Vi^dVu x^^{T^{x),x) 



50 



and vl>±^(x) := $i(T±,(x), x). 

Similar to (14.161) we define the pericentric tiypersurface 

Hi := TiiUi) (8.1) 
= {x eVi\ q{x) = si or {p{x), q{x) - Si) = 0} 

on which the pericentric time T° : Vi ^ M. vanishes. If > we smoothly 
extend T° by demanding that 

^^{x) ■=^T^{x),x)enu (xeVi). 

In analogy to (14.131) we introduce the diffeomorphism 

v^. : dV, X ^ (p, q) ^ I "^Xfl ' f: {- 5| ^ I (8.2) 

^^''^^ \ (x) > ^ > 

that permutes incoming and outgoing data and its analog ^ l i for the Kepler 
flow. 

We first consider the Kepler flow (with potential q ^ —Zi/\q — si\) with 
incoming, resp. outgoing coordinates 

^ ^ e&Dir\^E and (^1]:= ^ li ^ ^ 



q J \<i J \(1 

It is more convenient to work with the coordinates 



and w^:=^ -. (8.3) 



V2E 

Lemma 8.1 Ttie Kepler transformation is given by 

V+{v ,W } = ^ i_^2 — 

^ + (^- = 2m(1-m2(1+/3/2))i7- + (1-«2(1+/3))«;- 



(8.4) 



wliere 



/5 := and u = u{v^,w ) := — ^^ ' "-.J = —u{v^,w^). (8.5) 
CgE 1+2^ 
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Proof. For simplicity of notation we assume si = 0. Since the KS transformation 
(13.21) transforms the scalar product in the formula (14.151) 

u 



\P\' + \Q\ 



2 



into 

tr(Pg*) tr(M*) 



4V2E 2V2^ V2S 

l-Pp = Cg|tT~p = Cg(l + 13) and IQP = Cq, we see that the definition (18. 5p of u 
is consistent with the definition in (14.150 . 

Eqs. (18. 4 p are obtained by inserting (14.150 into the KS formula (13.20 for 
noticing that 

^ 2E 2E 
follows from the three-term-identity 

YXY = 2{X*,Y)Y - {Y,Y)X* (X, F e H) 

of quaternions, applied to p,q E ImH. □ 

As we assumed that the NC-condition holds (no three singularities on a line), 
bounded orbits must be scattered by an angle <(v'^,v^) < i!) which must be at 
least of the order of amin > (see (12.100 ). This can only happen if the mismatch 
between the initial velocity tT" and position is only of the order 0{1/E). We 
first show this for Kepler scattering: 

Lemma 8.2 For all energies E > Ef^ > 4Zmax/cg, angles 6 e [n ■ Etii/E^n], 
and initial conditions (p^, q^) E dVi fl S^; and m > with 



> (8.) 



Kepler scattering is in the forward direction, namely the total scattering angle is 
bounded by <(v'^,v^) < 9. 

Proof. Since the parameter m > from (18.50 equals 

{v-/\v-lw-),/TT^ _ ^ 

u{v ,w ) = <-{t^,w 
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so that 

1 - ^2 = (1 + u){l -U)>l-U>l 

we get from formula (18. 4p for v'^ that 



w 



(8.7) 



1 -H- 1 


up 1 ^-1 


Itr — V 


= :; :;\uv +w 




1 — 



< 



+ w' 



+ \v-\{{l-u) + \l-l/\v-\\) 



with p from (lO) . 

For E > £;th we have \(3\ < \Zi\/ {cqEtYy) < \, so that \v~\ = y/l + fJ G 
[|,2]. Thus inserting the last two estimates of (18. 7p and (18. 6p gives 



1^ < 2|/3| 



< 



Cn 



TT ■ ^th 

e 





1 - 



1 



_^ 2|Z,| 02 



E, 



Cg TT^-Eth 



e e 
< — + — 



2tt 2ti {2 



^ < —e 

47r 



This proves the assertion, since 



IT 


v+ 




V 


2 







< -{1 + <9. 



□ 



Later on we will study the linearization of the flow $*, using Poincare section 
techniques. Therefore we now calculate the linearization of the Kepler transfor- 
mation (18.41) . for tangent vectors {6p^,5q^) in the four-dimensional subspace 

%-,,^-)(9AnS£;), 

that is, for variations {5v~,5w~) of (18.30 meeting 

{^v-,5v-) = , {w^,5w')=0. (8.8) 
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Lemma 8.3 The linearization of the Kepler Transformation ( |g.4)) is given by 

u'^Ml - u'^) - u^f3)6v- - u(36w-] - (3[2uv- + (1 + u'^)w-]du 



and 

5w' 



with 



(1-^2)2 

(2n(l - u^) - u^(3)5v- + ((1 - u^) - u^f3)5w- 
1 -m2 

[(2(1 - + Pu\u^ - 3))v- - 2uf3w-] du 
^ (1 -n2)2 



du :- 



Proof. By differentiation of ([H3D. □ 
As the next lemma shows, there are two regimes of Kepler scattering: 



• hard scattering: If the scattering angle is larger than ci/y/E, then lin- 
earized scattering is basically a reflection combined with scaling. In partic- 
ular, for an energy-independent scattering angle the Liapunov exponent is 
of the approximate size E. 

• soft scattering: If, however, the scattering angle Aip is smaller than 
C2/VE, then linearized scattering is a perturbation of free motion. 

Lemma 8.4 1) For ci > and < 6 < ^ we consider initial conditions 

{p~, q^) E dVi HTje leading to a scattering angle Aip := <{v'^, v^) > ciE^^ . 
Then the linearized Kepler Transformation of Lemma \8.3\ is estimated by 



6v+ \ _ 4cgSin2(iA^) / R R\ / Sv 
5w+ ~ -Z, [ R R [ 5w 



6v 
6w 



■0{E^) (8.9) 



where R = R^+_^~ E 0(3, M) is the reflection by the plane perpendicular to the 
vector — v~ . 

2) If instead Aip < C2E-^ with | < 5 < 1, then 

Sv^ = (5{r+ 0{clE^~^^)-{\6ir\ + \6w-\), 

5w+ = 2u6ir + 6w- +2v- du + 0{clE^-'^^) ■ {\6ir\ + \6w -\). {8.10) 
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Remark 8.5 For ^ 7^ | the constants ci,C2 become irrelevant (for a threshold 
energy Eth > 1). 

Observe that the leading term in (18. 9p scales like c\E^^'^^ ^ 1, whereas the 
error term of (IHTigi) scales like clE^^'^^ < 1. 

Proof. 1) First we prove the identity 



,2 ^ 



with the eccentricity 

e = ^l + 2E\Li\yZf = l/sm{lAij) (8.12) 

of the Kepler hyperbola. By reinserting the (p, g)-coordi nates with the help of 
(18.31) . def. (18. 5p of u acquires the form u = —2 f \ (assuming without 
loss of generality that si = 0). We now insert the relationship 

(r,r>' = in' in' - ^n 2(e+ zz/c,)cj - 

into the expression for 1 — to deduce (18. lip . 

If > cE-\ 

1-m2 = 0(c2e2(^-i)) (8.13) 

so that 

^ JPf + _ /^(^g- + f -) ^ oiE") . (I^n + 1^.-1) (8.14) 

(1 — 1 — u"^ 

and 

2l3(v' + w~)du l3(6v' + 5w~) „/^n\ /le-^i le-^-ix 
=-^^7^ ^ 2 - + 0{E^)-{6v +6w . 8.15 

Using (ES) we get 



1 + 1/3 

The vectors w~ and v~ in the last formula are nearly anti-parallel. From the 



expression (EI]) for v+ and ( I8l3i) we see that = v' - ^^V'^a ^ + ^'(1/^) 



or 



^- = (l-n^)(^-y^ + 0(l) 
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so that 

5v^ 



2{v- -V+ + 0{1/E)) {v- -V+ + 0{1/E), 6v- + 5w-) 



f3{l + P/2) 



2{v - (v - 6v + 5w ) P{5v~ + 5w ) 
■ {\6v~\ + \6w"\) 



Extracting the square of the norm 



\v — irl 



\v^\ \v 



+ 0{1/E) 

= ^2(1 - cos(AV^)) + 0(1/ E) = 2 sin(iA^) + 0{1/E), 
we obtain the first estimate in (18. QD . since by (18. lip and (I8.12p 



P 



1 - m2 /3e2 
and the reflection equals 



4sm2(iA7/;) 



+ 0{E'), 



R^+-^- —11 — 2 



\v — v'^) {v — v'^l 



\v~ — v 



Estimate (lO) for 6w^ follows since the r.h.s. of (lg?15l) and (lg?14l) have the 
same form. 

2) We write 

5v^ = 5v~ + Ri + R2 , Sw^ = 2uSv~ + Sw~ + 2v~du + Ri + R^ 
with rest terms 



Ri :-- 



R2 ■.= P 



(3 



1 — V? 



6v~ + 



6v~ + 6w~ + 2 



[v + w )du 



6w~ v~du {v^+w~)du 



1 + u (1 + m)2 
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and 

-\- u + _ _ ^_ w~du {y~ + w~)du 



R3 ■.= P 



-6v - 6td + V du+ . , „ 
1 + u {1+uy 1-u 



2 



Now we assume that < C2E ^ with | < 5 < 1. 

So by the first part of the proof i?i = 0{4E^-^^) ■ {\5v-\ + By 
the same reasoning the (identical, up to sign) last term in R2 and R3 is of order 
0{1/E) ■ {\5v^\ + \Sw^\). The other terms in R2 and -R3 are of the same order, 
since they contain the multiplier p. □ 

We now use the results of Section |4]to compare the true motion with the Kepler 
motion of the last lemmata. 

Proposition 8.6 For pericentric initial conditions xq E Hi witti energy E : = 

H{xo) > Eth 

T^{xo)=Ti^,{xo) + 0{E-'/'), (8.16) 
and for {p^,q^) := (xq) and := ^f/xo) 

{p^/V2E, g±) = (pt/V2E, gt) + 0{1/E). (8.17) 
For\qQ-si\ = Oil/E) 

TtM = ± - J^2HEcJ\Zi\)^ + 0{E-'/'). (8.18) 

For arbitrary xq G Hi ttie r.ti.s. of ( Ig.jgj) is an upper bound for |T^^(xo)|. The 
diffeomorptiism 

onto /ts image wtiicti conjugates tfie maps = (S^)^^ o oEi) is C^-near 
to tlie identity in ttie sense ttiat in ttie Euclidean norm for the {v, w )-coordinates 

\Ei{x) -x\ = 0{l/H{x)), (x G m). (8.19) 
and the solution lj of the linear problem is C^-near to '^i in the sense 

\^i{x)-^L^ioEi{x)\ = 0{H-\x)) (8.20) 
\\Dmi{x)-DmL,ioEi{x)\\ = 0{H\x)). (8.21) 
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Remark 8.7 In particular the exit times (18.181) are independent of xq, up to 

Proof. Estimate (E171) follows directly from (140^!) . 

Estimate (14.391) says that for KS initial conditions Xq with 7r(Xo) = Xq the 
relative error 



\Qis,Xo)\ 
This translates into the estimates 



\qit,Xo) -q,{t,Xo)\ ^ , 

\qit,xo)\ ^ ' ^ ^ 

and 

\p{t,Xo)-pLit,Xo)\ _ 1 

^^^^-^ ) (t^O). (8.23) 

We know that inside Vi each trajectory t $*(xo) passes the hypersurface Hi 
only at xq, so that we can parametrize the curves (0,T;~^(xo)] 3 t ^-^ ^^(xq) = 
(p{t),q(t)) by their radius r(t) := \q{t) — si\ (similarly for negative times). 
The explicit formula for the exit time 

(with the smooth potential Wi defined in (I4.ip . and relative angular momentum 
Li in (15. 6p ) is obtained by integrating the inverse of the radial velocity {p,q) /\q\. 
Similarly 

inserting (18. 22^ and (I8.23p then shows the assertion (18. 16^ . 

(I8.18P follows by evaluating the explicit formula (15. 9p for the exit time of the 
Kepler flow: As Er'^^^ + Zirj^i^ — ^Lf = 0, the square root terms in (15. 9p vanish 
for the minimal radius. Further, by our assumption on go we have rmm < C\Zi\/E 
for some C > 0, so that 

< W- ^11 V = 2\q- smq- si\ + Zi) < ^^'^^^ + = 0{l/E). 
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Est. (Km follows from KTM and similarly (K7m . (KTH} follow from (\A7M 
and (I4.2ip . since on dT>L i the KS transformation from coordinates {P,Q) to 
the (iT, w )-coordinates (18. 3p takes the simple form 



V 

w 



Q*hP 
Q*hQ 



□ 



9 Long Paths Within the Interaction Zone 

Strong changes of directions only occur if the pericentric distance of the orbit 
from a nucleus is of order 0{1/E). 

To quantify this, we set for < < vr 



Hiie) ■.= {xen 



(with the pericentric hypersurface Hi defined in ( 18. ip ). 

Lemma 9.1 For C5 > large the following statements are true. Consider 
pericentric initial conditions xq E Hi of energy E := H{xq) > Eth and 6 G 
[C,/E,n]. 

• If Xq ^ Hi{9), then the directions p"^ := ]3(^f(xo)) of the orbit through 
Xq at the moment of entering, resp. exit from the ball Bi{cq) differ at most 
by the angle 

Aip = <{p-,p+) <9. (9.2) 

• Conversely if xq E Hi{6), then 

AiP > e/A, (9.3) 

and positions and momenta are nearly anti-parallel before scattering and 
parallel after scattering: 

< (qi^fixo)) - si, ±p^) = 0{{eE)-'). (9.4) 
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Proof. Because of estimate (I8.20p of Prop. 18.61 we need only choose C5 large 
and then show the following statements for the Kepler flow (which correspond 
to (Q, (Q resp. (IMD): 



A^l:=<(pZ>Pl) < ¥^ (9-5) 

sini^AipL) > 1.05 ■ sin(^/8) (9.6) 

and 

<{^-si,±pt)=Oii9E)-'). (9.7) 

with pf := Pi^t/xo)) and ^ := g(^ti(^o)). 

• We use the explicit formulae of Lemma 18.11 to compute 



1 - 



2 c2 



using the substitution (18. IIP . Now the eccentricity e of the Kepler hyperbola 
equals 

e=\l + 2Er^jZil (9.8) 

as one sees by comparing the definition ( 15.80 of the pericentric distance rmin with 
(I8.12P (formula (15. 8p is also valid if Zi < 0; then the energy is always positive). 
Thus 

u{l + l(3) 1 



< \2Er^,jZi + l\-' <-—-^p^^-— < smile) 

2 + sin(i6')sign(Z;) 

since u(l + |/3)/a/1 + equals the cosine of the angle between v'^ and and 
is thus smaller than one, see ( 18. 5p . Then ( 19. 5 p follows, since the sine function is 
monotone increasing on [0,7r/2]. 
• By (Km and (TO 



|1 + sign(Zi)sin(^/4)| 



^ sin(^/4) ^ sin(7r/4) 2sin(^/8) 



1 + I/V2 2 sin(7r/8) 1 + l/^/2 
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showing (19.61) . 

• With ([H3]), (jgUi) and (El) 



sin (< (a;o)) - si, ±p^)) = sin (< {w^, ±v^)) 



1 - 



w , w 



-) 



+ Oil/E) 



1 -|- 'lEVjoin 



0{1/E) 
/Zi 



2/(3+1 



+ 0{1/E) 



0{l/E) 



+ 0{l/E) = 0{l/{eE)l 



2cjZi + l/E 

showing estimate (19. 7p . □ 

Long trajectories in the interaction zone XZ must have close encounters with 
singularities of distance 0{1/ E). To show this, we now assume that the radius 
Cq of the balls Bk{cq) is so small that there is no straight line meeting more than 
two balls. To be concrete, we assume 



Ca < 



1 



(9.9) 



Proposition 9.2 For Ett large we consider a trajectory segment 

[0,T] 3 t ^ $*(a:o) = G V (9.10) 

with 

V:={xeP \ H{x) > Eth, q{x) G IZ}, 

wtiicti does not intersect any pericentric hypersurface 7ii(amin/2), (/ = 1, . . . , n). 
Tlien, in configuration space, it does not reenter a ball Bk{cq) after leaving it, 
and it does not intersect three or more such balls. 

Furthermore, the length of the time interval is bounded byT < ISR^h/V^E. 

Proof. • There is no subsegment [^1,^2] 3 t q{t) of the trajectory, lying 
in XZ{cq), which leaves and then reenters a given ball Bk{cq). Namely we may 
otherwise assume 

q{[ti,t2]) <ZlZ{cq) , \q{ti) - Sk\ = \q{t2) - Sk\ = Cq and > ^i- 
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Then there would be a time ^ ('^^1)^2) with maximal distance |g(to) ~ Sk\ > Cg 
from Sk- 

d (f 

^|g(to)-4r = and — |g(to)-4r<0. (9.11) 

But 



1^ 



> 2{E - Lo) - sup {VV{q), q- 4) > - Lq - Lii?,;^) 

ge2Z(cg) 

with Lq = snp^^jzic,) l^io)] and Li = sup^-.gj^(,^) |Vl^(g)|, so that (,|9JJJ) does 
not hold for E large. 

• We now prove that there is no orbit segment (IQ.lOp whose configuration space 
projection intersects three balls Bk^i^Cq), Bk^{cq) and Bk^icg) in succession. By 
going to a subsegment, we otherwise assume outgoing initial data 

m^dBk,{cg) . (p(0),g(0)-4o) >0, 

ingoing final data 

qiT) e dBk,iCg) , {piT),q{T)-Sk,)<0, 
and intermediate times 

0<r<t+<T with q{t^) e dBk,{cg) and ±{p{t^),q{t^) - Sk,) > 0. 

By the first part of the proof we know that ^ ki ^ /c2. Then the angle 7 
between the direction s'^o-'^i of the axis through the centres and p{t^) is bounded 
by 

7 < |amin. (9.12) 

Namely all straight lines intersecting Bk^^{cq) and Bk-^{cq) have a direction whose 
angle 7' with s''0''^i is bounded by 

sin(7') < jj— < |sin(amin) < sin(iamin), 



using ( 19.91) . so that (19.120 follows for _Eth large from Lemma [7~H 



We now consider the unit vector u = Xqs''^''''^ + X2s'^^'^^ , K > 0, which is 
perpendicular to s^o-'^a 
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Then by (19.121) the angle between p{t ) and u is bounded above by 

<{p{r),u) < <(s'"''*^\m) +7 = vr - Itt - a(A;i, A;o,/i;2) +7 

< |7r - amin + 7 < - amin) (9.13) 

(with cos{a{i,j, k)) = (^s^'\ P''^)). 

On the other hand for some time t E [t'^,T] we have 

because otherwise there is no trajectory between q{t^) and q(T). 

But by Lemma [7TT] in the interval t E the change of direction 

<{p{t+),p{t)) = 0{l/E), 

so that with (19.130 and (19.141) we must have a large scattering angle 

<(p(r),p(t+)) > ia^in. (9.15) 

within Bk^{cg). 

Our assumption that the orbit segment (19.101) does not intersect the pericen- 
tric hypersurface 7Yfci(«min/2) (see (19. ip implies by Lemma l9H that the deflection 
angle is bounded above by 

contradicting (19.151) . 

Thus a trajectory segment can be represented as the union of at most two 
segments inside balls Bk.{cq) and three segments in XZ{cq). 

• By (18.161) and (18.181) for large -Eth the time interval spent by the trajectory 
inside the ball Bk^Cq) has length T+ - T^" < ^ < By Lemma O with 

(17.11) the time interval of the segments in XZ{cq) are of lengths < 3i?vir/ v^2£'. 
So the total time T is bounded above by 

T < 2^+3^ = 134^^. □ 
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We now control the paths between successive close encounters with the cen- 
tres Sfc. For radii Cg > meeting (19.90 the balls Bk{cq), k e {1, . . . ,n} around 
the centres Sk do not intersect. 

By the NC condition 12.41 the distance dist(si, A'^'^), i =^ k,l between the 
centre and the axis 

A'''' := {tsi + (1 - t)sk \te[0, 1]}, {k + I) 

connecting Sk and si is bounded below by sin(ainin) ■ d^m- 
Thus the configuration space cylinders (see Figured]) 

Zlf := {qe ml I dist(g, A^-'O < Cy, \q- 4| > Cg < \q- Si\} 
have empty intersection with non-adjacent balls (2^/ H Bi{cq) = for i 7^ k,l) 

if Cy < Cg. 



The cylinders themselves are not mutually disjoint; 
in particular = Zj^. 

However, under the NC condition 

if the axes A''^^^' ^ A^^'^^ are parallel. 
For the choice 

of the radius we have in addition 

z^^^nz^^^ = {k^h). 

Finally there is a minimal nonzero angle between nonparallel axes. More 
precisely, for this choice of Cy a cylinder Z\l'^^ can only intersect a different 
({fci,/i} ^ {fc2,/2}) cylinder Z^^/^ if the angle arccos(| (s'^i-'i, s'^^.ia^ |) 
between their axes is larger than amin- The proof of this fact is a nice 
exercise in conic sections. 



Thus for 
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Figure 1: Configuration space projections of Poincare sections 



and E > Eth the phase space regions 



Z 



kJ 



|x = (p, q) 



e V 



p 



7Vm 



x] 



(with V = {x e P \ H{x) > Eth, q{x) G IZ}) do not intersect (Z'^i''^ nZ'^^-'a 

0for {h,h)^{k2,l2)). 

For k ^ I \Ne erect Poincare hypersurfaces Ti^'^ near the midpoint 

rnf''^ ■= |(4 + Si) e A'''^ 

by setting for Cq > 

n'JiE) := {qe M| | (g - m'^'', S'^'^) = 0, |g - m'^''\ < |cp ■ rf'^-'eo/i?} , 
so that n''J{E) C Z^/ for large and E > Eth, and 



P -^ki 



< IcpCo/E 



(9.16) 



so that H^'^ C Z^'\ too; see Figure m The are four-dimensional submani 



k,l 
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folds without boundary of 'Ze- Sometimes we work with their closures 

n'j^', (9.17) 

which are diffeomorphic to products of two closed disks (and as such formally 
speaking no manifolds). 

High energy orbits between near-collisions are nearly straight and thus move 
near some axis: 

Lemma 9.3 Let Eth snd the constant eo > in ( IP. J 6]) be large and consider 
for E > Eth trajectory segments of ttie form 

[-T,T] 3 t ^ x(t) = (p(t),glt)) := $*(xo) 

with initial values xq E T.e, starting and ending in the pericentric hypersurfaces 
(EH); 

but not intersecting a pericentric hypersurfaces in between: 

a;((-T,T))n7^K"min/2) = (/ = (9.18) 
Then the trajectory segment is contained in the phase space region 

xi[-T,T]) cZ''-^'+UVk_UVk^, (9.19) 
and there is a unique to E (— T, T) with 

^'"{xo) eHe''^^. (9.20) 
Conversely for ^ k_ the exit time 

yfc_,fc+,± . ^ ^± ^ ^k.,k^,±^^^^ _ ±inf {t > I <l>±*(xo) E dZ''-''+} 

from Z^-'^+ is of order 

\2E 

and the exit points {p^,q^) := <l>(T^"''^'^'^(a;o), xq) are estimated by 

{p^/V2E,q^) = + (9.22) 

with the intersection point n^'^ := si + Cqs'''^ between the axis A^'^ and the sphere 
dBi{cq), see Figure\^ 
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Proof. 

By straight line geometry and Lemma 17.11 trajectories t ^-^ q(t,xo) with xq E 



''^^ intersect dZ^^ near n^"^'^^ , and 



p^/^2E - s 



{p^/V2E-po) + {po-s 



0{1/E), 



again using Lemma [771] and the definition (19.160 of H'^ '^^ . This shows (19.220 . 
Since \m''-'^+ - n^^^''^ \ = \d^'^^+ - c,, (ICTl) follows from Lemma O 

Concerning the first statement of the lemma, we know from (19.40 of Lemma 
EI that for t± := (x(tT)) 



:= $ (t^x(±r)) G {q\^,)-\dB,,{c,)) 

have the property 



'2E 



p^ + 0{l/E). 



(9.23) 



Furthermore, by the assumption ( 19.180 and the first statement of Prop. [9^ the 
trajectory lies in XZ{cq) during the time interval [— T + t+, T — t^]. Thus (17. 2p 
and (O) imply that 

p+ p^ 



and 



V2E 

g+ - q- 



'2E 
p- 



\q- 



'2E 



+ 0{l/E) 
+ 0{l/E). 



(9.24) 
(9.25) 



Using ( ICT) . and then (ICT) . 
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- SkJ) + (g 



V2E 
2c,+ |g 



/2E 



p^ + 



so that 



p^/V2E-s 



fc_ 



A second application of (19.250 shows that the second term in (19.220 . too is 

0{llE). 

By Lemma Owe obtain (19391) and (19:201) . □ 
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10 The Poincare Map 



For E > Eiii consider the Poincare surfaces 



n 



'He := U n'J 



n% := He U dV% 



(10.1) 




with n)i defined in and 

VE:=Vn IIe, crD% := {(p, g) G SVe | ± (p, g) > 0} . (10.2) 
Then the return time to the Poincare surface Tig 

T^,:Ps^[0,oo)U{oo}, 
is defined by Te{x) := for x G dV'^ and 

Te{x) ■= inf {t > G } , (x G I)e \ 9©+). (10.3) 

Lemma 10.1 For E > Eth ttie Poincare return time Te is finite and 



Proof. Let xq G \ dV^. By Prop. O for a time h e [0, 13R^ir/V2E] the 
point xi := (xo) 

• either exits Ve, that is xi G <9Pg 

• or meets a pericentric surface near the A;th nucleus: xi G 7Yfc(«min/2). 

In the first case we are done. In the second case we iterate the above argument 



and find t2 G [ti,ti + ISR^ij-/ V2E] for which the point X2 := $*^(xo) either 
exits Ve {x2 G SP^) or meets a pericentric surface near the / ^ kth nucleus: 



Te = 0{1/Ve). 



X2 G 7Y;(«min/2). 



In the relevant second case we know from ( I9.20p of Lemma [975] that 



X3 := $*nxo) G 7^^'' 
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for some ^3 G (^1,^2)- D 
We shall analyze the Poincare map 

Ve-He^ n+, Ve{x) := ^Te{x),x). (10.4) 

When we work with the closures 



n 

He := U He , n% ■= "^e U &Di 



k^l 

(see ( igiTi) ). we write Te, Ve etc. 

Lemma 10.2 Ve is a bijection, and its restriction to 7i^\P^^(97i^) is smooth. 
Proof. 

Ve is one-to-one: \f Ve{xi) = Ve{x2) and Te{xi) < Te{x2), then xi = 
$*2-*i(a;2) and 

• either xi E dV^- Then X2 = Xi since by the virial inequality 

nSI)^ = for t>0, 

• or a;i G He- Then by definition of Te using the infimum, X2 = xi, too. 
Ve is onto: T"/me reversal 

TJZ -.P^P, {p,q) eT*M ^ {-p,q), xEP\T*M^x (10.5) 
is a smooth anti-symplectic transformation, with 

T7^ o $* o T7^ = (t G M) 

and 

TJZ{n^) = HI. 

Thus 

Ve^{x) = TRoVeo TR{x). 

Smoothness of Ve\u U := Tig \ Ve^^OT-Ce) follows from transversality of 
the codimension one 9-manifold U C S^; to the flow ^^ts^, since Ve\u maps 
to inner points of 7i^. □ 
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Symbolic Dynamics 

In order to use symbolic dynamics, we introduce symbol sequences 

k = {ki),^j e 

over the alphabet 

S := {1, . . . ,n}, 

where 

I = II := {ieZ\l<i<r} (10.6) 

for Z,r G ZU {±00} is a (finite, half-infinite or bi-infinite) interval, 
k is called admissible if ki 7^ /cj+i for all {i,i + 1} C /. 
For E > and (/cq, ki) admissible we set 

VE{ko, k,) := WE{ko, k,) := H'^''"', (10.7) 

and for • • • , ^0) admissible 

WE{k.^, . . . , fco) := WEik.i, fco) n VE{WE{k.m, • • • , (10.8) 

resp. for (/cq, • • • , ^m) admissible 

VE{ko, ...,kj:= VE{ko, fci) n V^'iVEik,, ...,km)), (10.9) 

m > 2. Then by Lemma [10.21 the iterated maps 

VE^ko, . . . ,km) : VE^ko, . . . ,km) ^ WE{ko, . . . ,km) 
VE{ko, . . . , km) ■= ^\vE{ko,...,km) 

are diffeomorphisms (we will show in Prop. 111.51 that the sets VE{k) are non- 
empty). Again, WEik) and VEik) denote the closures of these sets. 
We decompose the Poincare map VE{k^i, k^, ki) in the form 

PEik_i,ko,ki){x) = VE{k-i,kQ)oV'E{k-uko,ki)oVt{ko,ki){x) (10.10) 
= p-oP^oP+(x) (x G Fs(fc_i,A;o,fci)) 

with the diffeomorphisms 

VE-VE{k^i,k^) \4(fc_i,A;o):=P^(VE(fc_i,fco))caz'^-^''=« 
X ^ $ (^T^-^''«'+(x),x) , 
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and 

P^:\4(fc_i,fco,fci)-W^^(A:-i,fco,A;i) , x^^i{x) 
with ^+ defined in ([82])) for 

\4(fc_i, fco, fci) := VE{VE{k-u ko, ki)) C /co), 

are considered as perturbations of Poincare maps for the free flow, whereas 
Prop. l4.11 allows us to view V'e as a perturbation of Poincare maps for the Kepler 
flow. 

Adapted Coordinates 

We now define adapted coordinates in the Poincare sections He- This will 
simplify to introduce invariant cone fields later on. 
We complement the unit vectors 

= s^{k_i, ko, ki) with := s^'±i''=« 

pointing towards by a unit vector 

= i^{k^i, ko, ki) with = i~ perpendicular to span(s", s^) (10.11) 

(by the NC condition 12.41 the span is one-dimensional iff ki = A;„i), and set 

= u^(A;„i, ko, ki) := x t^. (10.12) 

Then and u^, considered as elements of Tg-]R| for q E Ti^/^ '^", form an or- 

thonormal basis of Tg-Tijy/ ' . We introduce adapted coordinates on Ve(A;_i,/co) 
and WE{ko,ki) by mapping 

VE{k_^,ko) ^ R^xM^ , (^-^ ^ (^-^ ^-), 

WEiko,k,) ^ R^xR' , + 

with 

and := |m*=o-±i - 4ol = Irf''*^'*^^ 
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Remark 10.3 The coordinates {y~,z~) on Ve(/c_i,A;o) depend through the 
vector £^ = i~{k_i, kg, ki) on the symbol ki. However, the coordinate systems 
for the symbols k{ and k{^ are related by 

y-'\ for'' 



where O E 50(2, M) is the rotation in the plane spanned by t~'^ and u^'^ which 
maps i^'^ onto (and correspondingly u^'^ onto u^'^^). 

Similarly we introduce coordinates in 

V;,{k.i,ko) ^ R^xR', (p,g-)H^(^'-,r-) 



with 



I'- ■= |n'=-i''^o _ I _ ^fc-i,fco _ /'+ := l^'^i.'^o - 4J = cg. 

The next lemma shows that arbitrary pairs z~ , z~^ of points in the configu- 
ration space projections of the Poincare surfaces are connected by a trajectory. 



Lemma 10.4 Choose a large enough constant cq in Def. ( 19. J 6]) of the Poincare 
surfaces. Then for Eth = -E'th(eo) large, all E > Eth snd 

z^ e n^f^'^E) 

there exists a x~ G VE(fc-i, ^o, fci) of the form 

x~ = {y~,z~) with ■.= Ve{x~) = {y~^,z~^). 

Proof. 



By (I9.16P in the {y , z )-coordinates Ve(A;_i, ko) = 71% ^''^'^ is of the form 

VEik^i,ko) = ByX (10.15) 

with the two-disks By in momentum space and B^ in configuration space 

By := {y e R' I \y\ < \cpeo/E} D B, := {z eR^ \ \z\ < \cpe^/E]. 

(10.16) 
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We thus consider initial conditions of the form 
and their images 

x'- = {y'-,z'-) :=Ve{x-) e \4(fc-i,M- 

The idea of the proof is that, as By has twice the size of B^, the family of 
trajectories with initial conditions in By x {z^} has an opening angle 0{1/E) 
that is large enough to hit a full (9(l/£')-neighbourhood of c^g. Then by near 
collision, the opening angle is nearly amplified to 2tt, if eo is chosen large enough. 

The threshold energy Eth depends on the constant cq. However, by Lemma 
17.11 and a compactness argument one sees that the error term 0{1/E) in the 
C^-estimates 

and an identical estimate for the linearizations 




(10.19) 
(10.20) 



of these maps do not depend on cq, that is, C > may be chosen fixed and 
then Co chosen arbitrarily large. 
The relations 

p = V2E^ ■ {ts^ + yfi^ + yfu^ + 0{E-'^)) 
q = + 1^ [zff^ + zfu^) 

with E^ := E — V{m^) show that the angular momentum relative to the k^-st 
nucleus before and after scattering equals 

Lu, = {q- 4J xp = ±V2E^l^ ■ {{y^ + zt)u^ - {y^ + zt)i^ + 0{E-^)) . 
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Thus |Lfcol = O(eo/v^), and 

\Lko\>\l^CpeQ/^/E for {y',z') with y'edBy. 



Estimate (110. 17p shows that similar statements hold true for L^,, oV 



E ■ 



Finally (changing to (j?, g)-coordinates for a moment) by Prop. 18.6] we know 
that the image 



X 



(p-"+,g-"+) := V^+(x'-) G W'j,{k_,,ko,h) 



of x' G V^{k_i, kQ, ki) is near to the image 

of a point x'^^ = H(x'^) in a (9(l/_E')-neighbourhood of x'~ in the sense that 

-p-'+|/v^+ -gl+l = 0{l/E). (10.21) 
For the Kepler exit map ip'^i the formula 



e = e{x'£) = ^1 + 2^^2/^2 
for the eccentricity of the Kepler hyperbola parametrized by 



Z{1 + e cos((/9)) 



shows that the maximal total scattering angle 2arccos(l/e) G [0, tt) can be 
increased to tt — e (uniformly in E) by increasing cq, so that the diffeomorphism 

By x{z-}3x- ^ <f+(Sop-(x-)) G dBk,{cg) 

onto its image covers the sphere dBk^^(cq) around Syfc„ except an ^-neighbourhood 
of the forward direction s^o + Cg ■ E dBko{cq). 
The same holds true for the diffeomorphism 

By X {z-} 3x- ^ q^iVEix-)) G dBk,{cg) 

onto its image. So by the NC condition 12.41 an ^-neighbourhood U C dBk^icq) 
of the point n^^''"^' = Sk^ — CqS^ is contained in that image for eo large. By 
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Def. (110.141) of the ^'''''-coordinates this corresponds to a ^-neighbourhood of 
z'^ = 0. Estimate (I9.22p then implies that on this neighbourhood 



y'+iV'j, o V^ix-)) = z'+{V'e o V-^ix-)) + 0{l/E) {x' e B, x {-}). 
Thus we conclude from (110.181) that the intersection 

VE{By X {z-}) = V+oV'EoVE{By X {z-})nBy X +} c W{ko, h) (10.22) 
is non-empty. □ 

Remark 10.5 In fact the point x^ E VE{k-i, ko, ki) is unique, as will follow 
from Prop. 111.21 below. 

Lemma 10.6 For all admissible {k_i, ko, ki) E the Poincare return time 
equals 



Te \v,i,.,MM) = 7^7^ ln(i^ .J, ) + 0{E-'/'). 

Proof. 



with?/ := ipl^o^{T'''^''"''+{x),x). The estimate thus follows from (EHl), ( ISlB 
and (18181) . □ 



11 Existence of an Invariant Cone Field 

Definition 11.1 A cone at x E He is the image of the standard cone 

{{u,v) G X I \u\ > \v\} 
w.r.t. an invertible linear map 

X ^ T^He- 

A cone field C \n U C He associates to every x E U a cone C{x) at x. 
C is called invariant if for every x E U with Ve{x) E U 

T,Ve{C{x)) C C{Ve{x)), 

and strictly invariant if 

T.Ve{C{x)) C \nt{C{VE{x))) U {0}. 
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Alternatively one may think of cones as subsets of RP^. 

In the next proposition we use the Poincare section coordinates {y^,z^) 
defined in (fmSj) . 

Proposition 11.2 For C > large and E > E^ii the cone field C in 

U:= U VE{k.i,ko,h) 

admissible 

C(x) := I V^ij(A;_i,fco,fci) | \5y - Sz\ < ^\5y + 5z\^ 

(11.1) 

for X = {y,z) G U, is strictly VE-invariant, and the linearized Poincare map 
equals 

'1 1 



T,Ve = f{k-i,k^,k^)E- 
with 



11 11 



[xeVE{k_^,k,,k,)), (11.2) 



f{k^i,ko,ki) :-- 



(11.3) 



Remarks 11.3 1) The linearized Poincare map (111.20 is invertible, although the 
matrix ( J J) is not. 

2) Although the coordinates {y~^,z~^) on WE{k_i, ko, ki) C 7i^^'^^ do not co- 
incide on their common domain WE{k_i, /cq, ki) n Ve(^0; ^i, ^2) with the coor- 
dinates {y~,z~) on VE{ko,ki,k2) C they are related by 



Oy 
Oz 



where O E 0(2, M) is the reflection in the plane perpendicular to s''"''''^, trans- 
forming the unit vector t^{k-i,ko,ki) into ^^(^0,^1,^2) and the unit vector 
u^{k-i, ko, ki) into u~{ko, ki, ^2), see Def. (110. lift and (110. 12p . Thus the above 
cone field C is invariant under that transformation. 



Proof. 

By mm 



5z' 



5z 



1 
(1-^)11 i^l 



5y' 
bz 



by 
bz 



(11.4) 

OiljE). 
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From (18. 3p and (110.140 we deduce the transformation formulae 

= dv'" and 



{5v-,i-) 



(115) 



from the ,6z' )-coord i nates to the (5^ ,6w )-coordi nates. 
After near-collision we have 



{Sv+,u+) 



5y 



and 



6z 



(11.6) 



The reflection Rs-_s+ by the plane perpendicular to the vector s —s^ transforms 
the unit vectors as follows: 

= —Rs+-s-is~) , = +Rs+^s-ii^) and n+ = -|-i?^+„5-(n"). 

Moreover, up to an error term 0{1/E), we may substitute Rg+^s- for the reflec- 
tion Rij+-i;- in (EM . since by estimate the differences = 0{1/E) 
(and since |s+ — s~\ > 0). In particular 



Alp = 71- a{k_i, ko, ki) + 0{l/E). 



(11.7) 



Thus putting together (fTOD . (IlTSl) . (TO . (ICTl) . (fTTBl) and (fTTTl) . we 
obtain 



]1 ]1\ / ]1 

1 a Mo (/'-/c,)! 



11 



+ 



4Z-cos2(|a(fc_i,fco,/i;i))E / 1 ]1 



1 ]1 



+ 



5f 

5z- 



5z 



5z~ 



(11.8) 



0(E" 



as Cq + /' 



5z 



Inserting ( 11X8!) in (110.2011 yields 



11 

'11 ]1 
n 11 



f{k_i,ko,ki)E 



6z 



6z 



0{l/E) 

OiE^), 
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as /- = ^d^-^'''°, proving ( ITr2il . 

Strict invariance of the cone field C then follows from (lll.2p . □ 

We now describe the domains Ve(^) and images WE{k) of the iterated Poincare 
map with more precision. So let 

VEik){z) := VEik) n {By X {z}) {z G B,) (11.9) 

and similarly 

WE{k){z) := WEik) n {By X {z}) {z G 5,) 

consist of the points with the same configuration space coordinate z. We also 
consider them as subsets of By, forgetting the fixed coordinate z. It turns out 
that they are diffeomorphic to two-dimensional disks. We call 

dia.my{VE{k)) '■= sup diam(VE(A;)(i')), 

resp. 

diamyCWEik)) '■= sup diam{WE{k)i^)) 

the y -diameter of Ve(^). resp. WE{k) (which is measured with the Euclidean 
metric in the y coordinates). 

Corollary 11.4 For all z^ G 5^ the map 

VEik^i, ko, ki){z-) 3y- ^ f ") G B, (11.10) 

is a diffeomorphism. Thus the domain can be represented as 

VEik.,, ko, k,){z-) = {y-{S-,z+) I ^+ G fij, 



with 



D,y-{z-, z+) = -11 + 0{E-'), (11.11) 



and 



y-{z-,z+) = -z- + OiE-'), (11.13) 

t/je error terms being independent of the parameter cq in ( I9..Z6|) . 
/n particular the y-diameter of V^l^-i, ^o? ^i) of order (9(eo/-E^). 
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Proof. From formula (111.21) for the derivative of the Poincare map we see 
(observing that f{k_i,ko,ki) ^ 0) that (111.101) is a local diffeomorphism, and 
in fact a diffeomorphism onto its image. 

Thus by the definitions (110. 90 and (110.81) of domain and image (111.101) is a 
diffeo. Estimates (111.111) and (111.121) follow by inverting (fTT^ . Defs. (EH]), 
( flUTl) and (110.131) show that the y-diameter of WEiko, ki) is of order O{eo/E). 
Thus (I11.12P implies that the y-diameter of Ve(^-15 ^o? ^i) is smaller by one 
order. 

By definition (19.161) . the domain H^'' has size 0{eQ/E), so that directly we 
get only y^(z~ , z^) = —z~ + O{eo/E) instead of (111.131) . However, we then 
see from (111.111) that we may enlarge the size parameter cq without enlarging 
the error, and thus prove (111.111) . □ 

Proposition 11.5 There exist Ce > 1 and 6E > 0, such that for all E > Eth, 
m > 1, ^ G X^^ and {yt,zt) ■= VE{k){y~,z-) C Tc'^e the following 
holds true: 

1. For VEik){z~) defined in ( IJJ.9)) the maps 

VEm^') ^y'^ zt{y-,z-) G [z- G B,) (11.14) 

are diffeomorphisms. 

2. _ 

D,y,iz-, z+) = -11 + OiE-') (z^ G B,) (11.15) 

uniformly in m. 

3. For all G Bz the vector field on B^ given by z~ ^ y^(z~ , z~^) points 
inside the boundary: 

yl{z-,z+)-z- <0 (z-edBz). (11.16) 

4. VE{k){z^) contains a disk of radius 

m— 1 

Ce'E-' H {\fik_,, ki, A;,+i)l " (E + SE))"' , (11.17) 



i=l 

and is contained in a disk of radius 

771—1 



CeE-' n (1/(^^-1' ki, ■ [E - 5E))-' . (11.18) 



i=l 
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5. For k,lG X."^^^ with ki = li for i = 0, . . . ,m but k^+i 7^ Im+i 



m—l 



dist {VE{k)), VeH))) > C^'E-' II fc„ ki+,)\ ■ {E + SE))-' . 

(11.19) 

Remark 11.6 Using time reversal, one obtains similar estimates for WE{k). 
Proof. 

• m = 1: 

Then VE(fco, ^i) = By x by (110. ISp . and the radius of these two disks 
is proportional tp 1/E, see def. (110. 16p . So 

VE{ko, ki){z) = VE{ko, ki) n (By x {z}) = By x {z}, 

showing ( ITlTtI ) and ( imsl l. 

\f k = (/co, ki, k2) and Z = (fco, ^i, h) with /2 7^ ^2, then VE{k) and VE(i) 
are both contained in the Poincare section VE{ko,ki) = 7i^'''^\ but are 
mapped to different sections in the next iteration. These have a distance 
from each other that is bounded below by an energy-independent constant. 
On the other hand by (18.91) . scattering by the nucleus ki only leads to an 
expansion of order E, which conversely implies 

dist {VE{ko, ki, k2), VE{ko, ki, k)) > C^^E'K 

• m I— i> m + 1: 

By def. mm of Ve{K) 

VE{ko, . . . , km){z) = {By X {z}) n V~E\VE{k^, . . . , fcj). 
Thus the induction step for (111.170 and (111.180 is provided by 



i.e. formula (111.121) . 

The additional control of the i'-dependence of the sets Ve{]^ needed for 



the distance estimate (111.190 is guaranteed by the coupling (111.131) be- 



tween the y- and the z -coordinate. 
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One application of the linearized Poincare map (111.20 maps the cone field 
with cones 

{{Sy,Sz) G T^g,g)VEik_i, k^, ki) \5y -dz >{]] 

into the invariant cone field (111. ID . In particular this applies to vectors of 
the form {5z , 0). 

The image of the invariant cone field C under time reversal (using TR, 
defined in (110. SD ) meets the equation 

C 

\5y + 5z \ < —\6y - Sz\. 



This then implies the uniform estimate (111.150 . 
By (imi) 

diamy{VE{ko, ki, k2){0) = 0{1/E) independent of cq. 

Thus y~(0,z^) = 0{1/E) independent of cq, too. So the estimate 
(111.151) of the first partial derivative gives. 

yl{z-,z+)=yl{0,z+)+ [ D^yl{tz-,z+)dt = -z- + 0{l/E). 

Jo 

Now the modulus of z~^ E dB^ equals \cpeQ/ E (see def. (110.161) ). or 

y , z+)-z= lcpeo(-|cpeo + 0{1))IE^ [z' e dB,). 
Enlarging the parameter cq if necessary gives (111.161) . 

To prove the iteration step for the first statement, we first show that that 
map (111.141) is onto: 

For all z,z^ E B^ there exists an y^ with z'^{y~ , z~) = z. (11.20) 

Namely for z^ E B^, k^ := {ko, ki, /C2) and k^^ := {ki, . . . , km) we con- 
sider the vector field 

F-.B^^R' , F{z):=yTn{z,z^)-yUz-,z). 
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— * 

If F has a zero, there exists an initial point {y~ , z~) G Ve{^{z~) meeting 
(111.201) . However, the vector field z ^ y^i,{z ,z^) points inside the 

boundary of B^, whereas by time reversal 

>0 (zedB,). 

Thus F points inwards, too, so that the degree of the vector field F on 
the disk is non-zero, which in turn implies that F{z) = for some z (see, 
e.g., Hirsch [H]], Chapter 5). 

Inspection of the derivative (111.2ft of the Poincare map shows that y 



z~j^{y ,z )\s injective, smooth and smoothly invertible. □ 

Lemma 11.7 For E > Eth, / < < r and k = {ki, . . . ,kr) admissible 

diam{WEiki, . . . , fci) n Vb(A;o, • • • , K)) < 

4 {disimy{WE{ki, . . . , ki)) + diamy{VE{ko, . . ., K))) . 

Proof. Set / := {ku ■ ■ ■ ,ki) and m := (fco, . . . ,kr). Let = {y\ z^) and 

A^^ = {y^^,z^^) be two points in WeH) n VE{m) and 

the i'-components of 'P£;(m)(A^) resp. VEiD^^iA^^). The disks 
D' := {x e VEim) n We{V) \ z^{x) = z''^} 

and 

D" := {x e VEim) f] We{V) \ z;{x) = z^^'~] 

intersect in a (unique) point A^^^ = {y^^\z^^^) of n''j^'^\ 

As \\A^^ - A^W < \\A^^ - A^^^W + \\A^^^ - A^l the lemma follows from the 
estimates 

diam(L'^) < Adiamy{WEiL)) and diam{D^^) < 4diamy(yE(m)), 

and by time reversal symmetry (110.51) it suffices to prove the first one. 
is mapped diffeomorphically onto We{L){z^) by 

X t-^ {y', z^) with the unique y' G By meeting zj{y\ z^) = z'[{x) 

(uniqueness of y' follows from Prop. 111.51 1). By the estimate (111.2ft on the 
linearized Poincare map this map increases distances by a factor less than 2. □ 
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12 Symbolic Dynamics 

In this section we analyze the set 

bE = bnEE 

of bounded $*-orbits of energy E. 

To that aim we equip the alphabet S = {1, . . . ,n} with the discrete topology, 
and for an interval (see (110.61) ) / = ![ C Z the space with the product 
topology. Finally, we introduce the topological subspace 

X[ := G 5^ I admissible} (12.1) 

of admissible sequences (that is, /cj+i 7^ ki), and use the abbreviations 

X:=X-^ , X+:=X- and := Xi^. 

The space X of bi-infinite admissible sequences is empty for n = 1, consists 
of two points for n = 2, and is a Cantor set for n >3. From now on we assume 
n>2. 

The shift 

(T : X ^ X , cr{k)i := ki^i {i G Z) 

is a homeomorphism on X. It is well-known that the topology on X is generated 
by the metric 

d{k,V) := 5^2-1*1 ■ (1 - e X). (12.2) 

For an admissible sequence = (k^, ki, . . .) G X+ we define 

VE{k+) := fl VE{ko, . . . , U C VEiko, k,). (12.3) 

mGN 

Similarly, for an admissible sequence k~ = {. . . , ko, ki) G X^ we define 

WE{k-) := fl WE{k.^, . . . , fci) C WEiko, k,). (12.4) 

meN 
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Lemma 12.1 For E > E^^ and keX+ the sets VE{k;^) and WE{k~) are the 
graphs of functions 

VE{k~^) : Bz ^ By resp. we{^) : B^ ^ By 

meeting the Lipschitz estimates 



E 



(12.5) 



resp. 

WE{k~){zi) - WE{k~){z2) - {zi - Z2) \ < C^-^-^—^^ {zi,Z2 G Bz 



for some C = C(i?th)- 

Proof. The sets in (112.31) are nested: Vsiko, . . . , km^) C VE{ko, . . . , fcmi) for 
m2 > mi. By estimate (111.181) . the y-diameter of these sets goes to zero as 

m — >• cx). 

We set / := (fco, • • • , km). For i*^ G dB^, the z -dependence of the curves 
{{yY{z,z'^),z) I z E Bz} in the boundary of Ve{1) is controlled by the m- 
uniform estimate (111.151) . showing (112.50 . Finally, the statements concerning 
WE{k~) follow by time reversal (110.50 . since 

WE{kr) = Tn{VE{ki, ko, ...)). □ 

Let 

A|:= U ^^(^'■) ■ ^E-= [j ^Eik-) (12.6) 
fe+ex+ fc-ex- 

and 

Ae := A^nA^. 

By restriction, we associate to A; G X half-infinite admissible sequences G 
X^*". Then we define a map J-'e '■ X — > A^; by 

J^E{k):=VE{k:')nWE{k-) (A; G X). (12.7) 

Note that in view of Prop.[II3]the disks yE(fc+) and WE{k-) in T^^^'^^ intersect. 
The Lipschitz estimates of Lemma [12.11 imply that their intersection consists of 
precisely one point, which we identify with an element of A^;. 
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On X we introduced in (112.21) the metric d. 

On the Poincare surfaces He = Ufc^i^s' use the metric 



dist(x, x' 



dist : He^'He 

1 , otherwise 



based on the , i')-coordinates (110.131) of a; and x' (Remark [10.31 showing that 
dist is well-defined). 

Lemma 12.2 There exist a > such that for E > Eth, J^e is an (a ■ InE)- 
Hdlder continuous homeomorphism, that is, 

dist{J'Eik),J^E{L)) < C{E) ■ {k,l e X), (12.8) 

for some function C > of the energy, conjugating the shift with the restricted 
Poincare map Ve ■= Ve \ae- 

J^j,oa = V^oJ^E. (12.9) 

Proof. Since C may depend on E, we can assume without loss of generality 
that the central blocks of k and / coincide, i.e. (A;_i, A;o, ki) = /q, h). Let 
I[ C Z be the maximal interval containing on which k and / coincide. Then 
by Lemma ULT] and ( ITLlSl) 

dist(^i,(A;),^,,(/)) 

< diam{WE{ki, ...,ki)n Vb(/co, • • • , K)) 

< 4 {diamy{WE{ki, . . . , fci)) + diamy{VE{ko, • • • , kr))) 

with /min := miui^j^k \ fihj,k)\, setting a slightly smaller than l/ln(2). □ 

For n > 2 we denote Te o Te '■ X IR+, with the return time Te defined in 
(110. 3p by Te, too. Being defined by composition of a smooth map with a Holder 
continuous map, Te is Holder. 

The continuous flow on Be is modelled as follows. 



85 



Definition 12.3 Given a roof function r E C°(X,]R+), we set 

:= X X R/ ~ 
where ~ is the equivalence relation defined by 

{k,t + r{k)) {a{k),t) ((^,t)GXxM). 

Then the r-suspension flow is given by 

a*:Xz,^Xs, [ik,s)]^[ik,s + t)] (t G M). (12.10) 

In the interesting case r = Te vje abbreviate X^ := X^^ and cr^ := cr^^. 

Definition 12.4 The set CT of collision times of a trajectory c : / — P is given 
by 

Cr{c) ■={tel \ qo c{t) G {si, . . . , Sn}} . 
The Morse index of a hyperbolic T-periodic trajectory c : [0, T) — > P is given by 

Index(c) := ^ dim (P„(c(t)) H Vertc(t)) + |Cr(c)|, 
te[o,T)\cr{c) 

where Es/u{c{t)) C Tc(t)P denotes the weak (un)stable subspace at c{t) and the 
vertical subspace Vert^^ C T^P at a; G T*M C P is the one annihilated by the 
linearized configuration space projection T^q. 

Remarks 12.5 1) The weak (un)stable subspaces Es/y_{c{t)) are the direct sums 
of the (2-dim.) strong (un)stable subspaces and the neutrally stable flow direc- 
tion (see |KH] ). Similar to Vertc(t), they are thus 3-dim. Lagrangian subspaces 

of Te(i)P. 

2) Actually Index(c) is always finite. We already know from our study of the 
(near)-collision process that CT{c) does not have accumulation points. On the 
other hand, the set of times t where the (un)stable subspace Es/u{c{t)) turns 
vertical is finite, too, since the kinetic energy term ^p"^ in the Hamiltonian func- 
tion is a positive quadratic form, see Duistermaat [Du]. 

3) The additional term |Cr(c)| in the definition of the Morse index is chosen in 
the only way that makes that definition invariant under small perturbations of 
the flow. 
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For < 6 < TV we introduced in (19.11) the pericentric hypersurfaces 

= e H, I k» - ?,| < '!;1(,;4) } (' = 1, . . . , «) 

(with Hi defined in (18. ID ) near the Ith nucleus. So the angle parameter 9 fixes 
the precise meaning of the term 'near-collision'. 

Definition 12.6 The set J\fCTg(c) of 6 -near-collision times of a trajectory c : 
/ — i> P is given by 

AfCTeic) ■.= {tel\ c(t) G u::=i7^K^)} • 

We say that the trajectory c -visits the nuclei k E X[ /n succession \f J\fCTg{c) ^ 
0, 

r = |ArCTe(c) n [0,oo)| , I = 1 - {NiOTe (c) f] {-oo,0)\ 

and 

for the enumeration J\fCTo{c) = {ti \ i G !'[, ti < tj+i, to < < ti} of 
near-collision times. 

Remarks 12.7 1) k is really well-defined, since the hypersurfaces 'Hi{9) C Hi 
do not intersect, and since the flow $* is transversal to Hi so that J\fCTg{c) is 
discrete. 

2) As long as we analyze bounded motion, the natural angle will be 6* = amin, 
so in that case we do not write it explicitly. 

Lateron, in the analysis of the scattering process, 9 will be a parameter, since 
scattering with small angles shows non-universal features. 

Theorem 12.8 1. For an NC configuration and E > Eth the map 

J^I-.Xe^Be, [{k,s)]^<^'{J^Em 

is a Holder continuous homeomorphism conjugating the suspension flow 
with the flow on the set Be of energy E bound states: 

oJ=\ = T\oa'E (tGM). 

Thus for 
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• n = 1 there are no bounded orbits, 

• n = 2 that set consists of one closed orbit: Be = 

• n > 3 is locally homeomorphic to the product of a Cantor set 
and an interval. 

2. All bounded orbits are hyperbolic. 

3. bE has measure zero w.r.t. Liouville measure \e on T,e. 

4. If a T -periodic trajectory c : [0, T) ^ S^; visits the nuclei (fci, . . . , km) G 

in succession, then its Morse index equals Index(c) = m. 

Proof. 1) By the virial identity (I2.12p all trajectories that leave the interaction 
zone XZ go to spatial infinity and thus do not belong to bE- Hence 

bE C Ve with Ve C T.E defined in ( ITO) . 

Lemma fTO.ll tells us that for bE the Poincare return time Te is uniformly bounded 
by 0{1/VE), so that a fortiori 

^TE{x),x)enE (xebE). 

Thus 

bE = HR,AE), 

showing that !F'^ : X^; — > fe^; is a bijection. Using Lemma 112.21 we see that 
is a Holder continuous homeomorphism, since the flow ^\Tie on the energy 
shell can locally be straightened out (see, e.g., [AM] , Thm. 2.1.9), and since the 
Poincare section He is transversal to the flow. 

Now for = 1 the space X of admissible sequences is empty, and for n = 2 

X= {(...12121... ),(... 21212...)} so that Xe = S\ 

Finally, for n > 3 the sequence space X is a Cantor set (a non-void compact 
totally disconnected set without isolated points, see, e.g. Katok and Hasselblatt 
[KH], A 1.) 

2) The existence of the strictly invariant cone field C (defined in (111. 1ft ) on A^; 
implies the hyperbolicity of the flow. The expanding (as well as the contracting) 
subspace of Tj-S^;, x E Ae, is two-dimensional, since it equals 



fl DP^CiV^^ix)). 
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3) The property of 6^; C S^; to have measure zero is defined without reference 
to a measure (since in every local chart of the smooth manifold T,e the image 
of Liouville measure A^; is continuous w.r.t. Lebesgue measure, see Hirsch [Hij . 
Chapt. 3.1). In fact, applying the Straightening Out Theorem ([AMJ, Thm. 
2.1.9) it suffices to show that A^; has measure zero. This follows for E > Et^ 
from estimate (112.81) . 

4) Setting ko := km, we denote by Si the intersection times with the Poincare 
surfaces, that is, c{si) E 'H^"^' ' {i = 1, . . . ,m). 

We assume w.l.o.g. that the T-periodic trajectory c : [0,T) T.e begins 
and ends in the Poincare surface 7i^°''^\ so that we have Si < ti < Sj+i and 
Si = 0. Additionally we set s^+i := T. Then we prove that Index(c) = m by 
showing that 

dim(E„(c(t)) n Verte(t)) + |Cr(c)n[si,Si+i)| = 1 (z = l,...,m), 

te[s„si+i)\cr(c) 

which means that every near-collision adds one to the Morse index. This claim 
is equivalent to 

(12.11) 



t&[si,Si+l)\{ti} 



where we used the coordinates c{t) = {p{t),q(t)). 

First we treat the case Zk^ < of a repelling Coulomb singularity. There we 
know that the trajectory does not touch the singularity (g(ti) ^ SfcJ. 

The simplest case is the one where along a neighbourhood of Sj+i]) the 
potential V coincides with the Keplerian potential — — s^. |, and the angular 
momentum L^. = {q{t) — s^J x p(t) of the trajectory relative to the nucleus 
vanishes. 

Then p(ti) = 0, which gives a contribution of one in (112.110 . since there the 
flow direction, given by the Hamiltonian vector field 



XniciU)) = i-VV{q{t,)),p{t,)) e Vert 



is vertical. 

However in that case there are no further contributions to the index. This may 
be seen as follows. The intersection of the (three-dimensional) weak unstable 
subspace at c(sj) with the tangent space to the local Poincare section T-C^ ''^^ 
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is a two-dim. subspace which lies inside the local cone field C(c(sj)) defined in 
(lll.ip . Thus we have dq-Sp>0 for a variation vector 

{Sp, Sq) e Eu{c{si)) n T,^s,)V!'-^'^' C Te(,,)P, (12.12) 
and this property is preserved by the forward flow, since 

j^mt) ■ m)) = -m) ■ D^vmmt) + m) ■ m) > o. 

For general repelling potentials and general trajectories c (112.111) equals one, too, 
by a continuity argument based on estimate (111.20 . 

Now we treat the case Z^- > 0, again starting with a Keplerian potential 
—Zi/\q — and vanishing angular momentum 1/^.. Here q(ti) = s^., so 
that we have to show that Eu{c(t)) does not turn vertical for t G [sj,Sj+i) \ 
{ti}. However, the configuration space trajectories g([sj,tj]) and (f([tj, Sj+i]) 
are straight lines, and the infinitesimal two-parameter family of diverging Kepler 
hyperbolae corresponding to initial conditions (112.121) do not have a conjugate 
point. A more formal way of seeing this is to use the Jacobi metric gE{q) = 
(1 — V{q) / E)g(q) discussed in (A Jl.lD of the Appendix. Then the energy E 
solution curves correspond to geodesies c in that metric, and the linearization of 
the flow $f corresponds to the Jacobi equation 

V^Y{t)+R,^t)Y{t)=6, 

see, e.g., fKl]. For a variational vector field Y orthonormal to the geodesic 
velocity vector c{t) the self-adjoint curvature operator Rc(t) has an ^-expectation 
value equal to the sectional curvature of the plane spanned by these two vectors 
([Ki], Prop. 1.11.3). The sectional curvature of has been calculated in (A J1.2p 
of the Appendix. 

For the geodesic under consideration the relevant plane contains the direction 
q — Sk of the singularity and is thus seen to be negative definite (setting qs = 
in (AOD). 

However, negative sectional curvature, together with initial conditions (112.121) . 
lead to absence of conjugate points. This proves (112.111) for the collision orbit. 

In the case of non-vanishing angular momentum 1/^.. the sectional curvature 
(A ll.2p ) is still negative in the plane perpendicular to L/,. . However, near the 
singularity (A J1.2p ) becomes positive for plane containing the vector Lfc. , leading 
to a conjugate point. So (112. lip holds for that case, too. 

Like in the case Z^. < 0, a continuity argument based on estimate (lll.2p 
shows assertion 4. for general potentials V and large E. □ 
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13 Fractal Dimension 



In this section we estimate the fractional dimension of the set Be of energy E 
bound states, for E large. This quantity, being of interest in its own right, 
governs the measure of those scattering orbits which have a large time delay. 



Dimensions: Definitions and Elementary Properties 

Besides the well-known dimension dim-;^ introduced by Hausdorff and Besicov- 
itch, there exist several other definitions of the fractional dimension. Of those 
we will only consider upper box-counting dimension dimg, since most dimensions 
take values between diiiiT^ and dimg, see Falconer IFal . 



Definition 13.1 Let {X,d) be a separable metric space and U C X , U ^ 
The diameter diam(f/) ofU is given by 

diam(f/) := sup{d{x,y) \x,y eU} and diam(0) := 0. 

For ^ C X, s > and (5 > let 

{oo 
^(diam(t/i))' | E C Ui^nUi,disim{Ui) < 6 

(setting 0° := 1 except for the case diam(0)° := 0). The HausdorfF s-dimensional 
outer measure of E, Ti^^E), is then defined by 

H\E) := suvUm = \^ni{E)- (13.1) 

<5>0 "5X0 

The Hausdorff dimension of E is given by 

<lim-u{E) := sup{s G M | W{E) = oo} = sup{s G M | W{E) > 0}. (13.2) 



The important property of is that it defines a measure on, say, the Borel sets. 
Therefore, Hausdorff dimension is not only monotone, that is 

EiC E2^ dim7^(Ei) < dim7^(E2), (13.3) 

but also countably stable: 

dim-H{Ui(=nEi) = sup dimniEi). (13-4) 
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Definition 13.2 Let {X,d) be a manifold X with metric d and E C X be a 
non-empty bounded subset. The (upper) box-counting dimension (or Minkowski 
dimension) dime(-E') is given by 

HNsiE)) 

dimg(i?) := limsup , 

£\o - \n{€) 

where N^iE) is the minimal number of balls of radius e needed to cover E. 

dims is monotone (see (113.31 )). and 

diiRB^Ei U E2) = max(dimB(Ei), dimf5(^2)), (13.5) 

but it is not countably stable, since dim^ is invariant under closure. Furthermore, 
for all (bounded) E, 

dimniE) < dimB(E). (13.6) 

Example 13.3 The triadic Cantor set c M has 

In 2 

dimH(^) = dimB(E) = 7-^ = 0.6309 ■ ■ ■ . 

Ill O 

On the other hand, the set Q C M of rational numbers has dimensions 

dim„(Q) = , dimelQn [0,1]) = 1. 

We have 

dimn{f{E))< dimniE) and dims{f{E)) < dimniE) (13.7) 

if the map / is Lipschitz. In general both dimensions are not additive w.r.t. 
cartesian products, but the inequalities 

dimn{Ei) + dimn{E2) < dimn{Ei x E2) < dimn{Ei) + dimf5(£^2) (13.8) 

and 

dims(Ei X E2) < dimB(^i) + dimB{E2) (13.9) 
hold true, see Mattila [Ma] and Tricot [Tr]. 
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Lower Estimate for the HausdorfF Dimension 

By (113.61) we need a lower estimate for the Hausdorff dimension and an upper 
estimate for the box counting dimension. 

Whereas any choice of coverings by e-balls leads to an upper estimate for 
dimg, a lower estimate for diniT^ is provided by the mass distribution principle: 

Proposition (see [FaJ). Let ;U be a probability measure on {X,d) and n{E) > 
0. Suppose that for some s > 0, C > and S > 

I2{U) < C ■ diam(f/)^ (diam(?7) < S), 

Then dim7^(_E') > s. 
Proof. If ^ C Ui(znUi, then 

oo oo 

i=l 1=1 

SO that < /i(-E) < Cl-Ll{E). Then the statement follows from the second 

expression for diniT^ in (113.21) . □ 

Thermodynamic Formalism for Dimension Estimates 

We base ourselves on the size estimate Prop. Ill.Sl for Ve{M) in terms of the finite 
geometric data encoded in / (see (111.31) ). 

Hence we consider the parameter-dependent {n{n — 1) x n{n — l))-matrix 
M(s) with double-indices in j) G 5 x 5 | i 7^ j} and entries 



M{s] 



i,j;k,l 



\f{i,k,l)\ ' ,i ^ j^k I and j = k 
, otherwise 



For s G M this is a matrix with non-negative entries, and for n > 3 centres all 
entries of (M(s))'" are strictly positive iff m > 3. Thus by the Perron-Frobenius 
(PF) Theorem M{s) has a unique eigenvalue Amax(s) of largest modulus, which 
is positive and of multiplicity one, and the corresponding eigenvector v(s) can 
be chosen to have strictly positive entries. 

In the case n = 2 the two eigenvalues of M(s) are given by ±Amax(s) = 

±l/(l,2,l)r^ 
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Lemma 13.4 For all E > E^h there is a unique solution d{E) := s of the 
equation 

A^ax(s) = E\ (13.10) 
and d{E) = for n = 2, whereas for n >3 



Finally 



Proof. 

We conjugate M{s)/ E^ with the {n{n — 1) x n{n — l))-matrix D{s) given by 

u{su,,.k,i - <^ Q , otherwise 
Then M{s,E) := D{s)M{s)D{sy'^ / E" has the non-negative entries 

M{s,E),,,,,-^ , otherwise ^^^'^^^ 

with 

fit, k, I) := v/#V^|/(^, kj)\>0 (t^k^l). 
The right Perron-Frobenius eigenvector of M(s,E) equals 

v{s) := D{s)v{s) with eigenvalue \^^^{s,E) := Amax(s)/-E'^ 

and (113. lOp corresponds to the implicit equation 

K.ME),E) = 1. (13.14) 

Although /(/, k, i) = f{i, k, I) by def. (fTOI) of /, for n > 2 the matrix M(s, 
is non-symmetric and even non-normal, leading to slightly more complicated 
estimates. 

We denote the left PF eigenvector by w{s), again assuming positivity of its 
entries. Then for E large \ma,xis,E) is strictly decreasing in s, since 

w{s),DiM{s,E)v{s] 

DiA„,ax(s, E) = ^ < 0, (13.15) 

{w{s),v{s)) 
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Di denoting the derivative w.r.t. the ith argument. 

Inequality (113.151) follows from the fact that all components of v(s) and w{s) 
are positive, and for E large all entries of DiM(s, E) are non-positive, and some 
are negative. So the l.h.s. of the equation Amax(s,-E) = 1 is strictly decreasing 
in s. Furthermore, 

Amax(0, E) = n-l>l and lim Amax(s, E) = 0. 

This implies that ( 113. lOD has a unique solution d{E) = forn = 2 and d{E) > 

for n > 3. 

More precisely we observe that the n — 1 non-vanishing entries in each row 
(ij) of M(s,E) are of the form (fE)-'. Setting 



/inin := mill /) and /, 



max 

:= max /(i, fc, /), 



we note that by a consideration of the largest and the smallest components of 
the eigenvalue equation Amax(s, = M{s,E)v{s) 



Hn - 1) ^ ^(^^ ^ Hn - 1) 



ln(/max-E) 



ln(/min-E) 



(13.16) 



This shows the estimate (113.111) . The implicit equation (113.141) for d{E) shows 
that 

-D2~Xrr,ME),E) d{E) 



d'(E) -- 

D,X^,ME),E) E-D,\^^^{d{E),E) 
So we have to estimate 



(13.17) 



D^\^,,{d{E),E) 



w{s),DiM{s,E)v{s] 



{w{s),v{s)) 



s=d{E) 



(13.18) 



more precisely in order to show (113. 12p . The matrix 



M := lim M{d{E),E) 

E^oo 



has the form 



M 



— 1) ,i ^ j, k ^ I and j = k 
, otherwise 
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since by (113.161) the quotients 

d(E) 



MidiE),E),,k;k,l/M{diE),E)i,,k';k',l' = [f{l\k\l')/f{l,k,l)^ 

= l + C(l/ln(E)) (13.19) 

of the non-zero coefficients converge to one, and since the PF eigenvalue of 
M{d{E),E) equals one. We now use the algebraic relation 

{{n - l)My + {n - 1)M = J" with J^i,r,k,i ■■= I {i^],k^l), 
which is approximately met for finite energies in the sense that by (113.19^ 



R{E) ■= (^{n~l)M{d{E),E)^ + {n - 1) M {d{E) , E) 

= J^ + C(l/ln(E)). (13.20) 

By (113.141) the PF eigenvectors fulfil the equations 

R{E)d{d{E)) = n{n - l)v{d{E)) , R{Efw{d{E)) = n{n - l)w{d{E)) 
so that in view of estimate (113. 20p its components are nearly equal: 



y,,.^,v{diE))ki 

mm, = ' ■ (1 + o{i/HE))) ^ j) (13.21) 

and similarly for w. 

Finally, for i ^ j,k ^ I and j = k 

D,M{d{E), E\,,u,i = - ln(|/(z, k, l)E\)/{n - 1) ■ (1 + 0(1/ HE))). 
Inserting that estimate and (113. 2ip into (113.180 we obtain 



DiA^ax(rf(^), E) = - \n{E) + 0{l). 
Putting that estimate and (113.161) into (113.171) gives (113.121) . □ 

Dimensions of the set of Bound States 

We now estimate dime(6£;) and dim7^(6£;) by the solution d{E) of the matrix 
eigenvalue problem (113.101) . 
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Theorem 13.5 For E > Et^ and n > 3 the HausdorfF dimension diiiiT^ and the 
upper box-counting dimension diiiig of the energy E bound states Be meet the 
estimates 

1 + 2d{E) ■ (1 - 0{{E\nE)-^)) < dimn{bE) < 

< dimisibE) <l + 2d{E) ■ {l + 0{{E\nE)-^)) (13.22) 

with the solution d{E) of lilS.lOj) . In particular they meet the rough estimate 

dimnibE) = 1 + ^^^^^ ' + O ((lnE)-2) = dimB(fe^) (13.23) 
For n = 2 centres dim-^(6£') = dimg(6£;) = 1. 

Proof. The second inequality in (113.220 is the abstract inequality (113.60 . Esti- 
mate (113.2311 follows from (113.221) by inserting (113.1111 . 

So it remains to prove the lower bound for dim^ and then the upper bound 
for dimg. For n = 2 centres bs consists of one closed orbit, whose dimensions 
equal one. So we assume from now on n > 3. 

1) Using the constant 6E from Prop. 111.51 and setting El := E + 6E for 
E > Eth, estimate (113.121) shows that 

d{E) =dL-{l-0{l/{E\nE))) for di := d{EL), 

so that the first inequality in (113.221) follows from an estimate 

dimnibE) >l + 2dL {E > E,^). (13.24) 

We show that this follows from 

dim^(As) > 2dL [E > Eth). (13.25) 

First we find e > such that the flow $ : M x P ^ P, restricted to (— e, e) x He, 
is a diffeomorphism onto its image. Thus 

U, := X Ae) C bE. 

satisfies 

dimniUs) = dimnii-e,e) x Ae) 

< dimn{{-e,e)) + dimn{AE) = l + dimn{AE) (13.26) 
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using (113.71) and (113.81) . By finiteness of the Poincare return time (Lemma 
110. ip . can be covered by finitely many time translates 

J max 

bE=[j^ej,U,) with j^,^ ■.= [e-'snp Te{x)] (13.27) 

so that 

dimnibE) = maxdimn{^{ej,Ue)) = dimniUe) 
j 

follows from stability of diniT^ and (113.71) . So we are reduced to show (113.250 . by 
employing the mass distribution principle. This will be based on a P^-invariant 
probability measure he on A^;, which is the image 

w.r.t. (112.71) of a measure /ix.E on X. yUx.s is defined through its values 



-dL 



(13.28) 



1=1+1 

on the cylinder sets 

Z{k) := {fc' G X I A;^ = ki^i G {/,... ,r}} {I <r eZ, keX^ 



with the left Perron-Frobenius eigenvector w^di) of M{dL, E^). 
Using the /^-normalization 

Y,Hs)^,k = l (13.29) 

of w{s) and the relations (113. 13p and (113. 14p we see that the definitions (113. 28^ 
are compatible and define a a-invariant Borel probability measure /ix,£;- By the 
conjugacy (I12.9p between the shift a and the restricted Poincare map the image 
measure he is then indeed P^-invariant. 

We now claim that the mass distribution principle applies. Namely we have 
for some C = C{E) > and 5 = 6{E) > 

I^e{U) < C ■ diam(f/)2^^ (13.30) 
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for all measurable U C of small diameter diam([/) < S. Instead of general 
such U we first consider balls B(S') C He of small radius 6' > 0, which are 
centered at a point x E Ag. The symbol sequence k := J^^^{x) G X of this 
point projects to the half-infinite sequences A;^ G (see (112.70 ). There are 
unique integers / < < r with 



Q'^l' n > (13.31) 

i=l+l 


but C,'E^' llifih^u h+i)EL)-^ < 6'. (13.32) 



respectively 



r-1 



Ce'El' llifih.i, h, h+i)EL)-^ > 5' (13.33) 



i=l 



but C,'E^' Hifih.,, h+i)ELy^ < 5', (13.34) 

i=l 

By the lower estimates (113.311) and (113.331) and (111.191) of Prop. \TTS\ 

VE{K,...,K)nB{6')y^il} only if {k'„ . . . , K) = {ko, . . . , K). 
and similarly for WE{k'i, . . . ,k[). However, this implies that 

B{5')nAECj'EiZ{ki,...,kr)). 

Thus by def. ( 113.281) of fix,E and the upper estimates (113.320 and (113.340 

MB{6')) < fIx,E{Z{ki,...,kr)) 

r-1 ^ 

= w{dL)kuki+i- n (^f{ki-i,ki,ki+i)EL^ 
1=1+1 

< c„ ■ (25' f''^ = ■ diam(S(5'))'''^ 

the constant 



Cn 



being independent of 5' and A; (remembering that /max = maxj^fc^z /(^, A;, /)). 
That is, the balls 5(5') satisfy the mass distribution estimate ( 113. 3O0 . 
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General sets U C He of diameter ^6' and f/ fl A^; 7^ are subsets of such 
balls of radius 6' centered at a point in U DAe- This proves the mass distribution 
principle ( 113. 30p in general. 

We see from (113. SOD that for any cover A^; C Ujf/j with diam(f/j) < 5, 



i i 

which implies H'^'^^^Ke) > 1/C and thus dim7^(A£;) > 2^^. 

2) Similar to the first case, to obtain the upper bound in ( 113. 22p for dimts^bE), 
it suffices to show that for 6E > from Prop. 111.51 

dimtsibE) < I + 2du (13.35) 

with 

du := d{E - 25E) and Eu := E - 5E. (13.36) 
Also ( 113. 35p follows from the estimate 

dimelAij) < 2du, (13.37) 

for then dimg(?7e) < l+2du (for the set f/^ = $((—£, e) x A^;)) is a consequence 
of (113. 7p and (113.91) . Relation (113.50 may then be used to determine the box 
dimension of the covering ( 113.270 . 

We now claim that for e > small we can cover A^; in the form 

Ae C U J^E^Zi^k)) (13.38) 



with index set 



l=-lr=l 



^ > C,E-' n {\f{h-i,h,h^,)\-Eu)-' > ^ . 

z=t-\-l 
r-1 

^ > CeE-'l[{\f{k.,,h,k.^,)\ ■ Eu)-' > ^ 

16 le/max-Ec 
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• To prove (113.381) . we choose an arbitrary bi-infinite sequence / G X, and 
show that there \s a ke with / G Z{k). 

Since we have chosen a small e, the right inequalities in the definition of 
are met for / = —1 and r = 1. We find k by setting ki := U and choosing 
r such that the right inequality would be violated for r + 1. Then the left 
inequality holds. Similarly we choose the minimal possible / < 0. 

• Next we claim that for kEle the subset J^E{Z{k)) of the Poincare section 
is contained in a ball of radius e. This holds true since 



• Finally, denoting by h := l/(16/max-E'(7) the constant appearing in the 
definition of le, the number Ns{Ae) of e-balls needed to cover is 



= WEih, ...,k^)n VEiko, ...,kr), 



see ( lT?7i) . so that using Lemma [HT] and (111.181) 



diam(JP-£;(Z(fc))) 
< 4 {diamyiWEik, ki)) + diamy(VE(fco, . . . , kr))) 
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bounded by 



Ne{AE) < \Ie\ 



< 



'OO oo 

ke U 



Z=-l r=l 



i=l+l 

r-1 



< 



u u^^ 



-1 r=l 



«=1 
r-1 



1=1+1 

— oo oo r— 1 

Z=-l r=l fceX[ i=Z+l 
oo 



u 



u 



with the matrix M{s,E) from ( 113.1311 . and 1 e M"("-i) denoting the 
vector whose components equal 1. Now for the arguments (113. 36p the 
Perron-Frobenius eigenvalue Xma.x{s,E) of M{s,E) is smaller than one: 

r , , ^ , XraUdu) XmUd{E-2SE)) 
Xra^^[du,Eu) - 

\^Ud{E-25E)) ^ 
^ {E - 26EYU - ^^3.39) 

In view of the near-constancy (113. 2ip of the left PF eigenvector v = 

v{d{Eu)) of M = M{du,Eu) we use the inequality 

1,M"]1) < 2n(n- 1)^- (men), 
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to exchange the vector 1 by v 

'v,M{du,EuTv 



oo 

m 

m=l * ' 

oo 



= {ehE/C^)-^''"2n{n - 1) ^ m ■ X^^du, 
< C'e 



m 

U) 



m=l 



with finite C", using (113.390 . Inserting this into the Def. 113.21 of the box- 



counting dimension shows (113.370 . □ 



14 Topological Entropy 

We shall now determine the topological entropy of the flow = ^^ts^ on the 
energy shell T^e- 

The estimate of htopi^E) based on Prop. 114.31 stating that the topological 
entropy of the flow is determined by its restriction to the set Be of bound states. 
It is then relatively easy to compute that topological entropy using symbolic 
dynamics. 

First we formally introduce the notion of topological entropy of T : X X, 
following the definition of Bowen. That definition is in a way more general than 
others in allowing for non-compact metric spaces {X,d) (see Walters |Wa] ). 

Definition 14.1 Let (X, rf) be a metric space and T : X ^ X a uniformly 
continuous map. Then for m G N, £ > 0, a subset F C X is said to {m, e)-span 
a compact C X if 

\/x e K^y e F : dm{x,y) < e, 

with the metric dm{x,y) := maxo<i<m-i '^(T'^x, T*?/). 

Let rm{£,T, K) be the smallest cardinality of an (m, e)-spanning set F of 

K, 

hr{e, T, K) := limsup — ln(rm(£:, T, K)), 

m— »oo Tfl 

and h,,p{T,K) := lim,^ohr{e,T,K). 
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Then the topological entropy ofT is 

/itop(T) := sup {htopiT, K)\K CX compact} . 



Note that hr{e, T, K) is monotonically increasing in e, so that its limit htop(T, K) 
exists (it can be oo). 

If d! is a second metric on X uniformly equivalent to d, that is, if the maps 

Id : (X, d) {X, d') and Id : (X, d') (X, d) 

are both uniformly continuous, then the topological entropies htopiT, d) and 
htopiT, d') coincide, see [Wa]. Therefore, for compact metrisable spaces X the 
topological entropy only depends on the topology generated by the metric. 

We are to estimate the topological entropy for the time-one flow T := on 
the energy shell X := Sg, and for any metric on Y^e uniformly equivalent outside 
the compact region projecting to the interaction zone to the metric induced by 
the Euclidean metric of T*(M'^— XZ). T is uniformly continuous since the motion 
is asymptotically free. 

In our case we have 

htopi^-/) = htopi^^E^): 

although in general htopi^E^,K) 7^ htopi^E^,K), since we have a symmetry 
w.r.t. time reversal. 

Prop. 114.31 below will show the importance of the set Be of energy E bound 
states. 

Definition 14.2 The non-wandering set of a continuous flow : X — 

X is given by 



:= <^ a; e X 



V neighb. f/9x,VT>0:f/n(|J I ^ 



\t>T 



For non-compact spaces X the topological entropy is in general larger than the 
topological entropy of the restriction to the non-wandering set. Nevertheless, in 
our case we have 
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Proposition 14.3 For E > Ei-^, the non-wandering set equals 

Q{<I>e) = Be, 

and the topological entropy of the flow on the energy shell is determined by the 
bound states, that is, htop{^E) = for n = 1 centre and 

V($l,) = (14.1) 

forn > 2. 

Proof. The proof parallels the one of Lemma 7.4 of [KK] . In particular 

htopi^E) = forn = 1, since then by Thm. lTTBl &F = 0. □ 

Theorem 14.4 1) For E > E^ii and n = 1 or n = 2 centres the topological 
entropy equals htop{^E) = 0- 
2) For n > 3 centres 

/.top($l,) = ■ ■ (l + + 0{1/E)^ . (14.2) 

Here is the unique solution (with eigenvector v) of the largest eigenvalue 
problem A(s) = 1 for the one parameter family of nxn Perron-Frobenius matrices 

M{s), 



v,Z)(l,v 

Ch...--= ^'th Z:=(Zi,...,Z„),]l:=(l,...,l) 

2{v,Uv) 

and the matrix D of distances d^'K 

3) Let Ne{T) be the number of closed ^^^-orbits of period smaller than T. 
Then there exists a (E-independent) constant C > 1 so that forT large 

p..ex£^ < N,m < (B > B.,) (144) 

l^top-L "■top-' 

with htop = htopi^E). 
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Remark 14.5 For configurations which form an equilateral triangle or a tetra- 
hedron, the trivial estimate 



'hop ^ 



ln(ra — 1) ln(n — 1) 



for the PF eigenvalue problem of the matrices Ai{s) is sharp. 



Proof. 1) For = 1, by Prop. I14.3l the topological entropy is zero for E large. 
By Prop. 114.31 for n > 2 we need only estimate htopi^Elbs)- 

For n = 2, we know from Thm. 112.81 that bE consists of only one closed 
orbit. Therefore, by Thm. 7.14 of [Wa] , the topological entropy htop{^E) = 

^top($U6^) = 0, too. 

2) The interesting case left is n > 3. Since /itop is a conjugacy invariant (Thm. 
7.2 of (Wal), we have by Thm. WZEl 



E\bp 



h 



top 



(14.5) 



for the time-one T^-suspension flow on (see ( 112.101) ). 

Generally, given a roof function r : X ^ R+ (see Def. 112.3 1). any cr-invariant 
probability measure yU on X induces a probability measure fir on X^ which is 
invariant under the r-suspension flow a*, namely the one obtained by normal- 
ization of the measure /i x da; on X x M, restricted to the fundamental domain 
{{k,t) I < t < r{k)} = X^. Conversely, any cr,*-invariant probability measure 
can be obtained in that way. 

Furthermore, the topological entropy is the supremum of the fCS'-entropies 
over all ergodic measures, and the same is generally true for the topological 
pressure 



/i T — invariant prob. measure 
/i ergodic w.r.t. T 



P{TJ) := sup< 


K{T) + 


fdfi 






1 X 


= sup < 


KiT) + 


fdfi 

1 X 



of a continuous map T : X ^ X of a compact metrisable space X and / G 
C(X,M), see Cor. 9.10.1 of [Wa| . 

Finally, any cr*-ergodic measure fi,. on X^, comes from a ci-ergodic measure 
jjL on X. 
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So we can apply Abramov's formula 



(14.6) 



(which is valid for any ergodic measure /i on X, see Theorem 2.1 of Chap. 3 in 
Sinai [Sin]), to obtain the formula 



htop{(rl) = sup !^h^{a) ^ J r dfi 



fi ergodic 



Eq. ( 114. 7p implies the scaling behaviour 



top 



(AG 



and the inequality 



htop{crlj > htop{crlj (ri < 



(14.7) 



(14.8) 



(14.9) 



Using Lemma [T0.6l for estimating the Poincare return time Te, we obtain for 
C large 



(1 - C/E) ■ r{E) < V2E ■ Te < (l + C / E) ■ r{E) {E > Eth), 



with 



(In F'\ 

So by applying ( llAgi) and (fTO) . we get 

/itop(4) = htopial^E)) ■ ■ (1 + (1/^)) • (14.10) 
The claim (114. 2p follows from (114. lOp and (114. Sp . if we can prove 

htopi^liE)) = hTop ■ (l + + 0{1/E)^ . (14.11) 

This is done by reformulating it as a question about the topological pressure 
P(cr, -s ■ r{E)) for the shift a. We know that for s = 

P((T, 0) = /itop(a) =ln(r2-l) >0, 
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since \n(n — 1) is the largest eigenvalue of the transition matrix of the shift, 
and since n > 3. On the other hand, for E > E^h the values r(E)(k) are 
larger than (imin/2, so that the roof function r = r{E) E C(X, M+) has mean 
Jj^rdix> (imin/2 w.r.t. to any probability measure fi. Thus P{a, —s ■ r) < for 
s > 21n(n — l)/(imin, and there exists an s > with 

P{a,-s-r) = 0. (14.12) 

By convexity of the map / i-^ Pier, f) (Thm. 9.7(v) of [Wa]) this s is unique. 

Moreover the roof function r is locally constant, so that there exists a unique 
equilibrium state /i for (114.121) (see [Ru]), i.e. 

h^{(y) = s I r dfi and hv{a) < s I r dv for 
Jx Jx 

Thus htopicrl) = s, using (114.70 (and ergodicity of yU, which follows from unique- 
ness, see Thm. 9.13 of [Wa] ). Since cr : X ^ X is an expansive homeomorphism, 
by Thm. 9.6 of [Wa| we can use the formula 

P{(yJ)= lim — In (p„(cT, /,«)), 
with the generating partition 

a = {Ai, . . . , An} with atoms A; := {A; G X | /cq = /} 

and 



Pmio; /, a) := inf < ^ sup exp{Smf{k)) 



(3 finite subcover of \J a 



m—1 
i=0 



Here Smf := T.T=o^ / ° Now for 5 := \n{E)/E 
exp(-s ■ Smr{E){k)) = 

exp (-^ (rf'^-^" - + (Z,„ - Z,Ji5)) ■ n M{s,6),,,,^^,, (14.13) 

with the n x n-matrix Ai{s,S) given by 
M{s,6)ij : = 



m—1 



i=0 



exp(-s(rf*'^-|(Z, + Z,))) 

,^ = j 
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So the function is constant on the atoms of the partition P := \/^^ T~*a. 
The first factor in (114.131) is bounded from below and from above, uniformly in m. 
The PF property of the symmetric matrix A^(s) then implies that P(cr, —s-r{E)) 
is equal to the logarithm of the largest eigenvalue \{s,5) of ^A{s,5). 

We now do perturbation theory around A4{s, 0) = A4{s), with A4{s) defined 
in (114.30 . and obtain from the condition A(s((5),5) = 1 

^ D2Xisi6),6) 
^ ' D,X{s{6),6y 

Since the PF eigenvalue is isolated, the derivatives exist, and since Ai{s,5) is 
symmetric, they equal at 5 = 

DiA(s(0), 0) = {v, DiMis{0), 0)v) = - {v, Dv) 

and 

D2A(s(0), 0) = {v, D,M{s{0), 0) v) = (tT, (l, . 

So 

s{5) = siO) + s'i0)6 + OiS') = /i- ■ (1 + C,,J + 0(5')), 
proving (114.111) . 

3) Finally, for fixed energy E > Eth the estimate (114.41) for the number of 
periodic orbits follows from Thm. 2 of the article [PP] by Parry and Pollicott (or 
alternatively by arguments based on the Renewal Theorem): If the Axiom A flow 
0"^ is topologically weak-mixing, then 

exp{htopT) 



NEiT) 



Otherwise all periods are known to be integral multiples of some Tq >, and the 
formula 

1 - exp(-/itopTo) /itopJ 
follows directly from Thm. 2 of |PP| . Here To is the largest positive eigenperiod, 
which is bounded by 



To < inf Ts(fc) = 0{l/^E). 



fcex 



As we just proved that /itop = 0{y E), htop ■ To has an i?-independent upper 
bound, and we may thus use an E'-independent C > 1 in (114.40 to bound the 
r.h.s. of (114. 14p from above and from below. □ 
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Remarks 14.6 1) Observe that, up to a small error term, htopi^^) indepen- 
dent of the charges Zi of the centres. Effectively an attracting potential speeds 
up the particle a bit, which then leads to the Z-dependent correction term in 

2) It is natural to ask whether one may improve (114. 4p to show 



iv.(r) ~ ';B!i!hfl, (14^15) 

However, the return time estimate of Lemma [TO. 6l for Te{x) is, up to the relative 
order 0{1/E), independent of the point x E Ve(A;_i, /cq, /^i). 

So we cannot exclude by that estimate that all return times are equal to a 
constant Tr for a symmetric configuration (e.g., for an equilateral triangle or 
tetrahedron with equal charges Zi). But by formula (114. 14p in that case Ne{T) 



is not asymptotic to any smooth function, so that (114.150 does not hold. In 



this sense our statement is optimal, given the, already quite precise, return time 
estimate. 

Any eventual improvement of Lemma [10.61 must be dependent of the addi- 
tional smooth component W of the potential V, and thus be complicated. 

Although the iterated Poincare map is certainly mixing for the measure of 
maximal entropy, by the above argument we cannot decide whether or not the 
flow is mixing. 



15 Characterization of the Scattering Orbits 

Up to now we were mainly concerned with the bound states Be C T.e- However, 
the topological entropy analyzed in the last section is an example for a quantity 
which, though a priori depending on the dynamics on the whole energy shell T.e, 
is determined by that subset of Liouville measure zero. 

When we now consider the scattering states se C T,e, again their description 
will be based on symbolic dynamics of the botvnc/ states. 

Our concrete question will be to enumerate the scattering orbits with given 
energy E > and asymptotic directions 9"^ E S^, that is, the subset 

iH,p-,p+r\E,e-,e+) 



of phase space (noticing that the asymptotic defined in (16.140 are constant 
on orbits). 

Before we come to the case of several centres, we first consider the simple 
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Example 15.1 : Keplerian motion. 

Without loss of generality V{q) = — j^. There are no undeflected orbits, i.e. 
only asymptotic directions 

{§-,6+) e{S^ X S^)\A, with diagonal A := {{6,6) \ 6 e S^} 

occur. The pericentric time T is a smooth function on the positive energy part 

P+:={x eP \ H{x) > 0} 

of phase space, see (15. 7p and the definition of the manifold P. In fact in that 
case the map 

{H,p-,p+,T) : P+ ^ M+ X {{S^ xS')\A)xR 

is a diffeomorphism. In particular for > and asymptotic directions 6^ 7^ 6^ 
there is exactly one Kepler hyperbola. 

If we allow for an additional smooth potential W, we need to exclude a whole 
neighbourhood of the diagonal A C S'^xS'^, since even for large energies forward 
scattering may be dominated by W and not by the Coulomb potential. 

Already in the case of bounded orbits we introduced the NC condition of 
Def. 12.41 which excluded collinear configurations of nuclei. Similarly we need to 
exclude certain asymptotic directions if we want to obtain statements which are 
independent of the precise form of the potential V. 

Thus for i!) E (0, vr] and 



i,k\ 



Qn^in : [0, vr] e^in(^) := min <{6, s' 

l<i^k<n 

(letting 6min(^) := vr for n = 1 nucleus) we restrict our interest to the asymptotic 
directions in 

AD{^) := {{6-, 6+) eS'xS^ I mill (<(r , ^+), e^in(r ), 0^i„(^+)) > i?} , 
and the ($*-invariant) subset 

se{^) ■■= {xEseI {p-{x),p+{x)) e M)ii3)} 

of scattering states, thus excluding near-forward scattering and scattering from 
or to any direction near an axis through two nuclei. 



Ill 



For n > 2 centres we have bound states which influence the scattering 
trajectories in their vicinity. As the high energy bound states can be described 
by symbolic dynamics, we introduce the set 



W := U 



fcGN 



of (non-empty) admissible words to enumerate the scattering states, see (112. ID . 

Every orbit within S£;(^9) will be uniquely characterized by its asymptotic 
directions and by the sequence of its near-collisions with the nuclei (^^-visits in 
the sense of Def. 112.60 . Conversely every admissible word will be shown to occur 



for some orbit. 

Namely we set the angle parameter controlling the asymptotic directions AD 
equal to 

^ = ^iE) := miniCj/VE, a^j2), (15.1) 

with Cj to be fixed in Thm. 115.31 

Similarly, the angle parameter controlling the near-collisions is fixed by 

e = 9{E) ■= min(C9//E, amin/2) with Cg := 4ci. (15.2) 

In Lemma WA\ the parameter ci controlled what the regime of hard scattering. 
It will be fixed in Thm. 115.31 too. 

For X e se let w{x) := if A/'C7^(s) = 0. Otherwise the orbit through x 
enters the interaction zone, i.e. intersects the hypersurface dVj^ at a unique 
point x' , and we set 

w{x) := k{x') G W, 

where the trajectory 1 1-^ 6'-visits the centres in succession k (we start at 

x' instead of x since we want the 9 — visits to occur for positive times so that 
kE'Kl for some r). 

Lemma 15.2 For Cg > large and E > _Eth every orbit in se{Cs/E) intersects 
the hypersurfaces dV^ exactly once, so that the entrance, resp. exit times 
on SEiCs/E) are uniquely defined by ^{T^{x), x) E dV%. 
The restrictions ofT^ and Tj^ to se{C^/ E) are smooth functions. 

Proof. If a scattering orbit in se does not meet the interaction zone, then 
by the virial estimate (12.151) there is a unique point x = (p,q) on this orbit 
with {p,q} = 0. By (12.131) the speed is bounded below by |p| > \/E. Thus 
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the direction p = p/\p\ differs from the asymptotic directions (defined in 
f 04l) ) only by 

using (16.61) . Choosing a large enough constant Cg, we see that 

<{p,p-{x)) + <(p,p+(x)) < Cs/E \fE> Eth. (15.3) 

Thus the orbits in se{Cs/ E) intersect dV^. 

For these scattering states we know from (12.150 that there is a unique en- 



trance resp. exit time. As long as the orbits intersect the C°°-hypersurfaces dV"^ 
transversally, smoothness of follows from smoothness of the flow $*. But 
this transversality can be enforced by further enlarging Cg. □ 

Theorem 15.3 We assume that V satisfies ttie decay estimates Ii6.16\) . Then 
for Cj > large in ( IJ5. 1\} . and E > Eth the map 

DiSe ■■ se{^{E)) AD{^{E)) xWxR 

X {p'{x),p^{x),w{x),Tj^{x)) 

is a difFeomorphism. 

Remarks 15.4 1) In particular the orbits in S£;(i9(-E')) are uniquely character- 
ized by their asymptotic directions and the sequence of near-collisions, and any 
sequence is realized by an orbit. 

2) In this section and Sect. [Hlwe exclude from our consideration small cones of 
aperture 0{l/\^) around the axes through two nuclei. 

This is indeed necessary, since there non-universal phenomena occur. Also, 
the order 0{1/\/E) of aperture is optimal, as we show now by example. 

We consider the purely Coulombic two-centre problem (see Appendix B) with 

:= Qj, S2 := (^V), Zi > and Z2 = -Z^. 

We claim that for energy E > _Eth there is no orbit colliding with si and hav- 
ing scattering angles 6{oo) w.r.t. to the negative 1-axis smaller than ^^/\Z^\JE, 
whereas there are such orbits for 6{oo) > 2^\Z^^\JE. 

Such an orbit would lie in a plane containing s\ and S2, say the 1 — 2-plane. 
By reflection symmetry we need only consider the part of the orbit after collision, 
and we denote by r the time of its last intersection with the gi = 0-plane. Then 
^(^(^)) = so that piij) = -a/2(E -p|(r)). So at time r the angle of the 
particle direction p(r) with the negative 1-axis equals 

e{T) = arcsin(|p2(r)|/V2E). 
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As we are interested in small scattering angles, we may assume |pi(t)| > |p2(t)|- 
Moreover 

as follows from the estimates of Sect. O Thus the angular momentum w.r.t. S2 
equals 

IL2I = Mr) - S2) X p{t)\ = \p2{r) - g2(r)pi(r)| = 2\p2ir)\ + 0{l/E). 

As the scattering angle A6' in a Coulombic potential of charge Z-2 meets the 
relation 



sinA0 |Z; 



2 



I+COSA6' y/2E\L2\ 

we obtain for l/\fE <ti 1^2(^)1 ^ \fE the total scattering angle 



e{oo) = e{T) + Ae-{i + o(i)) = {2Ey 



-1/2 



\P2 



J2\ 



\V2\ 



(1 + 0(1)), 



which is minimized by \p2\ ~ a/|^2|, with value 



^^(oo) = v/2|Z2|/E- (1 + 0(1)). 

Proof of Thm. [TOl 

As remarked before Lemma 18.41 there are two regimes for Coulomb scattering, 
hard and soft scattering. A large choice of the constant C7 in Thm. 115.?] separates 
the two regimes. 

1) In order to prove the theorem, we thus first have to fix the constants Ci,C2 
which appear in Part 1) and 2) of Lemma 1^341 for ^ = \- The lemma describes 
the linearization of the Kepler Transformation, but its estimates are also valid for 
the motion in the potential V , using estimate (I8.20p and (18.211) . 



We are interested in the effect of that linearized scattering transformation on 
cone fields of the form 

C{a):={{%)\\bu-bv\<a\bu^bvW (a > 0). (15.4) 

So we have C(6) C C(a) for 6 < a (note for comparison that in (lll.ip we used 
the cone field with energy-dependent a = C/E). 
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• e > is chosen small enough so that the (according to Lemma ITTT] nearly 
free) linearized flow within XZ(c^) maps the coneC(l + 2£:) into C(2), and 
the cone C(|) into C(l - 2e). 

For such a small £ > we then choose 

• ci > large enough so that the linearized scattering transformation maps 
the cone C(2) into C(i) 

• C2 > small enough so that the linearized scattering transformation maps 
the cone C(l + e) into C(l + 2e). 

All these choices can indeed be made (apply the relevant unperturbed matrices 



^g) and (/jO) from 1^ and (O), resp. (KWf on the vector 



By fixing the constants ci and C2, we have defined soft and hard scattering. 

2) Now all 6'-visits in the sense of Def. 112.61 (meaning that the orbit of energy 
E locally intersects the Poincare surface 'Hi{9{E))), are hard scattering events. 

Namely, as 9{E) = Aci/y/E (see (115. 2p ). the scattering angle Atp inside the 
ball Bi(cq) is bounded below by 

Atfj > e{E)/4 = ci/Ve, 

using (19.31) . It thus meets the criterion of Lemma l8T4] for hard scattering. 

3) We show now that, for large enough Cj in (115. ID {i.e. by excluding large cones 
of asymptotic directions p^(x) around the axes through two centres), Diff^; is 
well-defined, that is, that every scattering state x G SEi'&iE)) 6'-visits at least 
one nucleus, so that NCTq^x) ^ 0. 

In view of Lemma [T5.2l the orbit through x E SE{'d{E)) meets the interaction 
zone. If the directions at the times T^{x) of entrance resp. exit are denoted by 
pi resp. po, then like in (115. 3p . we have the estimate 



<{Pi,P~{x)) + <(po,p+(x)) = 0{llE) {E > Eth). 

Since <{p^ (x) , p^ (x)) > "d^E), we conclude that for large enough Cj the change 
of direction inside the interaction zone is at least 

<{pi,po) > ^^E). (15.5) 

Lemma 17.11 then tells us that the trajectory enters at least one ball of (energy- 
independent) radius Cq around a singularity s/. 
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We count the number of such occurrences. If > 3, then Prop. 19.2] says 
that the orbit «inin/2-visits at least one nucleus. Comparing with the definition 
([152]) of 9, this is a 6'-visit. 



But the same holds for N < 2, since otherwise 

• the changes of directions inside the N Cg-balls are each smaller than 9 (Lemma 
lOl . and 

• the components of the trajectory inside the interaction zone but outside 
the Cg-balls each contribute at most with C/E to the change of direction (Lemma 
O). 

• Adding these contributions and using the definition (115. ID of 9, we would get 



<{p^,Po) < ^ + 9{E) = ^ + ^ < = 1^{E) for large, 



in contradiction with (115. 5p if 

Cj > 4C9. 

4) We want to avoid intermediate scattering events with angles meeting 

^2 < < (15.6) 



E VE 

This can be done, too, by choosing Cj in (115. ip large enough. 

From 3) we know that the orbit 6'-visits the centres in some succession k = 
{ki,...,kr). There can be at most one visit of some ball Bk^Xcq) before and 
some ball Bk^^-^(cg) after this sequence of near-collisions (Prop. 

Moreover, for large Cj, these visits, if they occur, are soft scattering events 
in the sense of Lemma 18.41 We need only prove this for Bk^^{cq), as the result 
for Bk^^-^{cq) then follows by time reversal. 

So we consider the half-orbit $((—00, ti], x) through x with ti = mmMCTg{x) 
see Def. 112.61 This half-orbit ends in the Poincare surface near s^^, namely 
G 'Hkii9). More precisely, by def. (18. ID the configuration space distance 
is bounded by 

\q{ti,x) - SfcJ < ^ . 

CgV E 

The half-trajectory (f((— 00, ti], x) consists of three types of segments: 

1. The segment outside the interaction zone XZ Here the total change of 
direction is of order 0{1/E) (Thm.lO). 
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2. The (one or two) segments in XZ{cq). Here, too the total change of 
direction is of order 0{1/E) (Lemma [7.11) . 



3. The segment inside the ball i?fc„(cg), perhaps empty. We know that the 

half-orbit has no hard collision here, i.e. it misses Hko^O) (otherwise k 

would begin with ko instead of ki)). Thus the total change of direction 
inside Bk^(cg) is bounded by 



<e = 




Summing these contributions, for i?th \3rge the total change of direction on the 
half-orbit is bounded by 

2C 

<(p-{x),m)<-7= (-oo<t<r), (15.7) 



T being the time where the half orbit enters the ball Bk^{cq) 
On the other hand by our assumption x E SEi'&iE)) 



<{p-{x),s''^''^)>-^. (15.8) 



(v ,w ) = (j){T,x)/V2E,{q{T,x) — SfcJ/cqj are the coordinates (18.30 of the 
end point $(r, x) of the half-orbit. By Lemma 

<[V , -W ) < 

so that with (115.70 the relative position meets 

<iw-,-p~{x)) <<{p-{x),v-) + <{^r,-w-) < (^2C9 + ^^^^^^ ^ Ve. 

Together with (115.71) and (115.81) this means that for Cj large the distance of 
q{t,x) from the line L through and is bounded by 

C 

dist(g(t, x),L) > \—/= ■ |^(i, — Sk^ \ (— oo < t < r). 
V E 
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But then 



dist(g(t,x),4J > 



1 c, 



'-00 <t<T) 



SO that there is at most so /t scattering inside Bkg{cg), if Cy is large. Thus there 
are no intermediate scattering events with angles meeting (115.60 . 

5) The last remark implies that the symbolic sequence x i-^ w{x) is locally 
constant on So by Lemma 115.21 and Thm. 16.51 the map Diff^; is 
smooth. 

6) Our next task is to show that Diff^; is onto. It suffices to find for given data 

{e-,e+,k) eADi^) xW (15.9) 

an orbit in T.e with these asymptotic directions and ^^-visits. 

We proceed in a way similar to the construction of bounded orbits in Sections 
[11] and [121 If the symbol sequence equals k = {ki, . . . ,kr), then we erect 
incoming resp. outgoing Poincare surfaces V{6~,ki) resp. W{kr,6'^) with 



Vie, I) := {x^ip,q)EV 



-1^ 



< 



c 



10 



P 



2Cio 



< 



W{1,9) 



X = {p,q) 



W- si\ = c. 




9 G 5^, near qi := si — Cg9 resp. near qo := si + Cq9 (compare this definition 
with the one of 7^^'^ in ( l9lB ). 



For all q^ in the configuration space projection of V{9 ,ki) and all q 



in 



the configuration space projection of W{kr,9^) there exist p^ with {p^,q^) G 
V{9',ki), {iP~,q^) G W{kr,9+) such that the orbit through {p~,q^) 6'-visits 
the centres in succession k and then meets {fP~, q^). The proof of this assertion 
uses the estimates of Sect. [HI slightly modifying Prop. 111.51 
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We set Cio := C7/8 and choose a large value of Cj. As 6min(^'^) > '&{E) = 
Ct/^/E, by what we have proven in part 4), V{6^,ki) did not ^^-visit a nucleus 
in the past, and the orbits through W{kr,9^) will not ^^-visit a nucleus in the 
future. Instead by Prop. they leave the interaction zone in time ^0{\^). 

Using estimate ( 115.71 with t := r, we see that for every q~ there is at least 
one such that the orbit through {p^,(^) E V{9^,ki) has asyptotic direction 

lim p{t,ijr ,q~)) = 6~ . 

By standard arguments this family of orbits contains one 0-visiting the centres 
in succession k and then having limit 

lim p{t, {p~ , (T)) = 0^ ■ 

So Diff^' is onto. 

7) Cone field estimates based on Part 1) now show that Diff^; is one to one and 
smoothly invertible. Specifically we show that there exist invariant cone fields 
C{a) (see (115.41) ) along the orbits constructed in Part 6). 

The idea is simply that there is at least one hard scattering along the scat- 
tering orbit, making the family of energy E configuration space trajectories with 
initial asymptotic direction 6^ divergent for large positive times. 

So for given 6^ we consider the Lagrange manifold of points x E se having 
asymptotic direction p^{x) = 9~ . We claim that for Eth large the tangent space 
of this submanifold at (p^ , q' ) E V{9^ , ki) is contained in the (wide) cone C(2). 
Namely the half-orbit ending in consist of at most four segments: 

1. The segment outside the interaction zone IZ, controlled by Thm. 16.31 For 
large E^h this is in a cone C(l + e/2). 

2. The segment in XZ(cq) before entering the ball Bk^(cq), controlled by 
Lemma mi Here the linearized Lagrange manifold is contained in C{l + e). 

3. At most one segment in Bk^Xcq), controlled by the choice of the soft 
scattering constant C2 in Part 1) of this proof. So after this soft collision, 

we are in C(l + 2e). 

4. The segment in XZ{cq) after leaving Bkf^{cg) and before entering Bk^^Cq). 
By Part 1) here the linearized Lagrange manifold is contained in C(2). 
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Our choice of the hard scattering constant ci implies that after leaving Bk^{cq) 
we are in the narrow cone C(|). At least the same is true after leaving Bk^{cq), 
using the cone field estimates of Sect. [TTl Finally, we consider the positive half 
orbit starting at {p^.q^) E W{kr,9^) and get a similar sequence 

C(i) ^ C(l - 2e) C(l -e)^ C(l - ^e) 

of cone fields, using Part 1) again. For large energies this shows uniqueness of 
the orbit with data (fT5^ . □ 



16 The Differential Cross Section 



The scattering transformation S defined in (I6.15P contains complete information 



on the scattering process. As we have seen in the previous sections, it exhibits 
many aspects of irregularity if > 3. Nevertheless, the scattering transformation 
is not directly accessible in a (classical) scattering experiment. 

Firstly, one typically cannot fix the initial angular momentum of the test 
particle. Secondly (unlike in a quantum mechanical setting where interference 
effects exist) it is hard to measure time delay. 

What is accessible is the differential cross section -^(E,6^ ,6^). Informally 
speaking, this is the (density of the) number of particles per second scattered in 
the final direction 9^ , assuming a uniform flux of one particle per second and 
unit area of incoming particles of energy E and initial direction 9~ . 

One could expect to see some trace of irregularity in the differential cross 
section, and in fact for all systems considered up to now numerical calculations 
of the cross section indicated the existence of so-called rainbow singularities on 
a Cantor set of angles, see [Eel EJl IGal ISml [Te] . 

For the simplest case of a Kepler potential we obtain the so-called Rutherford 
cross section, see (116.61) below. It is remarkable that the differential cross sections 
for the cases n = 2 and n > 3 turn out to be very similar to the Rutherford cross 
section (see (116.80 for the statement, and Figure 12.1 of [KKJ for a numerical 



plot for c? = 2 dimensions). So the complicated structure of the time delay and 
the scattering orbits is not reflected in the cross section. 

The reason for that discrepancy is, roughly speaking, the following. For 
d = 2 the deflection functions 9^(E,9^ , L^) (depicted in [KKJ , Figs. 

10.2-10.3) are strictly monotonic w.r.t. the initial angular momentum L~ . It is 
clear from the definition of that extrema of the deflection function lead to 
singularities in the differential cross section. Since there are no extrema (except 
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for the degenerate situation L~ — > ±00), we have a nonsingular (except for 
the forward direction). 

Before stating our theorems, we shall recall a mathematically correct defi- 
nition of cross section. In the physics literature the cross section is sometimes 
introduced as a function, whereas it really is a measure. The difference is of 
some importance because in general that cross section measure is not absolutely 
continuous w.r.t. Haar measure. In our context, we shall show that under cer- 
tain conditions the cross section measure /'s absolutely continuous if one excludes 
the forward direction and the directions near s^'\ and that the Radon-Nikodym 
derivative, i.e. the differential cross section, is smooth. 

We denote by 

a : P+ — s> M"^ , (p, q) ^ q — {q, p) p 

the impact parameter. For phase space points x E P+ projecting to a singularity 
si or with X ^ (p, q) = (0, q) the impact parameter is defined by, say a(x) := 0. 
Now on the ±-scattering states C P+ the asymptotic impact parameters 

a^-.s^^M.^ , a^:= lim ao$* 

are well-defined continuous functions, since d(p,q) is continuous outside the 
interaction zone (that is, for |g| > R^i^iH {p, q})) , and there 

"'P-« = 2(77(M^17(M)- ^'''^ 

By Thm. 16.31 the r.h.s. of (116.11) has a limit on scattering orbits. 
For energy > we consider the maps 

A^:se^T*S^ , X ^ {a^{x),p^{x)) 

from the set of ± scattering states of energy E to their asymptotic data. These 
can really be considered as points in the cotangent bundle 

of the two-sphere, since 

• {a^ (x) , p"^ (x)) = , and 
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• the bilinear map (p, g) ^— {p^q) is the natural pairing between vectors and 
co-vectors (no Riemannian metric involved). 

The cotangent bundle N carries the canonical symplectic two-form uj^ and the 
volume form \n '■= \^^n A lun- 

The maps are constant on orbits, and for y E N the preimages {A'^)~^{y) 
consist of at most one orbit. 

Proposition 16.1 For energy E > the set A'^{se) C of asymptotic data 
of the scattering states is open, and its complement N \ A^{se) is a compact 
set of Xn -measure zero. 

If V satisfies the decay estimates ( 16. J 6)) . then the map 

Me: A-j,{se) ^ AI{se) , A'^^x) ^ A^^x) {x E se) (16.2) 
from initial to final data for energy E is a smooth canonical transformation. 

Proof. For energy > and radius r > R^i^(E) we consider the smooth 
Poincare surfaces 

Ul,r ■■= m e I \q\ = r,± {p, q) > 0}. 

These four-manifolds are transversal to the flow as {H,q'^} = ^{p'^,q^} = 
—2 {p, q) ^ 0. According to Lemma 8.2 of McDuff and Salamon [PS] , they are 
symplectic submanifolds of our phase space P, and the Poincare section map 

UE,r^SE^Ul,^nSE (16.3) 

which send x to the unique intersection point of the orbit x) with ^ is 
a symplectomorphism. 

We use on ^ the coordinates (a,p), which are maps 

Indeed they are diffeomorphisms onto their common image 

Nr := {(a',p) e N \ \a'\ < r}. 
Thus the map ( 116. 3p induces a diffeomorphism 
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which converges pointwisely to (116.21) asr — oo, using Thm. 16.31 With assump- 
tion (I6.16p . smoothness of is then implied by Thm. 1^751 

The scattering states s C P form an open subset (see Thm 2.3.3 of [DG]), 
and thus U^^ fl se is relatively open, too. All orbits in T.e which do not meet 
the interaction zone XZ{E) are scattering. Thus \ A^{se) is a compact set. 
Then using A%{s%) = N and asymptotic completeness Corollary |6.4||^ we see 
that Xn{N\A%{se))=Q. □ 

For As2-almost all 6~ G S"^ the map 

: TIS^ ^ S\ (pE,e-in :=p+(E,a-,r) (16.4) 

is measurable. Here the restriction of the asymptotic direction p+ : s+ ^ 5*^ 
(which is constant on orbits) to se is considered as a map 

see Prop. 116.11 

Definition 16.2 For E > and ^ G S*^ the cross section measure (t{E, 9^) 

on S"^ is the image measure 

\g- being Lebesgue measure on the cotangent plane T? S*^. 

Assuming a{E,6^) on 5^ \ {9~} to be absolutely continuous w.r.t. Haar 
measure A52, the differential cross section -^{E,6^ ,6^) is the Radon-Niko- 

dym derivative of (t{E, 6^). 



In Def. (116. Sp we by using Ag_ we normalize the flux through unit area in con- 



figuration space to equal one. 

Example 16.3 By radial symmetry, for the simplest case (16. ip of scattering by 
a Kepler potential with Z = 7^ 0, the (Rutherford) differential cross section 
depends only on the angle 

A9 := <{9-,9+) 

between the initial and final direction. Using formula (I8.12p for the eccentricity 
e = +1/ sin(|A6') of the Kepler hyperbola, we see that the modulus a of the 
impact parameters ar^ equals 

« = ^cot(iA^^), 
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so that 



da 
dM 



\Z\ 

AE sin2(lA^) 



We may assume that 9~ = (0, 0, 1) so that the 3-component of d~ = p~ x L~ 
vanishes and Ag_ corresponds to integration with the two-form dL^ A . 
Introducing polar coordinates (L,ip^) in the plane T^_S'^, and expressing the 

volume element at 0^ on S"^ in the form 

sin(A6') ■ d^p- A dM, 
we obtain the familiar expression 



da 
d§+ 



{E,9-,9-^ 



Ru 



L dL 



2E dA9 



AE sin' a A9) 



(16.6) 



for the Rutherford cross section. Note that it depends only on the modulus of 
the charge Z. 

For {9-, 9+) e Q;^M,keW let 

a, {e,9-,9 
with Diff^; defined in Thm. [1531 



lim a(Diff^^(t,r,^+,^)), 



Theorem 16.4 Let V be a Coulombic potential satisfying the decay estimates 

Then for energy E > E^h and i!) = rmn{cE~^ , a^in) with < 5 < ^, on 
Qmin(^) differential cross section is smooth, of the form 



da 
d9+ 



iE,9-,9+) = J2 detf 



da 



{d-,{E,t,9^)) 



and differs from the Rutherford cross section l\16. 6\) for charge Z 
only by 



da 



{E,9-,9+) 



da 



d9+ \d9+ 
uniformly in {9-,9+) G 



{E,9-,9' 



(16.7) 
(16.8) 



Ru 
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Remarks 16.5 1) In particular the relative difference w.r.t. Rutherford cross 
section is only of order 0{1/E) if one excludes cones of an energy-independent 
aperture 

2) In 2-dim. and purely Coulombic potentials the differential cross section is even 
smooth (up to the forward direction 9^ = 6^) for all positive energies [KK] . As 
shown in |Kn3] . this smoothness is rather exceptional. 

Also here for n > 1 centres we cannot have smoothness of the differential 
cross section (for no V and no E > Eth), if we add to AD^i!)) the neighbourhood 
of any direction s''^. This follows from the observation that (contrary to the 2D 
attracting case) in 3 dimensions any hard collision with a nucleus changes the 
degree. By a limit argument there must then exist points where the degree is 
zero. At these points the differential cross section diverges. 

Proof. Thm. [TS^l says that the orbits with data {E,6^,6^) are enumerated by 
W. So if the r.h.s. of (116.70 converges, then by the Transformation Theorem for 



Lebesgue measure (116.70 follows from the definition of the cross section measure 
in I^EB- 

Estimates (1617I) and (108l) of Thm. [63] imply that for |go| > -Rvir(-E), 
± (go,Po) > and multi-indices 7 := {a, (3) eN^xNl 



and 



d2,{p^{xo)-po) = o(^E-'-h-\^ 
dl{a^{xo)-a{xo)) = o(^E-'-h-\)^ 



the last estimate being obtained with the help of (116.10 . 

Similar statements are true for orbit segments in XZ{cg). 

So up to an error of order 0{1/E), all variations of the asymptotic data 
come from the single scattering processes within the balls of radius Cg around 
the singularities. 

We switch to (y , £)-coordinates. 

There are two types of such contributions: 

1. The ones coming from the hard collisions (6'-visits in succession k). These 
lead to factors of the form ( 111. 21) in the product formula for the lineariza- 
tion: 

T,VE = f{k.uko,h)E-(l 1]+0{E^) (16.9) 
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with 



f{k_i,kQ,ki) :-- 



ko 



2(5- 1^ 



As A'ip > cE ^ , the relative error in this estimate is of the order 0{E 

Visits of Cg-balls around some singularities, which are not hard collisions 
(6'-visits). By what we have shown, there can be at most two such events, 
one before the fc-visits, one after. 

So these visits meet the hypothesis (18. 6 p of Lemma 18.21 

'^— + w' >cE~^), 

i.e. with angle i3 = 0{E-^'), where 5' := 1 - 5 e [i, 1]. 

So formula (I8.10p of Lemma 18.41 says that the relative deviation of the 
linearized flow from free motion during these soft collisions is of order 



5ir^ = Sir+ oic^E'-^'' 

Sw+ = 2uSv- + 5w- + 2v-du + 0{c^E^-'^^' 
Now we see that the formal sum (116.70 converges: 

• There are exactly n ■ [n — 1)'^^ words A; G W of length / 

• For word length / + 1 of ^ G W the term 



{\5v-\ + \5w-\), 
{\5v-\ + \5w-\). 



det 



da 



{d-,{E,9-,e^)) 



in (fTO) is only of relative order 0{E-^'^-^'^) = 0{1/E^), compared to 
the term of the word shortened by one letter, as there is one extra factor 
E coming from (116.90 . Here d denotes the dimension, so = 3. 

So if > Eth and the threshold E^h is suitably chosen, this decay outweighs 
the exponential proliferation of words with given word length. 

The comparison in (116.80 with Rutherford cross section (116.60 for squared 
charge = J2^=i Zf follows by adding the contributions in (116.70 of word 
length one. \n d = 3 dimensions the leading errors of order 0{E'^^^^) come from 
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1) and 2) above, whereas neglecting the contributions of the longer words is only 
of order □ 

In d dimensions the Rutherford cross section equals 



^{E,9-,9+)] =(-— ^i^-— ) . (16.10) 

In 12.4 of [KK] we remarked that for d = 2 dimensions, Zi > and ^ 9~ 
the differential cross section of the n-centre problem converges to the d = 2 
Rutherford cross section as — > oo. The charge Z in (116. lOp . however, must 



be chosen as Z = ±Yl^=i^i (^nd not Zoo, as wrongly stated in [.KKJ). This 
result can be sharpened: 

Corollary 16.6 The analog of formula ( I J 6. g|) holds true in d = 2 dimensions if 
one sets Z := J2^=i 1^*1 ■ 



Proof. Up to error terms the formulae (lll.2p and (IS.lOp for the linearization of 



the flow which we used to derive the result (116. 8p \n d = 3 are invariant under 
rotations. Thus (116. 8p is true in = 2 dimensions, too. □ 



17 The Collinear Case 

in this section we show by counterexample that the non-collinearity conditions 
cannot be dropped in Thm. \TI^ and Thm. [1531 

1) We first consider the set 6^; of bounded orbits of energy E for attracting 
Coulombic potentials V which are rotationally symmetric w.r.t. some axis A C 
Thus in particular the nuclei are situated on that axis: 

si,...,s;eA (17.1) 

Conversely, that condition ensures that V is rotationally symmetric around A, if 
it is a purely Coulombic potential. 

By symmetry trajectories with initial conditions tangential to a two-plane 
F C ]R| containing A stay in F, and we may thus consider the restricted two- 
dimensional motion on F. 

The axis A is divided by si, . . . , Sn into n + \ closed intervals meeting only 
in their endpoints si. These intervals correspond for energies E > V^ax to 
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trajectories which are reflected by the nuclei. Thus two of these trajectories are 
unbounded and n — 1 bounded. 

These special trajectories are of course invariant under rotations around A. 
However, there cannot be any further energy E trajectories in F having this 
property. 

On the other hand it has been shown in |KK] that for n > 3 nuclei above 
some energy E^h there is a Cantor set of bounded trajectories in F. So in 
particular there is an uncountable number of trajectories in F which are not 
moving tangential to A and thus give rise to one-parameterfamilies of trajectories 
for the full motion in M. 

2) In general there are bounded orbits which do not lie in any plane F containing 
the axis A. We observe that by rotational symmetry the component (^L,,s^ of 
angular momentum in the direction s of A is preserved. We now indicate that 
for certain collinear configurations there exist two-parameter families of bounded 
orbits of a given energy, parametrized by that angular momentum component 
and its conjugate angle. 

For the sake of simplicity we consider a 3-centre potential 



V{q):=~Yl 



1(7 — 

1=1 



with Si := 0, S2 := d ■ s, S3 := —d ■§,§:= (0,0,1) and Z2 = which, 
in addition of being axially symmetric w.r.t. the 3-axis A = M • s, is mirror- 
symmetric w.r.t. reflection by the plane F12 := G M| | ga = 0}. We first 
consider periodic trajectories with angular momentum component L3 = in, say, 
the plane F = F13 := {g G M| | ^2 = 0} which are invariant under reflection by 
the plane F12: 

qi{-t) = qiit) , q2{t) = , q^i-t) = -q^it) {t G M). 

By symbolic dynamics arguments combined with Thm. 6.11 and Remark 11.2.2 
of [KKJ for all energies E > there exists a countable infinity of these or- 
bits, all being hyperbolic and having index (as orbits in the two-plane F13. 
Thus by invoking an implicit function argument one may show the existence of 
smooth family of energy-i? bounded orbits starting on the 1-axis near g(0) and 
parametrized by L3. Rotating these orbits around the axis A then yields a two- 
parameter family of bounded orbits. 
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3) For repelling axially symmetric potentials the situation is completely differ- 
ent. W.l.o.g. we again consider potentials V which are invariant w.r.t rotations 
around the axis A = M ■ s, with s = (0, 0, 1). But now we assume that 

{VV{q},q-{q,s)s) <0 (geM|\A). 

This condition is met, e.g., by repelling {Zi < 0) purely Coulombic potentials 
meeting (117. ID . 



Now consider a trajectory starting at (p(0),(f(0)) with 

m - (g(0), s)s + and (Vp(0), g(0) - (g(0), s)s) > 0. (17.2) 

Then 

(p(t),g(t)-(g(t),S)S)>0 (t>0) 
and is monotonically inreasing in t, since 

4 {P^Q- = - {yV{q),q- {q,s)s) + {p,p- {p,s)s) > 0. 

at 

But this means that the orbit leaves the interaction zone in finite time and thus 
is not bounded. The second of the conditions in (117. 2p is not restrictive, since 
otherwise we may consider negative times. 

We conclude that the only bounded orbits lie on the axis. Thus for E > \4iax 
there are exactly n — 1 bounded orbits, compared to the uncountable infinity of 
bounded orbits in the case of NC configurations and n > 3. 

4) Finally we consider scattering orbits. Already for = 2 nuclei and attracting 
purely Coulombic potentials we have one-parameter families of orbits of a given 
energy E scattering from a direction parallel to the axis A through the positions 
si,S2 and to the backward direction. These are obtained by rotating a given 
solution around A. There are infinitely many such families, as can be seen from 
the explicit Jacobi solution of the two-centre problem, described in Appendix B, 
or from [KK], Thm. 12.1. 



A Aspects of Geometry and Global Analysis 

This article on the 3-dimensional n-centre problem is based upon analytical per- 
turbation estimates. To the contrary the 2-dim. n-centre problem (and similarly, 
the 2-dim. periodic potentials of [Knl[|Kn2] ) was treated in [KK] using techniques 
from Riemannian geometry and global analysis. 
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In this appendix both approaches are compared. 

It is known that the trajectories of energy E > snp^V{q) generated by a 
Hamiltonian function H : T*M ^ M of the (local) form 

d 

k,l=l 

on a (i-dimensional Riemannian manifold {M,g) coincide (up to a time repara- 
metrization) with the geodesies in the so-called Jacobi metric qe on M, confor- 
mally equivalent to g: 

gE{q):={l-V{q)/E)g{q) (1.1) 

(here we assume for simplicity E > 0). 

In the simplest case covered by the paper V{q) = —Z/\q\ with Z > 0, i.e. 
the attracting Coulomb potential. There, using a formula from Spivak ( [Sp] , p. 
337), we obtain the expression 

'''^^^ 2E {\q\ + Z/E)^ ^ ' 

for the sectional curvature of the Jacobi metric in the 1-2-tangent plane at q. 
We can learn several things from that formula: 

• Setting gs = and thus considering planar motion, 

for positive E, and this expression is bounded below, by 

K,M > (1.3) 

This fact was used extensively in [Knl] , [Kn2] and |KK] to analyze planar 
motion by going to the smooth branched covering surface 



M:= \ {q,Q) eCxC 



1=1 



of the configuration plane, equipped with the lifted Jacobi metric (the 
branched covering being given by projection to the first factor q). 
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• For d = 3 dimensions the sectional curvature (11.20 is neither uniformly 
bounded in the q variable nor definite. In fact, for q = (0,0,^3) 

^i,2(g) ~ 4^ - +00 (g^O). (1.4) 

We thus consider here geodesic motion in mixed sectional curvature. Al- 
though the i?-dependence of (11.30 is quadratic, whereas (11.40 is only linear 
in E, negative curvature does never dominate positive curvature in our es- 
timates. Namely we have seen in Lemma IQ.ll that the minimal distances of 
the bounded orbits from the nuclei are of the order 1/E, so that effectively 
(n~4l) . too, goes like E"^. 

Mixed curvature dynamics is rather intractable in general. However, in the 
case considered here the motion near the singularities can be treated as a 
perturbation of Keplerian motion, and this allows us to control the motion 
in the high energy limit. 

Whereas the estimates of this papers are somewhat optimal in the high energy 
limit E > Eth, nothing much could be said about the energy region < E < E^i^ 
(the exception being Sect. E]). 

To the contrary, for two dimensions many results were shown for all positive 
energies, using the negativity of Gaussian curvature and the topology of the 
branched covering surface M (whose fundamental group is non-abelian for n > 
3). 

The branched covering M — > C globalizes the so-called Levi-Civita transform 
Q of celestial mechanics. So it is natural to pose the question whether 

there exists a similar globalization of the Hopf map 

^ {z, az) 

used in the Kustaanheimo-Stiefel regularization of a 3-dim. Coulomb singularity. 
As already mentioned, this was done in pHS] for n = 2 centres. However, a 
generalization to arbitrary n seems to be unknown. 

We expect that the corresponding manifolds, i.e. four-dimensional analogs 
of Riemann surfaces, should have interesting topological properties. 

One last aspect concerns structural stability. Both in = 2 and three di- 
mensions the compact set Be C T^e of bounded orbits is hyperbolic and thus 
structurally stable. 

Thus if we continuously move the singularities by suitable maps M*^, 
u t— > .si{u), we obtain a family of Coulombic potentials Vu and correspondingly 
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a one-parameter family of bound states 6^; As the parameter varies cyclically, 
we obtain a permutation of the bounded orbits in Be = 

For d = 2 this action of the braid group on n strands of by permutations 
is nontrivial in general (see Remark 6.12 of [KK]). 

Although in (i = 3 dimensions the manifold A/C of non-collinear configurations 
is not simply connected for n > 3, the action of the fundamental group 7ri(A/C) 
on Be is trivial. 



B The Two-Centre Problem 

Here we shortly discuss the purely Coulombic two-centre problem, i.e. 

— Z\ — ^2 



\q-sx\ \q- S2\ 



-1 






W.l.o.g. we assume that the two centres are at si := and := ^ 

As is well-known (see e.g. Thirring |Th]), the problem is analytically integrable, 
see |GKM] for an application to satellite motion and [5R] for an application to 
semiclassics of the hydrogen molecule. 

The motion is integrated using the prolate ellipsoidal coordinates (^,r],(.p) E 
M+ X [0,7r) X [0,27r) with 



■ qi\ I cosh(C) cos(r;) 
= ( <?2 J = I sinh(5)sin{r;)cos(i/j) 
sinh(^) sin(r7) sin{ip) 
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As in these coordinates H is independent of (p, the conjugate momentum 

Py, = q2P3 - q3P2 

is a constant of the motion, equal to the first component of angular momentum. 
We set li := p^{xo) for initial conditions xq. 

Then by going to extended phase space and using a new time parameter s 
defined by 

£ = 2(cosh^(0-cos2(r7)), 
the new Hamiltonian function separates: 

dt , 

H:= —{H-E) = Hi + H2 
ds 
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with 

HM,0-=Pl + ViiO with l^i(0:=-i|^-2Z+cosh(0-2Ecosh2(0 



;2 

H2{Pv, V) ■= pI + V2{v) with V2{v) ■■= . 2/ N + 2^- cos(r/) + 2^003^(77), 

sm [r]) 

where Z± := Z2± Zi. 

The motion on H^^i^O) coincides — up to time parameterization — with the 
motion on H^^{E). Setting 

K:=H^{x^) = -H2{x^), 

we have three generally independent constants of the motion H, Hi and li, whose 
values are denoted by E,K and li, respectively. 

The bifurcation set is then given by the set of values for which the mapping 
from phase space to the constants of the motion is not locally trivial (see [AM| . 
Sect. 4.5). The most interesting subset is the one for li = 0, i.e. two-dimensional 
motion. 

By inspection of the extrema of the Vi one sees that for /i = the image of 
{H,Hi) is the region in M? bounded by the curves 

KAE)-=l 3 ,0>E>-Z^/2 
\ -2{Z+ + E) ,E< -Z+/2 

and 

K (E) ■= J 2i > 1^-1/2 

\ 2i\Z^\-E) ,E<\Z^\/2 ■ 

The bifurcation diagramme (see Fig. |2]) is the union of K^, and the lines 

E = , K = and Ko{E) := -2{Z+ + E) 

inside the image of (H,Hi). 

The line Kq corresponds to the {H, ifi)-values of the closed orbit wandering 
between the centres si and 5*2, i.e. having coordinate = 0. 

The following relations are useful for the scattering problem. In the (1-2)- 
plane, i.e. for = 0, we have the polar coordinates qi = rcos0, q2 = rsin0. 
Then = cosh^(,^) — sin^(?7) and tan(0) = tanh(,^) tan(r7) so that in the 
r 00 limit rj coincides with (p. 
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E 




Figure 2: Bifurcation diagramme 

In the same limit p,, = qip2 — q2Pi — e^^(pi sin?; + p2 cos(?])) coincides with 
the angular momentum qiP2 — q2Pi of the (l-2)-plane. 

This suffices to relate the asymptotic data {p^,r]^) :— limg^g±{pjj{s),r]{s)) 
with the ones in the original system (the times s"*" > s~ being defined by 
lim^^^ ^{s) = oo). The constant K is then given by (p^)"^ + V2{r]^). Equalling 
the elliptic integrals 

d'l 

y/-K - V2{V) 



then suffice in principle to calculate analytically the scattering transformation, 
but the expressions become rather lengthy. 



respectively 



oo 



s'^ — s =2 1 — — with 
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Note added in proof. 

The anonymous referee informed me about the interesting related article: 
S.V. Bolotin, P. Negrini: Regularization and topological entropy for the spatial 
n-center problem, which meanwhile appeared in Ergodic Theory and Dynamical 
Systems 21, 383-399 (2001). 

Concerning the context of the n-centre-problem on M^, the authors succeed to 
construct a global regularization of n attracting singularities, based on the local 
KS transform. Furthermore, they prove that for n > 3 the topological entropy is 
strictly positive for all energies E >0, whereas the present paper is only dealing 
with all energies above a positive threshold energy. Additionally, the authors 
prove several results for configuration manifolds different from R^. 

I do not believe, however, that these topological methods could be used to 
substantially simplify the proofs of the analytic results given in the present paper. 
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